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This  report  documents  and  discusses  twenty-three  programs — 
written  for  the  Hewlett-Packard  HP-67/97  programmable  calculators — 
covering  a wide  range  of  problems  of  interest  to  staff  officers  in 
all  the  military  services.  Using  this  material,  the  user  may  quickly 
obtain  answers  to  specific  questions  arising  in  meetings,  at  the  desk, 
or  in  the  field.  Full  documentation  is  given  to  clarify  the  back- 
ground of  a topic  and  to  enable  the  programming  of  a subject  of 
special  interest  for  a machine  other  than  the  HP-67,  but  with  com- 
parable power. 

In  general,  the  report  avoids  the  "slide-rule"  type  of  topic 
where  only  a given  formula  is  to  be  evaluated.  Rather,  topics  are 
chosen  that  would  consume  too  much  of  a staff  officer's  time  to  pro- 
gram because  the  underlying  mathematics  may  be  obscure,  because 
approximating  techniques  must  be  sought,  or  because  the  programming 
itself  presents  problems. 

Several  programs  reduce  published  volumes  of  tables  to  one  mag- 
netic card. 

The  major  part  of  this  research  was  supported  by  The  Rand  Cor- 
poration from  its  own  funds. 


SUMMARY 


The  report  is  summarized  by  an  overview  of  the  topics  covered 
in  it . 


Part  I.  Geographic  and  Orbital  Programs 

1.  Geographic  Coordinates  to  UTM  and  Conversely 

Army  tactical  maps  use  Universal  Transverse  Mercator  (UTM)  co- 
ordinates. For  joint  operations  with  the  Air  Force  and  the  Navy, 
coordinate  conversion  to  geographic  coordinates  and  the  converse  is 
essential.  Accuracy  of  this  program  is  better  than  10  meters  in  the 
northing  (distance  from  the  equator)  and  1 meter  in  the  easting  (dis- 
tance from  the  central  meridian  of  a zone). 

2.  Sunrise,  Sunset,  and  Twilight 

The  times  of  sunrise  and  of  the  various  categories  of  twilight 
are  important  in  planning  many  types  of  military  operations  and  ac- 
tivities, although  adverse  weather  conditions  all  too  often  vitiate 
such  planning.  This  program  gives  twilight  times  for  any  day  of  the 
year,  at  any  latitude  and  longitude,  and  at  any  altitude.  Accuracy  is 
three  minutes  or  less,  except  under  special  conditions  such  as  high 
latitudes . 

3.  Geodetic  Distances  and  Bearings 

The  usual  formulas  of  spherical  trigonometry  that  are  programmed 
to  give  great  circle  distances  and  bearings  employ  a spherical  earth 
of  some  mean  radius.  Distances  can  be  in  error  by  as  much  as  20 
kilometers.  The  program  here  uses  formulas  of  the  National  Geodetic 
Survey  based  on  Bessel's  solution  for  the  geodesic  on  an  ellipsoid  of 
revolution.  Accuracy  is  good,  about  0.1"  or  3 meters. 

4.  Reentry  Trajectories 

The  program  uses  Sec.  20  (Fourth-Order  Differential  Equations). 
For  a body  with  zero  lift  and  a given  "beta"  entering  the  upper  atmo- 
sphere, find  the  subsequent  range,  altitude,  and  velocity  to  impact. 
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li. . Satellite  Orbtlal  Elements 

This  program  solves  two  of  many  possible  orbital  problems:  Given 
a satellite's  Injection  altitude,  velocity,  and  flight  path  angle, 
find  the  remaining  six  orbital  elements;  or  given  the  Injection  alti- 
tude and  the  altitudes  of  perigee  and  apogee,  find  the  remaining  ele- 
ments. Equations  are  provided  so  that  otiier  problems  may  be  programmed, 

b.  Satell  1 1 e Track  1 n^j 

Given  the  time  and  longitude  of  equatorial  crossing  of  a satellite, 
select  a ground  station.  Determine  If  the  orbit  can  be  viewed  on  tliat 
pass,  and  If  so,  determine  Its  range,  bearing,  and  elevation  from  local 
horizon  to  horlzori  as  functions  of  time. 

Pa_r_^  1 1_. M i.LlJiiry  >Um1  els 

T_. The  Deer  HinU  j Defense  less  Bonilijers) 

The  model  assesses  the  expected  outcome  of  a time-limited  battle 
in  which  a group  of  armament-limited  interceptors  engages  a group  of 
defenseless  penetrating  bombers.  A deer  hunt  is  the  paradigm. 

8^, A Bomber  Penetration  Model  (Defended  Bombers) 

In  this  model,  the  bombers  are  not  defenseless.  As  part  of  mis- 
sion planning,  bombers  divide  their  payloads  between  defense  missiles 
and  ground  attack  munitions  to  maximize  weapons  delivered  to  ground 
targets . 

9.  Damage  Probabilities,  PVN  and  J^VN  Targets 

The  program  gives  damage  probabilities  for  nuclear  weapons  of 
given  yield  and  CKP  applied  against  PVN  and  QVN  targets  at  the  optimal 
alrburst  altitude. 

10.  Four  Deuces  (Precision  4.2-lnch  Mortar  JFjL^’) 

This  section  is  an  example  of  data-table  replacement  by  func- 
tional fitting.  It  applies  to  the  4.2-lnch  mortar,  reducing  firing 
table  corrections  and  meteorological  conditions  to  formulas.  The  pro- 
gram yields  corrected  shell  charge  and  corrected  azimuth  and  elevation 
for  precision  fire,  and  permits  a difference  in  altitude  between  mortar 
and  target. 
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11.  A Laser  Equation 

The  equation  programmed  applies  to  propagation  In  the  atmosphere 
and  allows  for  blockage,  thermal  blooming,  and  jitter  factors.  Given 
any  two  of  the  three  primary  variables  power,  range,  and  average  In- 
tensity at  the  target,  the  program  finds  the  third  factor. 

12.  Shaking  the  Dice  (A  War  Gaming  Example) 

This  section  provides  an  example  of  how  random  numbers  are  used 
in  a flrefight  model  to  assess  outcomes  quickly  in  war  gaming.  The 
example  employs  a conceptual  mortar  round  with  an  on-board  heat-seeker 
sensor  that  causes  the  round  to  home  on  an  armored  target. 

ill-  - Optimum  Allocation  of  Resources 

The  title  promises  too  much.  This  is  a topic  in  nonlinear,  con- 
vex programming.  Military  applications  arise  in  search  planning, 
allocating  weapons  to  target  classes,  and  allocating  budgets. 

Part  III.  Cost  Programs 

14.  Log-Linear  Cumulative  Average  and  Unit  Costing 

These  programs  implement  the  basic  assumption  of  learning  curve 
theory  as  it  applies  to  production.  That  is,  each  time  total  produc- 
tion doubles,  the  cost  per  item  reduces  to  a constant  percentage  of 
the  previous  cost. 

15.  Time-Phased  Procurement  Costing 

Consider  a system,  weapon  or  otherwise,  witli  several  major  com- 
ponents. Each  component  has  its  own  lead  time  and  its  own,  possibly 
segmented,  learning  curve.  Specify  a delivery  schedule  over  future 
years,  and  find  the  New  Obllgational  Authority  by  fiscal  year  to  sup- 
port the  program. 

16.  Cost/Benefit  Streams 

This  model  deals  with  the  decision  to  spend  money  now  as  opposed 
to  later  during  the  life  cycle  of  a weapon  system.  For  example,  should 
engineering  development  money  be  spent  now  in  the  expectation  that 
future  operating  and  support  costs  will  be  lower?  The  yardstick  is 
the  present  value  of  a discounted  stream  of  cost  and  benefits  (savings). 
An  "Internal  rate  of  return"  is  calculated  to  provide  go-no-go  for  the 
decision . 
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Part  IV.  Mathematical  Functions  and  AlRorlttnns 

17,  The  Normal  Function  and  Its  Inverse 

The  normal  function  (probability  integral)  is  pervasive  in  military 
calculations.  The  program  is  frequently  used  in  conjunction  with 
others,  such  as  that  for  the  Q function. 

18 .  Tlje  J^  Function  (Offset  Coverage  Function) 

The  Q function  is  used  in  radar  detection  tlieory  and  offset 
bombing  calculations,  as  well  as  in  calculations  of  collateral  damage 
to  point  targets. 

19^  Linear  Programm inj^_ and  3^3  M.atrix  Games 

Many  models  may  be  stripped  in  a meaningful  and  transparent  formu- 
lation to  tliree  activities  as  a programming  problem,  or  to  three  own 
courses  of  action  pitted  against  an  enemy's  three  courses  of  action, 
in  order  to  make  a command  decision  by  game  tlieory.  This  program  uses 
the  pivot  method  and  has  some  interesting  Indexing  aspects. 

20.  Fourth-Or^er  [)if  f erent  ial  Kquat  ions 

This  program  supports  applications  to  reentry  trajectory  determina- 
tion, Lanchester  models  of  combat,  and  optimal  control  tlieory. 


21 . Cur ve  Families  and  Mach  Numbers 

Military  data  are  frequently  presented  as  sets  of  tables  or  as 
families  of  curves,  with  a parameter  naming  the  family  member.  This 
section  suggests  methods  of  representing  these  data  through  curve- 
fitting, using  elementary  functions.  The  methods  are  applied  to  the 
determination  of  best  Mach  number  for  the  A-yo  aircraft  on  long-range, 
constant-alt itude  cruise. 

22.  Te n- P o i n t C.auss ian _ 1 n tejj^r a t^ o n 

This  is  a utility  program  for  evaluating  definite  integrals  as 
they  arise.  Accuracy  is  usually  excellent.  For  example,  incomplete 
elliptic  integrals  are  computed  to  eight  decimal  places  by  this  method. 

23.  Truth  Tables 


A calculus  of  propositions  is  tailored  for  ready  implementation 
by  the  calculator.  The  program  systematically  solves  problems  in 


symbolic  IorIc,  consisting  of  a set  of  logical  conditions  that  the 
atomistic  propositions  must  satisfy.  There  are  real-world  applica- 
tions, usually  overlooked. 
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If  your  t omperamtMit  Is  like  that  of  the  author,  this  description 
of  the  psychology  of  progranmilng  will  souiui  familiar:  After  tl)e 
usual  time-consuming  process  of  getting  tlie  mathematics  of  a topic 
In  sliape,  the  urge  Is  to  program  as  quickly  as  possible  and  make 
Independent  cltecks  of  tlie  validity  of  tl>e  outputs.  It  works!  And 
we  move  on  to  sometliing  else. 

Any  program,  liowever,  certai.nly  including  lliose  In  lliis  report, 
can  he  improved — can  be  shortened  and  made  more  eU’gaTit  and  t ranspan-nt  . 
I'lte  result  of  this  product- improvement  effort  may  he  to  ri-duce  the  run- 
ning lime  and  to  find  program  and  storage  space  to  extend  tlie  progr.im's 
capability.  A rt'cxam  inat  Ion  of  program  h^gii-  is  part  of  this  effort. 

The  program  mav  he  maiie  more  robust,  minimizing  opi'rati'r  errors  that 
occur  when  com|iIex  input  operations  are  iitherwise  rt'iiu  i ri'ii . 

If  you  as  a user  are  Interested  in  a particuliir  topic  in  this 
report,  you  may  choose  to  m.ike  this  extra  effort,  which  will  lie  re- 
paid with  an  enriched  understanding  of  haiui  calculator  programming. 

Final  ly,  you  are  invited  to  communicate  to  the  author  any  errors 
you  detect,  errors  and  unforeseen  restrictions  being  inevitalile  in 
a report  of  this  nature.  You  are  also  Invited  to  send  to  the  author, 
for  possible  future  prograiimilng , descriptions  of  topics  lliat  you  feel 
may  be  of  Interest  to  some  significant  subset  of  the  staff  officer 
community.  And  by  ail  means,  copli'S  of  your  own  programs  wouKi  be 
we  1 corned . 


t-  ■ 


yi-.i..  W-JI. ! 


-XV- 


CONTENTS 


PREFACE  lii 

SUMMARY  V 

ACKNOWLEDGMENTS  xi 

TO  THE  USER  OF  THIS  REPORT  xlii 

INTRODUCTION  1 


Part  I;  Geographic  and  Orbital  Programs 


Section 

1.  GEOGRAPHIC  COORDINATES  TO  UTM  AND  CONVERSELY  9 

2.  SUNRISE,  SUNSET,  AND  TWILIGHT  22 

3.  GEODETIC  DISTANCES  AND  BEARINGS  30 

4.  REENTRY  TRAJECTORIES  42 

5.  SATELLITE  ORBITAL  ELEMENTS  52 

6.  SATELLITE  TRACKING  60 

Part  II:  Military  Models 

7.  THE  DEER  HUNT  (DEFENSELESS  BOMBERS)  75 

8.  A BOMBER  PENETRATION  MODEL  (DEFENDED  BOMBERS)  84 

9.  DAMAGE  PROBABILITIES,  PVN  AND  QVN  TARGETS  89 

10.  FOUR  DEUCES  (PRECISION  4.2-INCH  MORTAR  FIRE  96 

11.  A LASER  EQUATION  108 

12.  SHAKING  THE  DICE  (A  WAR  GAMING  EXAMPLE)  116 

13.  OPTIMUM  ALLOCATION  OF  RESOURCES  122 

Part  III:  Cost  Programs 

14.  LOG-LINEAR  CUMULATIVE  AVERAGE  AND  UNIT  COSTING  133 

15.  TIME -PHASED  PROCUREMENT  COSTING  140 

16.  COST-BENEFIT  STREAMS  152 

Part  IV:  Mathematical  Functions  and  Algorithms 

17.  THE  NORMAL  FUNCTION  AND  ITS  INVERSE  159 

18.  THE  Q FUNCTION  (OFFSET  COVERAGE  FUNCTION)  163 

19.  LINEAR  PROGRAMMING  AND  3x3  MATRIX  GAMES  172 

20.  FOURTH-ORDER  DIFFERENTIAL  EQUATIONS  186 

21.  CURVE  FAMILIES  AND  MACH  NUMBERS  196 

22.  TEN-POINT  GAUSSIAN  INTEGRATION  203 

23.  TRUTH  TABLES  207 

APPENDIX;  97/67  KEY  CODE  CONVERSIONS  221 


-1- 


r 


^ 'V-L  'J  ■ JliUU 


INTRODUCTION 


"The  general  who  wins  a battle  makes 
many  calculations  in  his  temple  ere 
the  battle  is  fought.  The  general  who 

loses  a battle  makes  but  few  calcula-  j 

tions  beforehand. " 

- Sun  Tzii  Wu,  The  Art  of  War,  ca.  500  B.C. 

Programmable  hand  calculators  are  little  more  than  five  years 
old,  but  they  are  already  in  their  third  generation,  the  gestation 
period  being  one  and  a half  to  two  years.  Up  to  now,  they  are  unique 

in  our  inflationary  world,  in  that  each  new  generation  has  much  more  j 

power  than  its  predecessor,  but  sells  for  much  less.  This  cost  trend 
may  reverse  should  these  calculators  become  more  competitive  in  power 
with  microprocessors. 

We  can  assess  the  impact  in  the  civilian  sector  by  noting  that 

the  number  of  user  programs  submitted  to  the  Hewlett-Packard  HP-67 

program  library  is  approaching  3000.  The  PPC  (Personal  Programmers 
* 

Club),  with  more  than  2500  members,  is  a nonprofit  worldwide  group 
of  people  who  own  and  use  PPCs  (personal  programmable  calculators) . 

The  monthly  club  newsletter  contains  a wealth  of  programs  and  imagina- 
tive programming  techniques. 

Remembering  that  modern  digital  computers  were  initiated  by  the 
military  under  the  pressures  of  World  War  II,  it  is  curious  that 
these  PPCs,  these  powerful  little  animals,  are  not  as  equally  wide- 
spread in  the  service  of  the  Department  of  Defense  as  in  the  civilian 
sector. 

The  PPCs  are  used,  of  course.  The  Joint  Technical  Coordinating 
Group  for  Munitions  Effectiveness  (JTCC/ME)  under  the  JCS  has  had  25 
HP-67  programs  prepared  for  mission  planning  by  squadron  ordnance 
officers  in  the  Air  Force,  Navy,  and  Marine  Corps.  The  Strategic  Air 
Command  uses  the  HP-65  in  bombing  mission  planning.  Some  System 

•ft 
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Projoct  Offices  (SPOs),  such  as  the  K-lb  SPO  at  Edwards  Air  Force 
Base,  use  the  HP-b? . Junior  officers  are  using  tlieir  own  funds  to  pur- 
chase PPCs,  which  in  some  cases  must  represent  a tradeoff  against-  a new 
TV  set.  But  tliere  is  no  recognizable  community  of  users  in  the  mili- 
tary sector.  There  is  no  mechanism — no  clearinghouses  like  those  in 
the  civilian  world — to  exchange  programs,  to  share  ideas,  and  to  state 

requirements  for  new  programs.  The  notion  of  a loosely  organized  ' \ 

"national  security  users  group"  to  achieve  these  implied  objectives 
naturally  comes  to  mind.  We  hope  that  this  report  may  have  some 

catalytic  effect  in  accelerating  such  a development.  j 

The  hand  calculator  is  particularly  suited  for  military  use  be-  i 

cause  so  many  applications  can  be  made  in  the  field  or  In  a meeting  | 

where  a senior  officer  wants  a quick  answer  to  support  a decision,  or  | 

where  a briefer  is  to  be  confounded.  But  for  field  use  the  calculator  j 

as  ciivvcntlij  designed  would  probably  not  meet  military  specifications.  j 

The  operating  range  for  the  HP-f)7  is  H)°  to  AO^C  (50°  to  10A°F)  and  j 

the  battery  pack  life  under  continuous  use  is  about  three  hours  before  ! 

recharging  or  replacement  is  required.  However,  current  machines  are  | 

j 

compatible  with  avionics,  producing  little  or  no  interference  with  j 

sensitive  electronic  circuits.  s 

But  powerful  as  they  are  in  their  domain,  the  PPCs  are  far  from  i 

a final  answer  to  personal  computing,  although  this  statement  depends 

on  their  future  evolution.  In  preparing  this  report,  many  instances  i 

occurred  where  much  more  storage  than  available  was  needed  and  where  i 

it  was  frustrating  not  to  have  available  a programming  capability  of  j 

more  lines  of  code  with  a higher-level  interpretive  language.  ; 

i 

Again,  the  civilian  sector  is  leading  the  way.  More  than  120  I 

companies  are  now  manufacturing  microprocessors  with  peripherals  for  j 

home  use,  and  more  than  900  home  computer  dealers  in  the  United  States  ) 

are  marketing  these  machines  at  relatively  modest  prices.  Memory  may  j 

be  added,  there  is  keyboard  input  and  cathode  ray  tube  display,  with  ^ 

_? 

BASIC  apparently  the  language  of  choice.  The  military  is  lagging,  ' 

even  though  it  is  a reasonable  bet  that  many  staff  officers  would  like  j 

to  be  freed  from  the  computing-center  bureaucracy  in  doing  their  daily  | 

jobs.  j 


J 
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But  once  more,  too  strong  n position  should  not  be  taken.  Mill-  1 

tary  computing  In  general  requires  large  main  fr;imes  to  support  ^ 

extremely  large  data  bases  and  progr.jms  with  a million  or  more  lines 
of  code.  One  would  certainly  hesitate  to  try  to  use  a microcomputer 
for  logistic  management  or  for  solving  three-dimensional  partial  dif- 
ferential equations. 

Nevertheless,  there  Is  a real  gap  In  the  spectrum  of  required 
computing  capability  to  meet  military  requirements,  a gap  whose  filling  ; 

this  report  can  only  adumbrate. 

A word  of  apology  is  in  order.  Recorded  program  cards  are  not  j 

1 

provided  with  this  report.  The  reasons  are:  i 
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• No  recipient  Is  likely  to  use  all  programs; 

• The  per-copy  cost  of  the  report  would  be  high; 

• It  requires  10  to  20  minutes  to  key  in  a program  and  check 
it ; and 

• Hopefully,  the  keyer  will  understand  the  program  and  be 
able  to  modify  or  tailor  It  to  his  or  her  desires. 

It  Is  recommended  that  users  step  through  the  Illustrative  prob- 
lems to  get  the  mechanics  straight.  And  It  is  always  a good  idea  to 
do  a problem  twice.  Krrors  In  keying  are  easy  to  make,  especially 
when  under  pressure. 

Finally,  what  Is  to  bo  said  to  the  staff  officer  who  wants  to  pro- 
gram his  or  her  own  problems  on  a PPC?  The  natural  question  for  the 
officer  to  ask  first  Is;  What  bounds  a problem  that  can  he  "fitted" 
to  the  machine? 

The  general  answer  Is:  If  the  problem  can  be  formulated  as  a 
c'htin  of  subproblems,  each  of  which  Is  within  the  machine's  coding 
and  storage  capability,  then  there  Is  In  principle  no  bound.  For 

A 

example.  In  the  prediction  of  tides  by  harmonic  an.ilysis,  17  con- 
stituents (cosine  terms)  each  with  three  constants  are  employed.  Since 
these  terms  need  only  be  added,  one  program  card  and  five  data  cards. 
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used  successively,  would  suffice.  As  other  examples,  six  linear 
algebraic  equations  in  six  unknowns  can  be  solved  using  both  sides 
of  two  cards,  and  a Star  Trek  battle  can  be  programmed  with  eight 
cards. 

For  problems  that  can  be  chained,  the  practical  limitation  is 
execution  time,  which  can  be  long  and  hardly  acceptable  if  many  prob- 
lems are  to  be  run,  as  in  tidal  prediction. 

But  not  all  problems  can  be  chained.  Operations  with  matrices 
of  order  higher  than  five,  and  solutions  of  partial  differential  equa- 
tions, are  usually  noncha inable . 

Even  if  a problem  s}ionld  fit,  it  is  frequently  liard  to  see  how 
to  make  it  actually  conform  to  the  calculator's  Procrustean  bed.  This 
could  be  because  tlie  underlying  matliemat ics , Including  approximating 
techniques,  is  beyond  one's  reacli.  The  lielp  of  a specialist  colleague 
is  then  essential.  Once  tliis  mathematical  hurdle  is  cleared,  program- 
ming— which  is  really  an  art  form  with  personal  brush  strokes — can  be 
exasperating.  Advice?  Read  and  understand  good  programs,  as  many  as 
possible — something  that  few  of  us  have  the  self-discipline  to  do. 

As  a postscript  to  this  Introduction,  an  as  yet  unexplolted  area 
of  the  military  application  of  PPCs  should  be  mentioned.  Two  or  more 
people  may  operate  their  calculators  in  parallel,  engaging  in  a co- 
operative, Interactive  exercise. 

For  example,  two  submarines  may  be  allies  in  a simulated  battle 
against  one  enemy  boat.  The  purpose  of  the  exercise  is  to  examine, 
by  repeated  simulation  runs,  the  tactical  utility  of  communications 
between  the  two  friendly  boats  during  the  battle — ranging  from  none, 
through  restricted,  to  complete  information  and  command  exchanges. 

Each  player  has  his  own  program  which,  by  sampling  from  probability 
distributions,  shows  the  output  of  his  sensor  systems  in  respect  to 
target  position  and  bearing,  and  the  damage,  if  any,  inflicted  by 
ordnance  launched.  Each  player  keeps  his  own  log  and  battle  plot. 

At  each  battle  increment  (say,  I")  minutes  of  real  time),  the  calcu- 
lators may  be  physically  exciianged  so  that,  as  appropriate,  Infornui- 
tion  can  be  entered  In  assigned  storage  registers,  and  the  calculators 
then  returned  to  the  right  boats. 
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As  another  example,  a War  College  seminar  may  be  examining  the 
cost  Implications  over  the  next  ten  years  or  more  of  various  possible 
strategic  postures.  Weapon  systems  may  be  phased  In  and  out.  New 
systems  require  research  and  development  monies  and  time.  In  general, 
each  weapon  system  has  cost  profiles  of  funds  required  for  RUT&K, 
procurement,  and  annual  maintenance  and  operating  expenses.  The  cost 
envelope  of  each  weapons  system  with  respect  to  time  is  calculated  by 
the  seminar  member  assigned  that  system.  All  programs  are  the  same, 
differing  only  in  their  cost  and  time  parameters.  The  seminar  leader 
totals  the  year-by-year  costs  of  all  systems  In  the  posture  and  checks 
for  feasibility  against  an  assumed  yearly  ceiling.  After  discussion, 
the  seminar  members  revise  phasing  or  numbers  procured  and  go  through 
another  iteration  to  see  if  the  celling  is  reached  or  exceeded,  and 
to  determine  if  the  posture  Is  balanced  in  regard  to  the  threat  and 
required  missions. 

These  examples  have  indeed  been  programmed  for  interactive  com- 
puting on  large  computers;  but  this  is  time-consuming  and  facilities 
may  not  be  re.uiilv  available.  Tlie  suggested  use  I’f  Tl’t^s  in  parallel 
is  an  option  that  c.ni  be  implementt'd  quickly  and  can  provide  a shake- 
down for  more  sopii  ist  ic.i  t ed  approaches,  whicli  in  some  cases  may  prove 
not  to  he  warr.inted. 
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L.  GEOGRAPHIC  COORDINATES  TO  UTM  AND  CONVERSELY 

1.1.  REFERENCES 

a.  Universal  Transverse  Mercator  Grid,  AMS  Technical  Manual 
No.  19,  Army  Map  Service,  Corps  of  Engineers,  Washington, 
D.C.,  1952. 

b.  Map  Projections , P.  Richards  and  R.  K.  Adler,  North-Holland, 
1972. 

c.  Map  Reading,  FM  21-26,  Department  of  the  Army,  October  1960. 

1.2.  DISCUSSION 

Tactical-scale  (1:50000)  Army  maps  use  the  Universal  Transverse 
Mercator  Grid  (see  Ref.  c).  The  map  borders  show  latitude  and  longi- 
tude ticks,  but  it  is  difficult  to  locate  the  geographic  coordinates 
of  a point  with  any  precision.  Conversely,  Air  Force  maps  use  geo- 
graphic coordinates  only.  Consequently,  in  joint  operations  such  as 
targeting,  coordinate  conversion  from  one  system  to  the  other  is 
essential.  FORTRAN  programs  exist.  The  program  used  at  The  Rand 
Corporation  has  132  lines  of  code.  Although  perfectly  adapted  to  the 
preparation  of  coordinates  for  a list  of  agreed  targets,  it  hardly 
meets  the  requirements  of  ad  hoc  field  use. 

The  UTM  system  covers  the  world  between  80°S  and  84°N.  Starting 
at  the  180°  meridian  of  longitude  and  moving  eastward,  the  globe  is 
divided  in  zones  6°  of  longitude  in  width,  numbered  1 to  60.  Each 
zone  has  a central  meridian  (CM) . The  following  formulas  relate  zone 
number  (ZN)  to  the  CM: 


ZN  = (CM  + 183)/6 
CM  = (6  • ZN)  - 183 

For  example.  Fort  Knox,  Kentucky  is  about  86°W.  Hence  the  ZN  is  the 
rounded  value  of  (180  - 86)/6,  ZN  = 16,  and  CM  = -87  or  87°W.  (See 
Ref.  c for  further  details  on  lettering  8°  zones  south  to  north,  and 
on  double-lettering  for  100,000  meter  squares  within  each  6°  x 8°  block.) 
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The  value  asslsned  to  the  CM  in  each  zone  is  501), DOl)  meters, 
called  tlie  j\iliSC  Hence  locations  in  a zone  west  of  the  CM 

Itave  an  less  titan  500,000  and  conversely.  The  is 

the  distance  from  the  equator  in  meters.  For  the  Southern  Hemisphere, 
the  equator  is  assigned  a ru't'thtfhj  of  10,000,000  meters  and 

numbers  decrease  southward. 

Tl>e  major  complication  in  coordinate  conversion  is  that  allowance 
must  be  made  for  the  earth's  oblateness.  Hence  the  equatorial  radius 
a and  the  polar  radius  b must  be  selected.  Actually,  a and  the  re- 
ciprocal of  the  flattening  f = (a  - b)/a  are  given.  For  the  Inter- 
national Spheroid, 


a = b 378  388  m , 1/f  = 297  . 

Since  f = 1 - ^l  - g",  where  e is  the  eccentricity, 

= 0.006  722  67  . 

Onfortunately , different  spheroids  (different  a and  f)  are  used 

for  different  areas  of  the  \corld,  for  historical  reasons.  For  example, 

the  Clarke  1866  spheroid  is  used  for  North  America.  The  other  spheroids 

used  are  Clarke  1880,  Kverest,  and  Bessel.  The  International  Spheroid 

2 

is  used  for  Kurope.  (Consult  Ref.  a.)  Consequently,  the  data  a,  1:“, 
n = (a  - b)/(;»  + b)  used  here  have  tc*  he  changed  for  certain  parts  of 
the  world. 

The  full  equations  for  the  conversions  (Refs,  a and  b)  are  quite 
lengthy  because  extreme  accuracy  is  desired  in  surveying  applications. 
For  military  purposes,  it  is  possible  to  dock  the  tails  of  these  formu- 
las and  still  get  accuracies  better  th.an  I meter  in  the  inat'  (K') 

and  better  than  10  meters  in  the  (N)  — the  distance  from  the 

equator  in  meters. 


T. 


iJPPii'JiTJ 


-11- 

l . L ■ C'-eo^rapl tic  Coord lj\a to  UTM  Grid  Coordinates 


(1)  + (ii)p‘  + (in)p^ 

(1) 

(iv)p  + (V)p'^  a 0,  K - 500  000  i E'  . 

(2) 

South  of  the  equator, 

N - 10  000  000  - N . (3) 

i 

The  slven  coordinates  are  latitude  (|)  and  longitude  \,  Then 

p = O.OOOl  • AA,  where  AA  Is  the  difference  of  longitude  from  the  CM, 

# 

in  t'iW'tuit'-.  K'  is  the  (positive)  distance  from  the  CM. 

(1)  » S • kjj,  whore  (A) 

S = Ai(>-Bsln2t}>  + CslnAi|’  • 

S Is  the  true  meridional  distance  from  the  equator  In  meters  and 
A “ all  - n + 0.7S  n'd  - n)  ) 


B • 1 . S a n ( 1 - n ) 

C “ 0.^1375  an“(l  - n)  , n « (a  - b)/(a  + b) 


kp  “ 0.9996,  the  central  scale  factor  to  reduce 
distort  Ion. 


k^^  sln‘  1"  • 10^ 

(11)  • -■  • u sin  2 vj*.  where 


(5) 


u 


2 2 

E*"  sin*"  is  tlie  radius  of  curvature  in 
the  prime  vertical. 
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(III) 


ky  siu^  1"  • lo‘^’ 


V ) 2 

• V)  sin  cos  (Jit'S  - ton  (Ji) 


k,  sin  I"  • lO**  • u cos  (Ji 


k^l  sin^  1"  • 10* 


U cos  ij'(2  cos'"  (Ji  - I ) 


L’vLvn.'_  . .I'T.M  O.r  IJ.  Coordinates  to  Ceojj^rapli  ic  C^oordinates 


ij)  “ <J>'  - 1 (Vll)  - (VllI)  q^l/JhOO 


AA  = 1(IX)  q - (X)  q 1/1600 


q - K'  • 10  i 0 


iji"  = N/Ak,, 


N/k^j  + B sin  2 «Ji"  - C sin  4 >Ji" 


u - a/  \ I - e"  sin'  (|i" 


(Vll) 


10  , 

2k^  sin  1"  U' 


(VI  ID 


. ijii’A  + A -i/yiA'A'A. 


24k|^  sin  l" 


ky  sin  i"  U cos  (J)' 


U)  . 1 -f  2 tan' 

6k|^  sin  1"  cos  (Ji' 


1.5.1  USER  INSTRUCTIONS 


1.1  GEOGRAPHIC  COORDS  TO  UTM  GRID 


INSTRUCTIONS 


INPUT 

OATA/UMTS 


KEY  LAT  (D.MS)  STO  D 

KEY  LONG  (D. MS)  STO  E 
(-  FOR  W.  LONG.) 


KEY  CM 


PRESS  A 


STO  7 


OUTPUT  IS  NORTHING  (M)  IN  D 


5 PRESS  R/S 


OUTPUT  IS  EASTING  (M)  IN  E 


INT.  GEOID  DATA  CARD 


° ;37S788.0M  8 

B.B0t?2267  I 
ko  e.99969$U  3 

lO^sin  1 " 4.84813888  4 

7888.888888  5 

508888.8888  6 

A 8787845.458  8 

B 18188.99888  B 
C 18.97888888  C 


ISTOl  1 D 1 

c „j  Lr_i 

[stoILe  J 

(1“  1LIZ1 


I ~ ir.TJ 

L \lz:a 

i:  jri:.] 

[ ~ ]i  I 

r;  jl  :"j 


L :i[  .:“i 
j LZJ 
I " 

L““iLrT.i 

rr_j  i:r;j 
r_iLi:] 
r~:]  rr.i 

I JIZ'I 
I „ ii„  I 


OUTPUT 

DATA/UNITS 


456820 
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1.6.1  GEOGRAPHIC  TO  UTM  COORDINATES 


STtP  M'i  { 


DEC.  DECS. 


A\  IN  SECS 


01  ! 

H,1S‘ 

16  J6 

Pi: 

STOP 

Ji  H 

oh: 

1 

»4 

or 

' 

-75 

Pis 

SIN 

41 

01? 

sar 

76  17 

oso 

-75 

0.'; 

Kcir 

Jt  14 

iv: 

% 

HC 

o:.' 

-75 

o:: 

SIN 

41 

O.T 

RCLS 

76  10 

% ^ • 
i « >. 

\ 

-75 

^ 

- 

-45 

'J? 

PfiP 

76  14 

o}? 

P'R 

16  45 

O.T 

KflO 

St  11 

001 

> 

-75 

* 

-55 

0-'.' 

SCI  7 

76  07 

y - » 

\ 

-75 

oj: 

sro? 

75  0? 

O’: 

FCLC 

76  14 

O': 

snv 

75  00 

OJS 

SIN 

41 

07 : 

Y1 

57 

07C 

m; 

76  01  ' 

Oil 

\ 

-75' 

oi: 

1 

01 

017 

- 

-45 

oil 

CHS 

01.1 

IS 

54 

o;. 

RCIO 

76  00 

01 : 

+ 

-74 

OIS 

; N 

57 

01.' 

^ros 

75  00 

pro 

FCLC 

76  07 

pji 

* 

07' 

or: 

-75 

PJ.7 

SIN 

41 

or; 

sets 

76  00 

or: 

-75 

1 

04 

EAT.  IN  DEC  DECS 


(1)  ' llDp' 


RADIANS 


NORTHING 


EAT.  D I LONG.  V,  E 


SkUP  HfcV  INIRV  kIVCOOfc  COMMtNIS  kev^NTHY  KEY  CODC  COMMENTS 


1.5.2  UTM  USER  INSTRUCTIONS 
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IBIS  PAGE  IS  BEST  QUALITY  PRAGtlCAlCil 
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1.6.2  UTM  TO  GEOGRAPHIC  COORDINATES 


STEP 

KEY  ENTRY 

KEY  CODE 

COMMENTS 

STEP 

KEY  ENTRY 

KEY  CODE 

COMMENTS 

001 

Ml 

tLBLR 

21  11. 

^sr 

NU 

003 

0CL6 

36  06- 

'058 

m 

43" 

,003 

KCL£ 

36  IS. 

E 

059 

X 

-33 

,004 

-45. 

060 

[088 

RCL8 

X M 

,00! 

MS 

16  31. 

881 

4 

04 

,006 

E£X 

-23. 

082 

y» 

3L 

.003 

6 

06. 

083 

+ 

-24 

.00$ 

-24. 

884 

RCL4 

38  *4" 

,0t3 

STOl 

33  4<_ 

q 

885 

+ 

-24 

010 

.010 

RCLD 

38  14. 

888 

RCL3 

38  08~ 

.011 

KL0 

36  11. 

087 

4 

•4" 

.013 

+ 

-24. 

088 

r» 

31' 

[0i3 

KL3 

36  03 

08? 

4 

-24" 

^814 

+ 

-24. 

070 

070 

2 

•2" 

015 

R«0 

16  46 

DEGS. 

071 

4 

04" 

,01S 

ST09 

35  09 

r (11) 

072 

4 

-24 

,81? 

4 

04 

073 

EEX 

-23" 

818 

X 

-35 

074 

3 

03" 

[818 

SIN 

41 

075 

0 

00" 

0^ 

030 

RCLC 

36  13 

~0?S 

X 

-35" 

031 

X 

-35 

077 

RCLI 

38  48' 

822 

CHS 

-22 

078 

4 

04" 

033 

RCL9 

36  09 

Q79 

Y” 

31" 

034 

3 

03 

080 

080 

X 

-35 

03! 

X 

-35 

081 

ST02 

35  02" 

(Sm)q^ 

036 

SIN 

41 

082 

RCL8 

38  08" 

027 

RCL8 

38  12 

"083 

TAN 

43" 

038 

X 

-35 

084 

RCL8 

38  08" 

039 

♦ 

-55 

"085 

X* 

53“ 

030 

030 

RCLD 

38  14 

088 

•r 

-24" 

031 

RCL3 

36  03 

"087 

RCL4 

38  04" 

033 

4 

-24 

"088 

■f 

-24" 

033 

♦ 

-55 

"08S 

RCL3 

38  03" 

034 

RCLR 

36  11 

090 

080 

X* 

53' 

03! 

4 

-24 

081 

•r 

-24 

036 

R*D 

18  48 

DECS 

082 

2 

02 

037 

ST09 

35  09 

i>'  (12) 

"883 

•r 

-24" 

038 

SIN 

41 

■084 

EEX 

-23“ 

039 

N* 

53 

085 

1 

01" 

D40 

040 

RCLl 

38  81 

08o 

8 

08 

041 

X 

-35 

887 

X 

-35 

043 

1 

01 

■088 

RCLl 

38  48“ 

q 

043 

- 

-45 

088 

X» 

53“ 

044 

CHS 

-33 

too 

100 

X 

-35“ 

04! 

fX 

54 

'101 

RCL3 

38  02 

046 

RCL0 

38  00 

"102 

- 

-45" 

047 

4 

-24 

'103 

CHS 

-23~ 

848 

l/x 

32 

"104 

RCL5 

38  05“ 

049 

ST08 

33  86 

I'  (13) 

"105 

,+ 

-24 

OSO 

050 

RCL9 

38  09 

108 

RCL8 

38  09 

051 

TAN 

43 

107 

♦ 

-55" 

852 

X* 

33 

108 

*HHS 

18  35“ 

D.MS 

853 

3 

83. 

108 

STOD 

35  14“ 

854 

X 

-35. 

110 

110 

R/S 

51 

DISPLAY 

055 

5 

85. 

111 

RCL9 

38  0* 

IISL. 

•f 

112 

-JtL 

iL. 

REGISTERS 


° a 

' e2 

'(2III)q^ 

^ k 
^0 

‘'l0®5in  1 ■' 

" 3600 

® 500000 

^ CM 

6 

1/ 

SO 

SI 

S2 

S3 

S4 

S5 

S6 

S7 

S8 

S<> 

* A 

® B 

O 

n 

" E,X 

' q 

2.  SIINKISK,  SIINSKT,.  IWllJCA'l 


2.1. 

a.  Russell,  nugau,  .uul  Sti'wart,  •li!,' c. ''Ii ’n;/,  lHiin  v**  Co.,  Ni’w  York, 

H45. 

b.  'Plu'  .ViJkm'ocj/  Mr'kinu.'  for  tiu'  Yoor 

U.S.  Covoriimont  I’rlntiug  iM  f loo,  W.isliinutoii,  D.C.,  l‘>7(i. 

i‘.  Fj.'pl.vuiti'i>n  Fupi'i I'rnont  to  tlw  Art  >'o>:<'n;i' o.i ! F]'>u''’!ot'i  o.  Her 
Majesty's  .Stationery  Off  lee,  l.ondon,  l^bl. 

- .'I'AAAAA'AA 

Charts  are  prepared  and  issued  for  e.ii'h  maj»'r  military  opi'r.it  Ion 
or  operational  area  jtlving  sunlight,  moonlight,  and  lid.il  dal.i. 

Nautie.al  twilight  (sun's  zenith  angle  from  102“  to  9h”)  provides  enough 
illumination  for  most  types  of  ground  aelivlly,  although  bomb  Iv'adlng 
and  repair  work  require  artlfieial  light.  Civil  twilight  (sun's 
zenith  angle  from  ‘lb®  to  90”)  pi-rmits  normal  day  activities  such  as 
observed  artillery  fire  and  visual  bombing.  Sunrise  occurs  when  t hi' 
sun's  upper  limb  h.as  a zenith  angle  of  90”.  This  makes  the  zi'uith 
distance  of  the  sun's  center  90”'j0'  (90.8  3”),  allowing  14'  for  hori- 
zontal refraction  and  lb'  for  the  sun's  semid i amet er . For  some  air- 
craft applications,  a correction  of  l'.17»*li  is  .added,  where  h is  the 
altitude  in  feet.  If  H is  in  kllofeet  and  decimal  degrees  are  used, 
this  correction  Is  0.bl7/H  degrees. 

The  fundamental  relation  (Ref.  c,  p.  401)  is 

cos  h • - tan  iji  tan  6 + sec  i{>  sec  iS  cos  z , (1) 

where  h and  6 are  the  hour  angle  and  declination  of  the  sun  at  the 
time  of  the  phenomenon,  il*  Is  latltvide,  and  z is  the  zenith  angle. 

(The  correction  of  the  declination  from  ephemerls  noon  to  approximate 
' Islng  or  setting  Is  at  most  O.l",  a refinement  we  will  neglect.) 
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Neglectlng  \1  - e'^  (=  0.99986),  and  solving  for  0 j'lelds 


sin  9 = 


sin  a 


lie  cos  a ’ 


(4) 


where  the  + branch  is  used  from  the  vernal  equinox,  21  March  (Day  89), 
to  the  autumnal  equinox,  23  September  (Day  275). 

From  (2)  a first  approximation  to  a is  simply  = 2nD/365.  This 
is  refined  by  using  a = + A in  (2)  to  get  the  correction 

sin  a 

A = 77; (radians)  . (5) 

60  - cos  a 

o , 

Finally,  using  direction  cosines,  the  declination  becomes 

sin  6 = - cos  8 sin  U , (6) 

where  U = 23.44°  is  the  inclination  of  the  earth's  axis.  The  sign 
becomes  + between  the  two  equinoxes  as  the  earth  passes  through  the 
summer  solstice,  21  June. 

Next  a formula  for  the  Equation  of  Time  (EOT)  is  required.  The 
EOT  is  the  difference  in  hour  angle  of  the  sun  and  the  fictitious  mean 
sun  used  for  ordinary  time.  The  difference  owes  to  two  causes;  (1) 
the  variable  motion  of  the  sun  because  of  the  eccentricity  of  the 
earth's  orbit,  and  (2)  the  obliquity  of  the  ecliptic. 

The  figure  on  p.  147  of  Ref.  a suggests  that 

EOT  = - A sin  (0  - a)  - B sin  (20  + b)  . (7) 

We  get  from  that  figure  the  approximate  values  A = 8,  a = 5.92°, 

B = 10,  b = 4.73°.  Here  Day  0 (D  = 0)  is  25  December,  since  the  EOT 
is  0 on  that  date. 

The  extrema  of  the  EOT  are: 


i 
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12 

February 

1)  = 49 

KOT  = -14.29 

min 

14 

May 

D = 140 

EOT  = 

+ 3.72 

min 

2b 

duly 

1)  = 213 

EOT  =■ 

-b . 4b 

min 

3 

November 

1)  = 313 

KOT  = +lb.41 

mi  n 

Roplaoo  A by  A + AA,  etc.,  substitute  in  (7),  take  only  first-order 
terms,  and  four  linear  equations  in  the  unknowns  AA,  Aa,  AB,  Ab. 

These  are  solved  quickly  by  Program  7 of  tlie  Hewlett-Packard  Math 
Pac  i.  The  corrected  values  of  the  parameters  are 


A = 7.4A47  , a = 5.935  , B = 9.894  , 


4.941 


Tlie  resulting  mean  absolute  error  throughout  the  year  with  respect 
to  the  tabulated  values  of  the  HOT  is  24  sec. 

Rising  and  setting  times  are  now  computed  by 


Rising  = 12  - HOT  - h 

in 

Setting  = Rising  + 2h  . 

These  are  local  me.in  times  witli  respect  to  the  central  meridian  (t^N) 
of  a given  time  zone.  To  correct  for  other  longitudes,  subtract  4 
min  for  each  degree  east  of  the  CM,  since  the  sun  is  earlier,  and 
add  4 min  for  each  degree  west  of  the  CM.  The  correction  is 
progr.'immed . 


(1)  The  day  number  I)  (1)  = 1)  + 3)  for  a given  date  is  needed, 
counting  from  Christmas  as  Pay  0.  Subtract  I from  the  month  number 
and  multiply  by  30.42,  the  average  number  of  d.ays  in  ;i  month.  Take 
the  Integral  part.  Kor  month  numbers  1,  8,  9,  10,  11,  12,  add  b; 
for  months  3,  4,  5,  b,  7,  aild  5;  and  for  month  2,  add  7.  Pinisli  by 
adding  the  days  of  the  date. 

(2)  Kor  the  longitude  correction  (f  l.Bh  0)  add  ibO"  if  Wi'st 
longitude  (entered  neg.itive).  Then  obtain  the  correct  central 
meridian  by  checking  whether  the  fractii'nal  part  I'f  the  longitude 
divided  by  15  is  less  than  or  greater  th.an  0.5  (1/2  hr). 


I 

I 

) 


[i 
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( ))  Al  tiatos  wiu'n  llu'  sun's  cior  1 inat  Inn  is  o U>so  ti'  0,  t’»>rnuila 
(•i)  I'an  yiold  a numht'r  vury  sli(;litly  ^roator  In  absoluto  vahu*  t l\an  1. 
This  wauld  nonoratu  an  Krror  si>;nat.  Siii'li  a niiniluT  is  roplai-oil  l>y  1 
in  r l.BI,  S. 

(4)  Tliis  program  riHiniroit  220  stops.  Tln'ro  was  lu't  sp.ioo  to 
program  rising  and  sotting  in  tiu'  Soutticrn  liomisptu'ro.  To  do  tliis, 
prooood  as  folli’ws: 

• Kind  tlio  sun's  doclination  on  tiu‘  di'siroil  dato. 

• t".l\ai\go  tlto  sign  and  uso  I'ormula  (h)  to  gt't  a now  0 and  a 
now  dav  numhor  ( ibS  • 0/ 5b())  , dit'l'oring  hv  about  b nitmtlis. 

• Kind  t lio  rising  t imo  lor  t ho  lattor  d.iy  in  t ho  Northorn 
lli'misphori'. 

• Add  K.OT.,  - K0T|. 

K.x.unpl_o.  Kiuii  sunriso  I'U  S May  at  38".^  (.Contral  Mi-ridian).  Run  pro- 
gram witli  +i8'’.  Red.  2 to  got  iS  “ IS. 48“  (tho  t ruo  doolination  is 
\b.lS“).  Rl'l.  I ti'  got  KOTj  - .ObOO  l)r . By  (b)  with  6 ■=  - 1 S . 48 . 

0 ■ 47. d or  )12.1.  Cliooso  tho  lattor.  0 = 117.  or  7 Novombor.  Rvin 
program  with  dato  11.07.  Rising  t imo  is  b h 32  m (b.Sii,  KOT , = i).2b. 
Rising  t imo  is  b.S3  + 0.2b  - 0.0b  = b.73  or  b h 43  m.  Tho  valiu'  for 
this  oxamplt'  givon  on  p.  Sbb  I'f  Rof.  (b)  is  b h 4S  m. 

(5)  I’rogr.tm  running  t imo  is  .ibout  2S  soo . Compared  with  tho 
values  t.tbulati'd  in  Rof.  b (pp.  414lf),  orrors  aro  0 to  3 min  with  (' 
and  I tho  most  likolv  v.iluos.  Tlu'  orrors  can  bo  groator.  howovor, 
at  high  l.ititudos  at  dates  oloso  to  tlioso  with  twiliglu  Lasting  all 
night  . 


Kx_amj>lo.  Ast  rommi  ioa  I twi  1 ight  , S0"N,  14  .Inly:  Tlio  oi'inputod  v.iluo 
is  0 h 31  m versus  tho  aotual  0 32  for  beginning,  and  23  3d  versus 
23  32  for  end  of  twilight. 
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2.5  USER  INSTRUCTIONS 


2.  SUNRISE,  SUNSET  AND  TWILIGHT 
ZENITH  A,  DATE  B,  LAT.  C,  LONG.  D PRESS  A,  R/S 


STEP 

INSTRUCTIONS 

INPUT 

OATA/UNITS 

KEYS 

OUTPUT 

DATA/UNITS 

1 

LOAD  BOTH  SIDES  OF  PRGM  AND  DATA  CARDS 

1 11 

1 . 

1 M 

I 

2 

ZENITH  DESIRED  STO  A 

1 II 

I 

SUNRISE  (UPPER  LIMB)  90.83 

90.83 

1 IIENTI 

90.83 

CIVIL  TWIUGHT  96 

1 II 

I 

NAUTICAL  TWILIGHT  102 

1 M 

I 

ASTRONOMICAL  TWILIGHT  108 

1 11 

I 

1 11 

I 

TO  CORRECT  FOR  ALTITUDE  IN  KFT 

1 II 

I 

ADD  0.617  V H TO  THE  ABOVE  (40  KFT) 

3.90 

1 1 1 

I 

94.73 

isToii  A 

I 

94.73 

4 

DATE  AS  MM.DD  STO  B 

1 M 

I 

(N.B.  OCT  9 IS  10.09) 

10.09 

ISTOll  B 

I 

10.09 

1 II 

I 

5 

latitude  (NORTH  ONLY) 

1 M 

I 

DD.N-...'  STO  C 36°  22*  N 

36.22 

1 STOl  1 C 

I 

36.22 

1 II 

I 

6 

LONGITUDE  + DDD.MM  STO  D 

1 II 

I 

(-  FOR  W.  LONG)  121°  8‘W 

-121.08 

iSTOll  D 

I 

-121.08 

1 11 

I 

7 

PRESS  A TO  GET  RISING  HH.MM  (5h  45m) 

1 II 

I 

5.45 

1 II 

I 

8 

PRESS  R/S  TO  GET  SETTING  HH.MM  (I7h  57m) 

1 II 

I 

17,57 

1 11 

I 

9 

DAY  NUMBER  OF  DATE 

1 RCL 1 1 0 

I 

291 

RCL  0 AND  SUBTRACT  9 

9 

1 11" 

I 

282 

1 II 

I 

10 

EQUATION  OF  TIME 

1 II 

I 

DSP  4,  RCL  1,  g H.MS 

1 RCL  1 1 1 

I 

0.2128 

TO  GET  .MMSS  (12m46s) 

1 9 IH.MSI 

.1246 

1 11 

I 

11 

DECLINATION  OF  SUN 

1 RCL 11  2 

I 

-5.9663 

RCL  2,  g H.MS  to  get  + DD.MMSS  (-5° 58') 

1 5 1 H.MS 

-5.5759 

1 II 

I 

12 

RISING  IS  STORED  IN  3 

1 RCLl  1 3 

I 

5.7632 

SEHING  IS  STORED  IN  4 

1 9 IH.mS 

5.4548 

(g  H.MS) 

1 11 

I 

1 II 

I 

13 

ERROR  MEANS  TWILIGHT  LASTS  ALL  NIGHT 

1 11 

I 

iL 
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' I-.  JlLSllANCKS  ANJ)  ,B1;AR  1 N^S 

t 

) . I . . 

I a.  Krwin  Schmid,  T)‘ta>uiuLit!0>:  i-on  Without 

TabuLlttomi,  unpublished  Ms.,  National  Ccodetic  Survey, 
duly  1972. 

■ b.  r.  A.  Smitli  and  H.  C’..  Massey,  A AAbb  ft-')'  r/.'i’  b'odiU'!^' 

, hn\'yije  Cor!p:ituiu,'ti,  The  Rand  Corporation,  l’-4930,  January 

j 1973. 

t 

i 

A-A-. . 

The  programs  usually  written  for  great  circle  distances  and  bear- 

A 

* ings  .issume  a spherical  earth  of  some  mean  r.ulius  and  employ  the  ele- 

t 

I mentary  formulas  of  spherical  trigonometry.  The  results  can  be  in 

‘ error  bv  as  much  as  20  kilometers.  I’o  tiu'  geodesist  (whicl>  1 am  not), 

, suit)  progr.ims  are  to  be  anathematized.  . Hut  beyond  the  evident  demands 

of  surveying,  there  ia*t’  applications  in  the  military  and  international 
' domains  where  much  greater  precision  is  required. 

^ In  1828,  F.  W.  Bessel  gave  the  general  solution  for  the  geodesic 

on  an  ellipsoid  of  revolution.  Two  differential  equations  must  be 
solved  as  part  of  Bessel's  rigorous  procedure. 

Unfortunately  for  rigor,  the  geodesist's  w'orld  is  neither  perfect 
nor  static.  Periodically,  the  semi-major  and  semi-minor  axes  of  the 
geoid  are  chaoge4-sjinu'yhat  in  value.  We  recommend  using  WCS  72  (.World 
Ceodetic  System  1972),  now  rather  geiTeriaJ  ] v^ accept  ed  , which  assumes 


a = h 378  135.0  m b = b 35b  750.233  m , 

and  hence  a squared  eccentricity  of 

= 1 - (b/a)‘  = 0.00b  b9A  407  . 

(The  reciprocal  of  the  flattening  f is  298.25b.) 

The  fotVAird  problem  in  geodesy  takes  the  latitude  and  longitude 
of  a station,  an  azimuthal  bearing  measured  clockwise  from  north,  and 


a il  ist  aiii'i'  or  ni'inli'i  lo  sokiiumiI  , aiul  asks  lor  tli»‘  K^'>*Ki'apli  i ca  I ooonll 
natos  ol'  tiu'  tiM'inlnus  ol  llu>  sf^nu'nt  . 'I’lu'  i ntU’ftU'  iirohlt'in  wants  t In- 
m'OvU't  U'  illstaiu't'  bftwoiMi  two  sla(i(<ns,  p.lvi'u  tliolr  >;fOKra|>li  lea  I oo- 
oiallnatt's,  as  w»‘ll  as  l ho  l)i*arlnKa  ol’  I'aoh  station  I rom  I lu“  otiior. 

ktillAr.U'NS 

Tho  is|uat  Ions  lor  I ho  solution  ol  tlu'so  I wv'  prohloms  aro  ratlu'r 
lon>tlliv.  To  oonsi>rvt>  spaoo  In  this  roport  , thov  ar*-  not  roproilnooil 
lu't'i'.  Thov  aro  ti>nnil  li\  Rol  . h,  aval  lahli*  I rom  Tho  Kami  (Tirpi'i'al  Ion 
PuhlloatLons  Dopartmont  . Santa  Mi>nloa,  C^  4040h.  I'or  thoso  roailors 
who  tlolvo  Into  this  rotoronoo,  not«-  that  In  I ht>  solution  pro 

^’ranmioil  horo  thort'  Is  a roplaoomont  In  lv>rmnla  (H)  ol 

* * - hv  I - k . 

I + k‘/A 

J -6 

sliioo  k'/A  is  li'ss  than  I ''  10  , ami  thofi'  is  a ilolollon  In  tormnla 

(10)  ol  tho  torm 


1 

k • oos  I)  S'  ‘sin  1 AS'  , 


which  is  ol  tin'  orilor  1 v 10 

This  Is  ilono  to  savo  soari't*  program  spai'o  ami  to  aooo  li-rat  t*  oxi" 

out  Ion  sli>thtlv-  In  ooinparlnn  tho  rosnlts  with  tho  moro  I'xaot  valnoj 

louml  hv  tho  National  Oooilot  lo  Survoy,  tho  rosultlns  latltmh'  inav  ho 

/ 

In  orror  up  to  0.2"  ((>  iiu'tors),  and  tho  louKitiulo  by  oi'ns  lilt>rabl  v 
loss  than  this, 

1 . •'* . PKlUIRAM  Ji‘\'>ll'i'‘LA*''.'’AWARl)  SOI.IIT  ION) 

(1)  Tln-ro  aro  throo  plaoos  lit  tho  prv>>iram  whoro  oari'  nuist  h«' 
taki'n  ti>  onsuro  that  tho  arotannont  niv«'s  an  an^li'  In  tho  oorroot 
<piadrant . Tho  'n  tan  *'  I'nnot  Ion  vnt  tho  lU'-h7  prodnoos  att^los  onlv 
In  quadrants  I and  IV,  but  tin'  roi'l  annu  lar-t  o-po  1 ar-oi>ord  I nat  o o»'n- 
vorslon  providos  tho  oorroot  qviadrant  I'or  +/  + , +/-,  and  -/■¥.  K» 


Price  to  private  Individuals:  $i.00  postpaid, 


in  the  y value,  the  numerator  witli  i'fa  sln«.  I’ress  KNTKR.  K»‘v  in 
the  X value,  the  iiei\i>minator  with  its  slun.  I’ress  -*•  I’.  Press 
h X y to  net  the  angle. 

(2)  Stai'k  manipulation  is  used  in  two  places  to  hold  valxies  in 
stack  storage  until  they  can  he  placed  in  prlm.iry  storage,  avi>iding 
the  extra  steps  Involving  't  I’  •»-*■  .S'  if  secondary  storage  were  us»*d. 
It  is  good  practice  in  such  progr.xmming  1 1'  use  'SST'  in  run  mode  and 
get  successive  t rai'es  of  the  stack  C4>ntents  hv  'g  STK ' , (This  is 
/■e.’-hugglng  rather  t Ixan  <;’i  -bugg  ing . ) 

(1)  llecause  of  program  size  (218  steps),  tlie  user  is  asked  to 
make,  if  necessary,  a simple  final  correct  it'n  1 1'  the  illsplayeil  liutgi- 
tiule  (8  .tnd  under  l.S.l,  User  lust  riu-t  lints)  . It  is  impintant  to 
use  ' IPO,  iTiS,  h 11.NS+'  rather  th.in  st  r.i  iglit  f itrwaril  suht  r.ii't  iint  in 
step  8. 

Kxan^tlos.  Tlie  fol  hawing  comparisons  with  National  (lei'det  ic  Survey 
values  use  their  get'ld  constants  I'f 

a = h )78  I'tS.OO  m h = h ISh  7‘j'').7h  m 

= 0.00(1  54 ■)  . 


Moretiver,  their  i ni'ryiie}  solutions  (geoid  ilistances  as  si'gnunit  lengths) 
are  used  to  chock  solutions.  (Note:  I'o  i-inwert  kilonit'tiMs 

to  nautical  miles,  divide  by  1.852  exactly.) 

I.et  the  station  of  origin  bo  tlte  munlclp.-il  airpi'rt  .ti  I'airbanks, 
Alask.i— tiA"40'08.05",  - 1 47"  5 1' 5 1 . lb"  . Tlte  n we  find 


Place 
Az imut  h 
K.i  li'ineters 
l.at  1 1 tide 
A lat 
Long  i t utie 


Los  Angeles 
I 15"I2’ 18.42" 
3‘)(i  1 . b444 
14"02' 5‘).t)  I" 
0.07" 
-I  18"I  r5‘».‘»(i" 


Jakarta 
28l"25'57.0b" 
-I  I I44.  I480 
-ii"  10 '00.  lb" 
0.  lb" 
I0b"4‘)”i<).y  I" 


Te 1 Av i V 
r>7"45'2b.  12" 
-0250.2287 
12"04' 50. 7')" 
0.21 " 


I4"4b'0000" 
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If  tiu'  a/.imutit  (boarlan)  Is  ^roalor  tliaa  180”,  tlu*  dlslaiu-o  Is  oiUorinl 
as  a minus  nuanfity.  Ttu'  A's  sl\i'w  il\o  il  I f f ori'iiou  witli  ri>sp»>i'l  ti’  tin' 
NlIS's  muro  accurate  value  (0.1(1"  in  lailtiule  Is  ab»>ut  10  feet). 

K 5.  _ rROl'^R/\M  NpTKS  X,1NV1':^RSK  SOl.U'riON) 

(1)  In  tills  pro^;ram  alse,’  tlii'  2-p;iraniet «*r  arct.'in>;ent  proceilure 
must  he  folluwetl  to  net  .mules  In  the  j-errecl  ipi.ulrant. 

(2)  Ag.'iln  bec.iuse  of  progr.'im  sp.lce,  we  negU'i-t  terms  involving 

) 2 
k and  except  tor  one  cas»'  those  involving  k . 

())  On  tlu-  first  K/.S  if  the  number  .ipiie.ir  i ng  (the  dilfi-rence 
In  r.'idians  of  the  longitudes  of  the  second  .iiul  I irst  stations)  is 
gre;it»’r  th.in  ti,  then  k»'v  Oll.S,  hn,  2,  x,  t,  R/S.  Otlu-rwlse  we  get 
tlu-  greater  of  the  two  geodesic  dist.inci-s  lu-twei-n  the  stations.  Of 
cours<>,  this  is  rea^lily  pr«'gr.immed  but  we  do  not  havi-  remaining  .ivail- 
.ihle  sp.ice. 

We  now  re\>eat  the  aluwe  i-xamplcs  using  the  invi-rse  solutii'u.  ^ 


Place 

l.os  Angt'les 

.lak.iiT  a 

Te  1 Aviv 

l.at  it  tide 

lA^or 

-h” 10' 

12"0S’ 

l.ong  i tilde 

-ll8'’l/«' 

10h‘’')0' 

14  "4h’ 

K 1 lometers 

I'lnl  .(v'lOti 

1 1 144.  1418 

82S8.222I 

A (meters) 

1.8 

h.2 

h.h 

A/,  limit  li 

ri'S'’12'  18. b2" 

78"  14'0  1.00" 

02"  14’  U.hS" 

A A/. 

0.20" 

O.Oh" 

0.0  1" 

ll.ick  A/. 

ri8".Vi'02  .Oh" 

l')'S‘’07’  18.00" 

i 78"'i2’  18.  1 1" 

A 11 /AX. 

0.08" 

0.01” 

0.01" 

in  this 

program,  .i/.  Imutlis  ar*- 

less  than  180"  and 

give  the  angle 

from  N to  the  eastl>ound  port  Ivui  of  the  g,eodesic. 

The  accur.ii’V  should  be  .iccept.ible  to  .ill  but  (he  mo^t  deni.inding 
user.  We  pav  for  this  accur.-u'V  in  the  usual  I'oin  t't  liiero.ised  com- 
puting time.  Tile  re.ider  mav  wish  ti'  progr.-im  using  a me.in  e.irth  r.idius 
.iiul  the  formul.is  of  splierical  t r Igonomet  rv , as  nu-nt  ioned  In  .Sec.  1.2, 
to  i>ersu.ide  liimselt. 
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BiHMuse  of  lack  of  program  spac'o,  tlie  usor  must  tako  aicount  ol 
tl\roo  proi'.oJuros: 

• If  hotli  given  latitudes  are  negative  {Si>uthern  Hemisphere), 
ehange  signs  of  both  latitudes  and  I'oth  longitudes. 

• If  both  stations  t»ave  tlie  same  longitude  (meridional  .ire), 
the  program  will  show  an  error  at  step  044  (1/tan  0).  In 
this  ease  add  0.01'  to  one  longitude  (0.01)  anil  run  from  the 
beginning,  although  the  solution  is  unstable  for  sueh  small 
differenees  in  longitude. 

• Registers  A,  R,  0,  0 are  used  both  for  initial  and  inter- 
mediate storage.  If  a lU'w  iiroblem  is  to  be  run  with  the 
same  first  station,  its  eoordinales  must  be  restored.  You 
ean,  however,  store  the  eoord I n.i t es  also  in  SO  and  S7.  being 
eareful  to  preeede  and  follow  by  f 1’  S. 


t 

? * 
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3.6.1  PROGRAM  USTING  ^ 


STEP 


t20 

130 

140 

150 

' 

160 

11' 

PCff 

38  03  ' 

^ISf 

f 

-24 

1I-! 

* 

-55  ■ 

1*: 

17C 

CHS 

-22 

« ♦ ■■ 

3T0!> 

35  99 

Act 

! 7! 

PCLS 

36  09 

1 1 : 

PfLg 

36  96 

Iff 

SIN 

41 

ir 

♦ 

-55 

HUB 

^I’f 

RC1.3 

36  03 

Ilf 

EItTt 

-21  ■ 

STACK 

r- 

2 

02 

11.- 

ENTt 

-21  ■ 

MANIPULATION 

|jSB 

Iff 

' 

-35 

♦ ? ' 

FCL€ 

36  96 

1 1,'S 

COS 

42 

!•  * 

* 

-55  ■ 

4 t 

-35 

Iff 

fTCf 

35  03 

2 a 

Iff 

r 

-55 

4 

y^v 

-41 

IB 

lf.« 

RCL8 

36  08 

li- 

STOe 

35  06 

°2 

ISC 

2 

02 

Iff 

5 

02  ■ 

♦ r ‘ 

T 

-24 

Imi 

f’^Ff 

35-24  03  ■ 

a 

isf 

I 

01 

'■CLt 

36  06  ' 

ISf 

+ 

-55 

' 

SIN 

41  ■ 

1 1 ; 4 

PCI  3 

36  09 

13? 

Ef-Tt 

-21  ■ 

TO  GET  X^,  o-j 

4 -7C 

\ 

-35 

Iff 

RCLE 

.56  tft> 

IN  SAME  QUAD 

Icc 

•f 

-55 

IJl 

fOS 

42  ■ 

— 

* *7 

FC1.5 

36  85 

Iff 

PCLf 

36  05  ' 

lES 

COS 

42 

' 7' 

:as 

42  ■ 

ISS 

\ 

-35 

13-‘ 

-35  ■ 

190 

ISO 

FCLS 

36  08 

Iff 

fP 

34 

!.«: 

X 

-35 

A v’  » 

.v?v 

-41 

i.’f 

4 

04 

4 

EWIT 

-21 

17? 

f 

-24 

Iff 

"CLf 

36  07  ■ 

1.°-! 

36  00 

13- 

- 

-45  ■ 

l.’f 

<#2 

53 

i-if 

STOE 

35  15  ■ 

AX 

l?i 

X 

-35 

u: 

k:y 

-41 

ir 

1 

01 

u: 

ST07 

35  07 

i?f 

FCLl 

36  01 

I 

RCLE 

36  06  ' 

1.=.’ 

- 

-45 

l-S' 

fOS 

42  1 

200 

'>fi  *» 

FC1.5 

36  05 

1-T 

PCL5 

36  05  ' 

fc: 

COS 

42 

Nf 

SIN 

41  ■ 

7^0 

A 

-35 

i4r 

-35 

PCLS 

36  09 

Nf 

flN--' 

16  41 

SI 

X 

-35 

14.- 

fTC4 

35  04  ■ 

7 fr 

- 

-45 

Iff 

’’NN 

43  ■ 

fCf 

FCLE 

36  15 

4 r 1 

PfLl 

36  01 

mu 

fC' 

♦ 

-55 

Iff 

f 

-24  1 

mm 

fc: 

FCl-E 

36  12 

1 

TUN-' 

16  43 

fCS 

HIS* 

16  36 

154 

R-D 

16  46 

2 

210 

fif 

C-R 

16  45 

Iff 

■•HMS 

16  35 

fi; 

+ 

-55 

1ft 

P S 

51  ■ 

DSP  LAT 

2*2 

■Si 

16  46 

* 

FfLS 

36  99  1 

fif 

*Hns 

16  35 

If: 

“1 

02  1 

fl4 

FTN 

24 

15.’ 

-35 

fif 

• LSLe 

21  00 

Iff 

!IN 

41  ] 

fi: 

.”:y 

-41 

I:'.* 

Pfff 

36  03  1 

21  ^ 

- 

-45 

Iff 

4 

04 

fif 

FTN 

24 

* 

y 

-35 

lf4 

CPS 

42  " 

2?0 

Itf 

A 

-35  j 

Ic":* 

PCL8 

36  08 

io7 

,\ 

-35 

1.<S 

4 

04. 

COMMENTS 


A 

B 

C 

D 

E 

0 

a 

b 

c 

d 

e 

1 

0 

1 

2 

3 

4 

2 

1 

7 

e 

, ■! 

9 

3 

DSP  LONG 


SET  STATUS 


FLAGS 


TRIG 


(MSP 


ON  OFF 

0 D □ 

1 □ □ 

2 D □ 

3 D n 


OEG  □ 
GRAD  □ 
RAD  □ 


FIX  a 
SCI  n 
ENG  □ 


3.2  GEODETIC  DISTANCE  BETWEEN  TWO  STATIONS, 

AND  THEIR  BEARINGS  Z 

(INVERSE  SOLUTION) 


INSTRUCTIONS 


lAT.  OF  1ST  STATION,  D.MS  STO  A 


- IF  S.  LAT 


LONG.  OF  1ST  STATION,  D.M.S.  STO  B 


- IF  W.  LONG. 


INPUT 

DATA/UMTS 


64.4909  1 (STO II  A ] 

dlJCZJ 

[zz][z:j 

-147.5151  I LSTOl  ( B I 


LAL OF  2ND  STATION,  STO  C 

LONGr0^2ND  STATION,  STO  D 

jw/ss^~ 


ON  R/S,  IF  > TT  (3.1416), 


KEY  CH^  Htt,  2,  x,  +,  R/S 

"OKT NEXT  R/S,  SEE  DISTANCE  (KMS.) 


ON  NEXT  R/S,  SEE  BEARING  STATION 


1 TO  STATION  2 ( < 180,  FROM  N TO 


INBOUND  PORTION  OF  GEODESIC)  D.MS 

ofT lact~r7s,  see  bearing  from  statT^ 

2 TO  STATION  1 (AS  ABOVE)  D.MS 


-118.1359 


FOR  S HEMISPHERE  SEE  NOTES 


DATA  CARD 


< STO  4> 


0.081  819  356 


_VTr«(!_  STO  1 0.996  647  176 


STO  2 6 356.750  23 


(WGS  72 


OUTPUT 

DATA/UMTS 


3961 .64 


135.1218 


158.4502 


[ 


THIS  PAGTE  IS  BEST  QUALlTT  m«lCAteI 
nOM  OOlPY  rURMiiiilit)  10  , 


S iP  KEV  entry  key  code 


3.6.2  PROGRAM  USTING 

COMMENTS  STEP  KEY  ENTRY  KEY  CODE 


LABELS 


SET  STATUS 

FLAGS  TRIG 

(MSP 

DEQ 

T1 

FIX 

n 

QRAD  LI 

SCI 

u 

RAD 

tn 

ENG 

tj 

n 

.... 
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I 

i 

j 

I 

\ 

I 

j 

1 

t 
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4.  REKNTRY  TRyUlXTORIKS 


4a.  ^ .REFERENi^S 

a.  M.  M.  Moo,  "An  Approximation  to  the  Re-Entry  Tra  jei'torv." 

ARi}  Aoutnial,  January  19b0,  pp.  S0-5J. 

b.  R.  Blum,  "Re-Entry  TraJ  eot  I'l"  ies  : Elat  Earth  Approximation," 

ARt'  Jou}‘>u.}L,  April  19b2,  pp.  blb-b20. 

hIREUSSION 

For  a body  witli  zero  lift  and  s'ven  we  ighl -to-draj;  ratio  (beta), 
entering;  the  upper  atmosphere  at  time  0 with  assigned  altitude,  veloeity, 
and  patli  an>;le,  find  at  desired  time  intervals  tlie  subsequent  ranne , 
altitude,  velocity,  and  path  angle  to  impact. 

An  "exact"  solution  (an  analytic  solution  is  not  possible)  re- 
quires the  numerical  integration  of  two  second-order  differential 
equations,  botli  highly  nonlinear.  The  tecluiique  is  to  merge  a program 
for  tlie  functions  of  the  differential  equations  witli  a slightlv  re- 
tailored version  of  Program  20,  wliich  solves  fourth-order  differential 
equations,  also  providing  on  a separate  card  an  initialization  program 
to  determine  values  for  the  variables  at  the  first  time  interval,  t = 1>. 
Tills  merging,  retailoring,  and  initialization  is  a useful  model  for 
similar  applications  in  other  work,  wlien  a particular  set  of  equations 
is  to  be  solved  frequently. 

The  equations  adopted  are  those  for  a nonrotating  round  earth. 

There  are  three  assumptions: 

1.  In  the  we iglit-to-drag  ratio  B = mg/t'j^A,  the  drag  coefficient 
is  held  constant,  altliough  it  actually  varies  witli  Mach 
number ; 

2.  Surface  gravity  g is  used  uncorrected  for  altitude  bv  (r^/r)"”, 

where  r is  the  radius  of  the  eartli  aiiii  r is  the  distance  to 
e 

the  body  from  the  earth's  center; 

3.  riie  density  of  the  atmospliere  is  approximated  bv 

p(y)  ^ 0.002  )7  exp  {-y/24,000}  slug/ft  \ where  y = r - r^^ 
is  the  altitude  in  feet  (Ref.  a). 


t 

) 

i 
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Wt'  luivo  at  Rand  an  on-lino,  t inu--shar  in>;  program  in  .UiSS,  wiitton 
hv  0.  Ki'pluirt  , whloli  port'i'nns  t ho  nuiiiorli-al  int  t'nrat  ion  n.sin>;  t ho 
t I'ui  t It-ordor  Rnnjio-Rutta  t’l'ruuilas  and  a t iiiu'  intorVvil  ot  I so^- . Himu  i- 
a ov'nvoniont  haso  of  oomparison  with  t ho  simp  lor  and  ooarsi'r  proil  iit  I'r 
ooMi'oti'r  approaoh  v'f  I'ro^ram  .’I'  is  avMilahlo. 


I'ho  hasio  variahlos  aro: 


r = distanoo  from  oarth's  oont«'r  1. 1 t 1 

= polar  an>;lo  from  initial  vt'otor  to  onrii'nt  vootor 
V = spood  v'f  t ho  hoilv  U t/st'i') 

0 = an^lo  moasurod  posit ivoly  downward  t rom  t ho  looal 
luH'izontal  to  tho  volooitv  vootor. 


In  thoso  variahlos,  tho  oipuit  ions  of  motion  .ir»' 


1 d 
r dt 


0 


in 


d~r 

) 

dt" 


>;  + s i n 0 , 
m 


• f-") 


whi'ro  tho  dr.ig  foroo  dno  to  .lir  is  I'j^  » p t v )V"l'^^A/ J . Wo  .ilso  h.ivo 
I'idv^i/dtl  = V oos  0 , dr/ilt  = -V  sin  0 

in 

> > > 

V“  = fdr/dt^“  + (r  d^'/dt)-  . 


1 

i 


i 

i 


ILsiiif;  tho  v.iri.ihios 


\ « r vf  (rani-o)  v = r - r faltitudo) 

t'  1' 


11 
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l 111'  svstom  ot'  tour  1' i rst -ordor  li  i t t'oroiit  i.i  1 oqu.Uiiius  is 

X ' = u V ' = V 


V + r jP 

L ■ o 


(4) 


r 

o 


SP(v)V 

S - -Vg-v 


wlioro 


■1 

V 


i,l  + v/r  )“  [.r  + v"  . 
o 


As  a matter  ot  interest,  it  will  lu>  found  by  eompul  ing  dV^^/dt 
and  dO^^/dt  that  for  some  reentry  altitudes  V ran  initi.illv 

and  0 ean  fpiieh-up).  An  example  will  illustrate  these 

phenomena . 

Turn  now  to  tlie  determination  of  Xj,  Vj,  U|,  at  t = h,  the 
first  time  inerement,  whieh  are  values  needed  for  the  method  iM'  I’ro- 
gram  20.  I’his  will  he  programmed  to  initialize  the  numerieal  inte- 
gration. We  h.ive 


Xi  = u^^h  + u'  ir/2  V,  = + v^,h  + v’  ir/2 


i') 


Ui  = Uj,  + u;^  h 


Vi  = + v’  h . 


where  h is  the  interval  of  integration  (10  see  seems  suitable). 

The  equations  for  the  flat-earth  .ipproximat  ion  .ire  found  from 
(4)  by  putting  r^^  = The  interested  re.uier  may  wish  to  t rv  pro- 
gr. miming  this  e.ise  along  the  lines  of  this  seetion  to  eonip.ire  the 
answers  with  those  for  a round  earth. 

4.4.  JIKOCRvW  NOTKS 

Program  20  is  modified  as  follows: 
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1.  At  the  beginning;  ot  l.Bl.  A,  4 is  stored  in  registers  H and 
K,  replacing  and  u^j  used  during  tlie  initialization. 

2.  Most  of  LBl.  4 is  deleted,  since  v.c  are  working  with  4 
equations  and  we  want  to  compute  velocity  and  patli  angle 
for  display. 

3.  The  subroutine  programming  of  u'  and  v'  is  st  raiglit  forward , 
but  we  have  had  to  use  program  steps  for  g/2  .ind  g because 
storage  registers  are  not  available. 

i 

j^x^nitj^le.  Reentry  conditions  at  t = 0 are 

= 250,000  ft  = 10,000  ft /sec  B = 1000  6^^  = 3°  . 

I'he  tabulation  below  shows  v. lines  at  intervals  of  50  sec,  although 
tlie  step  size  used  was  10.  I'he  numbers  in  parentheses  are  the  outputs 
ol  Kephart's  more  refined  integration  procedure.  His  density  function 
was  used  , 

0.0027  exp  {-y/ 2 3, 500}  , 


instead  of 


0.00237  exp  {-y/24,000}  , 

whose  constants  are  stored  on  our  data  card  (in  S8  and  S^). 

After  the  "knee"  of  the  trajectory  the  altitudes  are  in  error 
by  about  1000  ft  low.  The  reason  is  that  10  sec  is  too  great  an  inter- 
val in  this  region.  If  b.  = 2 sec  is  used,  the  maximum  error  in  altitud 
reduces  to  about  80  ft  and  the  iither  values  are  st'iisiblv  exai-t  . .And 
if  li  = 1 sec,  the  value  used  in  Kephart's  calculation,  tlie  maximum 
difference  in  altitude  is  3.5  ft. 

This  is  entirely  of  theoretical  interest  as  a comparison  of  two 
metliods  of  numerical  integration:  the  fourth-order  Runge-Kutta  and 
the  extended  Hamming  predictor-corrector  metliod  emploved  here.  In 
the  practical  sense,  tlie  variation  of  tlie  drag  coefficient  with  M.ich 
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TIME 

(sec) 

X 

RANGE 
(n  mi) 

y 

ALTITUDE 

(ft) 

V 

VELOCITY 
(f  t/sec) 

0 

ANGLE 

(deg) 

20 

97.33 

(97.37) 

199  751 
(199  753) 

29  993 
(29  987) 

4.60 

(4.58) 

50 

242.03 

(242.07) 

132  807 
(133  228) 

28  692 
(28  730) 

4.05 

(4.00) 

80 

363.79 

(364.40) 

80  636 
(81  490) 

19  009 
(19  176) 

4.27 

(4.19) 

110 

422.90 

(424.67) 

46  342 
(47  349) 

6 657 
(6  786) 

8.26 

(8.06) 

140 

441.84 

(444.03) 

19  889 
(20  860) 

2 109 
(2  218) 

23.19 

(22.36) 

160 

446.25 

(448.64) 

3 474 
(4  328) 

1 154 
(1  207) 

43.38 

(41.82) 

170 

447.37 
(— ) 

-4  333 
(— ) 

921 
(— ) 

55.10 
(— ) 

number,  the  departure  of  the  atmosphere  from  the  ideal  one  used  in 
the  calculation  and,  above  all,  the  probably  asymmetric  ablation  of 
the  nose  cone  can  vitiate  such  accuracy.  Taking  h = 10  sec  is  a 
reasonable  choice  that  minimizes  computing  time  while  yielding  accept- 
able accuracy  for  the  purposes  to  which  this  program  should  be  put. 


4.5  USER  INSTRUCTIONS 


4.  RE-ENTRY  TRAJECTORIES 


INSTRUCTKmS 

INPUT 

DATAAIMTS 

KEYS 

OUTPUT 

DATA/UMTS 

LOAD  DATA  CARD  (BOTH  SIDES) 

1 II  1 

0.00 

1 ir  1 

TIME  INTERVAL  h,  STO  D,  STO  1 

10 

1 STOl  1 

0 1 

1 STO)  1 

1 1 

/3  STO  B,  V STO  D,  9^  STO  E 

1000 

1 STO|  1 

B 1 

1000.00 

(9  IS  POSITIVE  DOWNWARDS) 

30000 

ISTOl  1 

D 1 

30000.00 

5 

1 STOl  1 

E 1 

5.00 

Vo  STO  6,  STO  7 

250000 

1 STOl  1 

6 1 

250000 

1 STO)  1 

7 1 

250000 

LOAD  INITIALIZATION  CARD 

1 II 

1 

250000 

1 II 

1 

PRESS  E y, 

1 II 

E 1 

224358.95 

1 11 

1 

LOAD  PROGRAM  CARD 

( (( 

1 

1 II 

1 

PRESS  A 

1 II 

A 1 

1 II 

1 

SEE  t (h  PAUSE) 

1 M 

1 

20.00 

1 II 

1 

SEE  RANGE  n.mi.  (f-x-) 

1 II 

1 

57.32 

1 II 

1 

SEE  ALT  ft  (f-x-) 

1 II 

1 

199752.44 

1 II 

1 

SEE  VEL  ft/SEC  (f-x-) 

1 II 

1 

29989.44 

1 II 

1 

SEE  ANGLE  (f-x-) 

1 II 

1 

4.60 

1 II 

1 

IF  MORE  TIME  NEEDED  TO  RECORD.  USE  R/S 

1 II 

1 

30.00 

1 II 

1 

(NOTE:  VALUES  DO  NOT  AGREE 

1 II 

1 

146.00 

1 11 

1 

WITH  EXAMPLE  IN  TEXT  BECAUSE 

1 11 

1 

176225.76 

1 II 

1 

A DIFFERENT  DENSITY  FUNCTION 

1 11 

1 

29870.78 

1 II 

1 

IS  USED  HERE.) 

1 11 

1 

4.41 

1 11 

1 

1 

1 II 

1 
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4.6  PROGRAM  LISTING 


STIP  Ker  ENTfiV  KEV  CODE  COMMENTS  STEP  KtV  ENTRY  KEV  CODE  COMMENTS 
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4.6  RE-ENTRY  TRAJECTORY  PROGRA^\  (ROUND  EARTH) 


•TER 

Ktv  imthy 

REV  COOE 

COMMENTS 

STEP 

KEV  ENTRY 

KEY  COOE 

COMMENTS 

saM 

ti-BLn 

DSP' 

4 

SJ02 

sroE 

j 

STOI 

11 

-t:  92 

04 

J3  14 
sr  15 
03 
JS  4e 

(FROM  NOTE  1) 

•5.’ 

*50 

959 

869 

961 

962 

963 

DSZl 

BCLi 

4 

0S3I 

RCLi 

4 

X 

18  35  4t 
36  45 
-55 
It  35  46 
36  45 
04 
-35 

wmm 

feiMta 

OM 

iroc 

.V  13 

864 

865 
966 
867 
96S 
869 
979 
871 

972 

973 

974 

975 

976 

977 
878 

er? 

4 

5 

T 

ISZI 

ISZI 

ISZI 

STOi 

RCLC 

1 

STOD 

X-9? 

6704 

ISZI 

ISZI 

6T01 

-55 
05 
-34 
it  26  46 
16  26  46 
16  26  46 

35  45 

36  14 
01 

-45 
35  14 
lt-43 
33  04 
16  26  46 
16  26  46 
33  01 

M!) 

tie 

til 

tl2 

tl3 

tl4 

tl5 

tit 

ttr 

tie 

tie 

tse 

tzi 

t22 

t23 

«Lei.c 

6S6i 

1521 

SICi 

Rcie 

k 

t 

D$2t 

DS2I 

ecu 

4 

is:: 

RCLi 

DS2] 

21  13 
23  45  , 
It  26  46 

35  45 
00 

-35  ■ 
03 
-35 
le  25  46 
It  25  46 

36  45 
-55 

it  26  46 
36  45 
It  35  4t 

(SEE  FROM  20) 

010 

070 

■■■ 

MiriMI 

DtO 

HHB 

e24 

STDi 

35  45 

oao 

080 

SL6L4 

21  04 

(FROM  20  STEFS) 

t2S 

K4 

-31 

081 

RCL£ 

36  IS 

077 

t26 

1S2I 

16  26  46 

982 

STOD 

35  14 

t27 

isz: 

io  fc't  'tO 

983 

6SB2 

33  03 

t28 

ISZl 

ib  43b 

984 

PCLl 

3b  bi 

080 

829 

STOi 

35  45 

985 

FSE 

16  51 

t 

090 

t3B 

RCLC 

36  14 

986 

RCL3 

36  93 

082 

831 

1 

01 

987 

b 

06 

832 

— 

-45 

986 

£ 

00 

(TO  GET  N.MI.) 

833 

SWD 

35  14 

989 

8 

08 

834 

lt-43 

o«o 

999 

0 

00 

B35 

CToe 

33  00 

991 

X 

-34 

836 

ISZI 

It  3t  40 

992 

PR7X 

-14 

083  X 

837 

tSZl 

ib  £b  4b 

993 

RCL7 

36  07 

089 

JS 

CTOC 

£3  13 

994 

PRTX 

-14 

__y 

839 

tLBLe 

31  00 

(SEE  FROM  20) 

995 

6SB9 

23  0? 

B«0 

849 

ecu 

36  15 

896 

STOD 

35  14 

V 

941 

STDv 

35  14 

997 

PRTX 

-14 

842 

ftae 

36  00 

998 

F2S 

io'Si 

843 

SUl 

35-55  01 

999 

RCL5 

36  05 

V 

944 

tSB2 

23  &2 

100 

199 

prs 

lt-51 

945 

3 

03 

191 

RCLD 

36  14 

846 

STOI 

35  46 

192 

•f 

-34 

947 

croi 

33  01 

193 

CHS 

■-22 

948 

3i  01 

(SEE  FRGM  2^ 

1 

104 

SIH-’ 

16  41 

e 

949 

CSBi 

33  45 

105 

FRTX 

-14 

060 

850 

ISZI 

16  26  4t 

lOt 

6706 

22  11 

851 

BCLi 

36  45 

103 

•L6L3 

31  03 

852 

♦ 

-55 

198 

P2S 

lt-51 

(SEC) 

853 

RCLO 

36  90 

199 

RCLI 

3t  01 

854 

X 

-35 

110 

lie 

P.^S 

16-51 

(PRO 

855 

Z 

03 

111 

R7N 

34 

u 

956 

X 

-35 

ML. 

•I.BL7 

31  or 

REQISTtRS 


0 

) 

? 

3 

4 

s 

8 

9 

so 

Si 

SI 

S3 

S4 

ss 

S6 

S7 

S8 

S9 

A 

B 

C 

0 

E 

t 

4fi‘  'PI'.  Pf 


i ...  IIJIIV'l'. 


’TzrrrzT 
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\ 

1-. L-_  _ KKFERKNCES 

a.  R.  H.  Frick,  W.  1.  Rumor,  and  K.  H.  Sharkey,  ii>ui 

Orbit  Plottrr  I urtructicn  Mumnil,  The  Rand  Corporation, 
R-A18-1’R,  October  1963. 

h.  Opai\:  Ouidc,  USAF  Air  Force  Systems  Command,  1 July 

1965  (For  Official  Use  Only). 


The  orb i t a 1 el ement s are 


1 . 

Injection  altitude 

( n m i ) 

»1 

2 . 

Injection  velocity 

(ft /sec) 

"l 

3. 

Injection  flight  path  angle 

(deg) 

^1 

4. 

Per lod 

(min) 

T 

5. 

Kccent  r ic i ty 

(dimensionless)  C 

6. 

Semi-major  axis 

(n  mi) 

A 

7. 

Perigee  altitude 

(n  mi) 

11 

p 

8. 

Apogee  altitude 

(n  mi) 

11 

a 

9. 

True  anomaly 

(deg) 

^’l 

(at  injection  measured 
from  .apogee) 

The 

program  below  solves  two  of 

m.'iny  possihlt' 

problems 

A.  Civen  1,  2,  3 

f ind 

tlie 

rema ining 

element 

B.  Civen  1,  7,  8 

f ind 

tlio 

remaining 

element 

Reference  b uses  a sequence  of  nomograms  to  solve  these  problems. 
Some  of  these  nomograms  are  difficult  to  read  with  any  precision  be- 
cause of  close  Interval  spacing  for  relatively  large  increments. 


5.3._  KqUAXlONS 

In  Fig.  5.1  distances  arc  measured  in  units  of  earth's  radius 

(3437.9  n mi)  and  velocity  in  units  of  /lAg,-,  (25,943  ft/sec),  the 

L U 
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surface  orbltax  velocity,  and  the  surface  gravitational  accelera- 

tion (32.174  ft/sec^). 


r = R/Rj,  . 


C =*  RjV'j  cos  Yj  * KV  cos  Y . 


V/r^ 

c » c/v/r^ 


K .x:  "o5 

^'0  ’ 2 R 


E = E^/R^,«q 


2 

P = 


The  orbital  period  is 


T„  = 2tt  l-2Er^/" 


0 


Other  pertinent  relations  are 


- = ^ - p cos(0  - V ) 
r 2 I 

c 


tan 


Y = -r  p sin(0  - v^) 


JL  = i _ p 

r 2 
a c 


1 1 


p c 


+ P 


where  r^,  r are  the  apogee  and  perigee  distances. 


r - r 


e = --——-1?  (eccentr  ic  it v) 
r + r 
a p 


(1) 

(2) 

(3) 

(4) 


(5) 


(b) 


17) 


(«a,8b) 


(0) 


r + r 

a ” — ^ (semi -major  axis) 


(10) 
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5,4.  PROGRAM  NOTES 

Flag  F3  (set  by  digit  entry,  cleared  by  test)  is  used  to  direct 
the  program  to  Problems  A or  B.  The  program  flow  for  Problem  A (equa- 
tion numbers  in  parentheses)  is: 


1. 

LBL  A,  GTO  LBL  0 (1),  r^ 

2. 

LBL  B,  GTO  LBL  1 (1),  v^. 

3. 

LBL  C,  GTO  LBL  2, 

4. 

LBL 

e (2),  1/C^,  (3),  |e|  (total  energy/non-dim.) 

5. 

LBL 

D (8b),  Hp 

6. 

LBL 

E (8a),  H 

7. 

LBL 

a (5)  Tq 

8. 

LBL 

b (9),  e 

9. 

LBL 

c (10),  A 

10. 

LBL 

d (7),  Vj 

This  program  can  be  appreciably  shortened  by  a judicious  use  of 
subroutines.  A suggested  exercise  is  to  rewrite  it  to  see  liow  many 
problems  other  tlian  A and  B can  be  packed  on  one  card. 
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A:  H, 
B:  Hj 


5.5  USER  INSfRUCTlONS 


SATELLITE  ORBITAL  ELEMENTS 


STEP!  PROBLEM  A INSTRUCTIONS 


KEY  Hi(n.  mi.),  PRESS  A,  SEE  r. 


KEY  V,  (f.p.s.),  PRESS  B,  SEE  v 


KEYy,  (degs),  PRESS  C,  SEE  y 


PRESS  f e 


PRESS  D,  SEE  Hp  (n.mi.) 


JW^S_ E,  SEE  Ha  (n.mi.) 

PRESS  f o,  SEE  Tq  (mins) 

PRESS  f b,  SEE  e ^ 

PRESS  f c,  SEE  A (n.mi.) 

P^SS  f d,~lE7  *^^gs) 


INPUT 

OATA/UNITS 


OUTPUT 

OATA/UNITS 


rejss  C , SEE  y, 
PRESS  B,~  SEE  Vi 
PRESS  R/S 
PRESS  f g,  SEE  Tg 
PRESS  f b,~  SEE  g 
PRESS  f c,  SEE  A~ 
PRESS  > d“  SEE  V. 


DATA  CARD  o.coodeseee 

— e.  cmeeeee 

?.  emeeeee 

C. 090986066 
0.009060000 

0.000000000 

(?.  000000080 

3437.900000 

3Z933. 32088 

3=. S4260000 

0.000000000 

0.000000000 

-0.000000000 

0,000000000 

0. 000000000 

0. 000000000 

I 


rz  j zn 

ZZILZZ 
[ZZ][ZZ 
zz  zz 
r*i  rzzi 


46.534 


405.442 


92.866 


1.175 

2.273 

26047 

.454 

97.62 

.066 

3787.9 

39.203 


PRICXUNO  PiOS  auunc 
i 
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6.  SATKLLITK  TRACKING 


b . I . RKFKRKNCK 

a.  R.  Honsou,  "Comput  lt  izetl  Sati'Ilito  Trai'king, " 7?  rkiijaiv:  nt’ , 
I'obriiary  1^)77. 

■k 

Amatoiir  radio  operators  make  extensive  use  i>t'  OSCARS  (Orbiting 
Satellites  Carrying  Amateur  Radio)  for  long-range  eommun ieat ions . 
i'he  OSCARS  are  in  near-circular,  sun-synchronous  (lolar  orbits  at  an 
altitude  of  about  1500  km.  OSthXR  comes  within  range  of  any  given  spot 
on  earth  twice  each  day,  local  morning  heading  south  and  local  eve- 
ning heading  north,  and  may  be  within  range  for  as  mucli  as  25  min. 

An  ephomeris  provides  time  and  longitude  (if  orbital  equatorial 
crossings  (KQX).  The  problem  is:  t’.iven  the  latitude  and  longitude 
of  a ground  station,  determine  whether  or  not  a given  orbit  can  be 
viewed  and,  if  so,  find  the  time  the  bird  rises  over  the  horizon  and 
from  then  on  determine  its  range,  bearing,  and  elevation  at  desired 
times  until  it  disappears  below  the  horizon.  Figures  b.l  and  b.2, 
prepared  by  W.  B.  Graham  (KbQB)  of  The  Rand  Corporation,  show  an  ele- 
gant solution  for  his  station  in  Pacific  Palisades,  California.  (See 
examp  I es  be  low. ) 

The  problem  is  clearly  of  militarv  ;(S  well  as  amateur  interest. 

equations 

Notat ion : 


A = station  latitude  (north  only) 

B = station  longitude  (-  if  west) 

S = orbital  inclination  measured  c('>unterc lockwise  from 
the  equation  (•-90°  posigrade,  '•90°  retrograde) 

A 

For  details  on  present  and  future  OSCARS,  see  a series  of  articles 
In  QST  magazine  (ARRL)  starting  with  the  Januarv  1977  issue  (Vol.  bl. 

No.  1). 


Elevation  (deg.  above  horizon) 


($t.  tni 


IIIKR  ■ 
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H = orbital  altitude  (n  mi) 

T = orbital  period  (min) 

a * latitude  of  subsatellltc  point  (SSP)  at  time  t 
Y = longitude  of  SSP  at  t 
t = time  from  EQX  (min) 

Yq  = longitude  of  EQX  northbound  (-  if  west,  t = 0 at 
this  node) 

R„  = earth's  radius  (3437.9  n ml) 

K 

Dp  = arc  from  station  to  ascending  node  (EQX) 

Dj^  = arc  to  point  of  tangency  of  suborbit  witli  latitude  B 
l)j^  = arc  to  SSP  of  closest  approach 
D = arc  to  SSP 

t^  = time  of  closest  approacl\ 
t^  + At^  = time  to  rise  above  (go  below)  the  horizon 
0 = bearing  from  station  north 
<{>  = elevation  above  horizon 
R = range  from  station  to  satellite  (n  ml). 

By  the  law  of  sines, 

oi(t)  = sin  j^sin  B * sin  (360  t/T)J  . (1) 

By  the  law  of  cosines, 

Y(t)  =*  cos  ' l^cos  (360  t/T)/cos  rt(t)j  “ Yp  * (-) 

where  t/4  is  the  correction  for  the  eartli's  rotation.  Also  by  tlie 
law  of  cosines, 

D(t)  - cos~*  [sin  A • sin  a + cos  A • cos  a • cos  (B  - Y)]  • ( 1) 


The  equations  for  0,  <)),  and  R are: 


0(t)  = cos 


(j)(t)  = tan 


-1 


In  g - sin  A • cos  D 
cos  A • sin  D J ‘ 


cos  D - 1/(1  + H/R^) 

hi 


sin  D 


(4) 

(5) 


R(t)  = (H  + Rg)  sin  D/cos  ())  . 


(6) 


The  above  equations  are  essentially  those  of  Ref.  a,  except  for  (6). 

Critical  times  are  dealt  with  by  the  following  equations  (which 
are  somewhat  in  error  because  their  derivation  assumes  the  suborbital 
trace  of  the  satellite  is  a great  circle,  which  is  not  true  on  a 
rotating  earth).  For  the  time  t^  of  closest  approach. 


T 

360 


tan 


-1 


(cos  D^/cos  Dq) 


(7) 


and  the  incremental  times  to  zero  elevation  are 


where  D is  the  arc  length  to  the  SSP  at  zero  elevation  and 


cos  D = 1/(1  -t-  H/Rg) 


Finally,  the  period  in  minutes  is 


(8) 


(9) 


(10) 


6.4.  PROGRAM  NOTES 

1.  045,  046:  Change  8 to  first  quadrant  to  get  latitude  of 

orbit  tangency. 
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2.  115,  116:  For  retrograde  orbit,  subtract  In  (2), 

'3.  117,  118,  119:  These  steps  produce  the  sgn  of  t (i  1),  since 

the  HP-67  does  not  iiave  this  signum  function. 

4.  120,  121:  For  t 0,  south  of  the  equator,  im)ve  time  back- 
ward so  that  posigrade  and  retrograde  orbits  reverse  in  getting  SSP 
longitude  by  (2). 

5.  LBL  corrects  longitude  if  outside  (-180,  +180),  y > 180 

becomes  -360  + Y -180  becomes  360  + y* 

6.  l.BL  1 gets  bearing  clockwise  from  station  nortli. 

Kxamp 1 e 1.  In  Fig.  6.1,  the  ground  station  is  at  34°03'N  and 
118“33'W.  Tlie  satellite  is  retrograde  at  an  altitude  of  790  n mi 
and  has  an  inclination  of  102°.  Track  the  satellite  when  in  view  if 
the  EQX  is  llO^W  at  time  0. 

5^  1 uT  ion . 

34.03  f H (34.05)  STO  A;  118.33  f H (118.55)  CHS 

STO  B;  790  STO  C;  102  STO  D;  100  CHS  STO  K; 

-1  h STl. 

A:  T = 115.11,  t^  = 11.68,  -29.74  (will  come  in  view). 

R/S:  RCL  5,  Y = -100.39;  RCL  4,  a = 1.31  (SSP) 

R/S:  RCl.  3,  t = 0.043;  RCI.  7,  0 = 148.72,  RCL  8,  4)  = -1.27 

Time  t for  first  appearance  is  slightly  small.  Try  t = 1 
STO  3,  PRESS  B,  RCL  8 and  get  (^  = 0.51.  Try  t = 0.8  and  get 
-0.12,  which  is  good  enough.  The  following  table  is  prepared 
by  successively  storing  t in  3 and  recalling  7,  8,  9: 


iil). 

II 

B 

1.1515  stat  mi) 

0.8 

148.40 

-0.12 

2842 

4 

144.56 

11.47 

2154 

8 

132.20 

34.11 

1380 

12 

62.63 

60.68 

1014 

16 

1.62 

31.70 

1435 

22 

348.65 

2.69 

2654. 

‘ values 

are  plott 

ed  as  solid  bullets  in 

Fig.  6.1. 
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Example  2.  All  values  are  those  of  Example  1,  except  the  EQX  is  70"E. 
Hence  the  satellite  will  approach  Pacific  Palisades  from  the  north. 

^JLut^qn.  Proceed  as  before.  Note  T/2  = 57.5  min  and  STO  3,  Press  B, 
RCL  7,  0 = -124.39  is  the  descending  node  (KQ.X,  moving  south).  Time 
will  now  be  measured  positive  to  the  south  of  tlie  equator.  STO 
-124.39  in  register  E.  Change  the  sign  of  ttie  inclination  3 to  minus 
and  STO  D.  Now  go  back  to  square  one  and  redo  the  problem  from  LBL  A 
on.  The  results  iire  tabulated  below  and  again  plotted  as  solid  bullets 
in  Fig.  6.2  for  comparison. 


t (corr . ) 

0 

R 

-22 

35.5 

19.14 

-.19 

2847 

-18 

39.5 

23.88 

14.95 

1991 

-14 

43.5 

36.65 

41.33 

1241 

-12 

45.5 

61.85 

61.46 

1008 

-10 

47.5 

136.24 

65.43 

981 

-6 

51.5 

181.04 

28.90 

1506 

_2 

55.5 

188.51 

8.23 

2323 

+ .5 

58.0 

190.72 

-.47 

2866 

(stat  mi) 


Ephemeris  EQXs  are  given  in  ZULU  time  (CMT  or  UTC  - Universal 
Coordinated  Time  for  the  radio  amateur).  Correct  for  local  station 
t ime . 


6.5  USER  INSTRUCTIONS 


INSTRUCTIONS 


STO  DATA  IN  Rq,  Ri,  R^ 

STATION  LATITUDE,  STO  A (CONVERT  TO 
DEC  PEGS  BY  fH  IF  NECESSARY) 

STXtION  longitude  ^O  B 1 - FOR  ^ LONG,' 
CONVERT  TO  DEC  DECS  BY  m ^ ' 

ALTITUDE  (iu  mi7j,  STO  C ^ ' 

ORBJ^TAL  INCLfr^TIONTlT'O  D (CC  FROM  ' 
EQUATOR)  ’ y ^ 

EQX  LONGITUCrE,  STO  E ' 


FOR  POSIGRADE^I  ST  ^ 
FOR  RETROGRADE,  -1  h ST  1 


INPUT 

OATA/UMTS 


OUTPUT 

DATA/UMTS 


PRESS  A,  FIRST  f-x-SHOWS  T 

. SECOND  f-x-SHOWS  TIME  OF  CLOSEST 
APPROACH 

10  ON  R/S,  IF  NUMBER  IS  POSITIVE'^TELLITE 
WILL  NOT  COME  IN  VIEW. 


11  PRES^R/S.  ON  HALT^RCl^STO  G^ 
APPROX  LONG  OF  SSP^N 'first  _ 
APPEARANCE.  RECc5rD^  RCL  4 TO  GET 
LATITUDE  OF  SSP. 


12  PRESS  R/S.  ON  HALT,  RANGE  IS 

DjSPLAYED".  RCL  3 TO  GET  TIME , RCL  T_ 

TO  GET  BEARING,  RCL  8 TO  GET  ELEVATION 


■ n: 

r 

-24 

■ 143 

KCL2 

36  92 

' 143 

X 

-35 

' 144 

CHS 

-22 

!4r 

RCL5 

36  95 

■ Nc 

* 

-55 

■ ir 

ST05 

35  95 

' NS 

RTN 

24 

■RT" 

irsr 

' 15S 

RCLB 

36  12 

151 

RCL5 

36  95 

15' 

- 

-45 

15  j 

COS 

42 

15- 

RCL4 

36  94 

155 

COS 

42 

I5t: 

X 

-35 

* 

FCLH 

36  11 

15J 

COS 

42 

15" 

X 

-35 

ise 

RCL4 

36  94 

Icl 

SIN 

41 

If? 

RCLA 

36  11 

15' 

SIN 

41 

15-' 

X 

-35 

155 

♦ 

-55 

leo 

CCS-' 

16  42 

157 

ST05 

35  96 

155 

RTN 

24 

FLAGS 


ON  OFF 

0 □ o 

1 □ □ 
2 D n 
son 


SET  STATUS 


TRIG 


DEG  □ 
GRAD  D 
RAD  D 


FIX  □ 
SCI  □ 
ENG  rj 
n 
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7.  THE  DEER  HUNT  (DEFENSELESS  BOMBERS) 


7.1.  REFERENCES 

a.  C.  H.  Builder,  The  Penetration  Integral  and  Tables,  The  Rand 
Corporation,  R-1257-PR,  June  1973. 

b.  N.J.J.  Bailey,  The  Elements  of  Stoahastia  Processes,  John 
Wiley  and  Sons,  New  York,  1964. 

c.  A.  T.  Bharucha-Reid , Elements  of  the  Theory  of  Markov  Pro- 
cesses and  Their  Applications,  McGraw-Hill,  New  York,  1960. 


7.2.  DISCUSSION 

Reference  a deduces  an  integral  and  constructs  tables  to  assess 
the  expected  outcome  of  a one-sided,  time-limited  battle  in  which  a 
set  of  armament-limited  interceptors  engages  a set  of  defenseless 
penetrating  bombers.  The  formulation  also  applies  to  a set  of  vessels 
transiting  a minefield  where  each  mine  has  a fixed  number  of  warheads 
of  some  description.  There  are  undoubtedly  other  military  applications. 

Builder's  original  paradigm  is  preferred.  At  the  beginning  of  a 
hunt,  there  are  A deer  and  B hunters,  each  of  the  hunters  armed  with 
m rounds  of  ammunition.  Encounters  are  at  random  with  parameter  X. 

A hunter  can  expend  only  one  round  on  each  engagement,  with  kill  prob- 
ability p.  The  hunt  lasts  T units  of  time. 

In  Ref.  a,  the  parameters  adopted  are 

j = XpBT  k = pmB/A  , (1) 

which  provide  a bridge  to  this  discussion.  The  parameter  j is  the 
potential  number  of  lethal  encounters  per  deer  during  the  hunt,  based 
upon  the  expected  encounter  rate,  while  k is  the  potential  number  Of 
lethal  encounters  per  deer  based  upon  the  total  hunter  armament . The 
parameters  X and  p are  perhaps  the  more  natural  ones  to  employ. 

A much  simpler  approach  to  finding  the  expected  values  than 
that  of  Ref.  a is  adopted  here.  This  approach  is  adequately  illus- 
trated by  the  cases  m = 1,  m = 2. 
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For  m = 1, 


da 

dt 


-Apab 


dt 


= -Xab  , 


(2) 


where  a and  b are,  respectively,  the  number  of  deer  and  the  number 
of  armed  hunters  remaining  at  time  t.  Division  gives  a first  inte- 
gral a = pb  + a,  where  a = A - pB.  Whence 


A(A  - pB) 

A - pB  exp  {-(A  - pB)  At}  ’ 


(3) 


and  a(t)/A  is  precisely  the  of  Ref.  a if  (1)  is  used. 

The  a and  b found  by  (2)  were  called  the  expected  values.  Ac- 
tually they  are  the  deterministic  values.  Using  the  methods  of  Ref. 
b (p.  118),  which  set  up  a partial  differential  equation  for  the 
moment-generating  function. 


du 


10 


dt 


-Apu 


11 


01 


dt 


= -Ay 


11 


(4) 


where  Pji^qj  are,  respectively,  the  means  or  expected  values  of  a 

and  b,  and  y^^  is  the  correlation  between  a and  b.  Hence  for  (4)  to 
agree  with  (2)  we  must  have 


'11 


10 


01 


or  Exp  (ab)  = Exp  (a)  • Exp  (b) . This  is  not  true  because  the  popu- 
lation sizes  are  mutually  dependent.  (See  Ref.  c,  p.  184.)  Conse- 
quently, we  proceed  with  the  understanding  that  we  are  dealing  with 
deterministic  values  rather  than  expected  values. 

Turn  to  the  case  m = 2.  At  time  t the  values  are  a,  b^^,  b2, 
where  bj^  hunters  have  one  round  left  and  b2  still  have  one  round 
pouched.  The  deterministic  equations  are: 


2 • 


= -Apa(bj^  + b2)  , 


= -Aabj^  + Aab2  , 


= -Aab 


(5) 
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Two  integrals  are  found  immediately: 


a - pb^  - 2pb^  = A - 2pb 


>1  = -N  ^-N/b 


Whence, 


dT  " -Xb,(A  - 2pB)  . (7) 

For  particular  values  of  the  parameters,  (7)  is  readily  integrated  by 
Program  11  of  the  HP-67  Math  Pac  1.  The  results  agree  exactly  with 
those  of  the  table  for  m = 2 in  Ref.  a. 

In  principle,  for  m > 2 numerical  integration  is  possible.  In 
practice  this  would  be  very  time-consuming,  if  not  infeasible,  for 

A 

the  HP-b7.  The  next  section  provides  a completely  different,  al- 
though hcui'it^tic; , approach. 

7.3.  EQUATIONS 

Again  take  m = 2.  (The  symmetric  equations  to  be  derived  are 
readily  extended  to  a general  m.)  On  each  encounter,  A decreases  on 
the  average  by  p.  Hence  a(n)  = A - np  is  the  expected  number  of  deer 
just  after  the  nth  encounter.  The  average  time  to  the  next  encounter 
is  At  = l/(Xa(n)  • b(n)),  where  b(n)  = b^(n)  + b2(n)  is  tlie  remaining 
number  of  armed  liunters  after  the  nth  encounter.  Then 

b,(n  + 1)  = b (n)  - 1 , b,,  (n  + 1)  = b 

11  2 n 


with  probability  b^(n)/b(n),  and 


bj^(n  + 1)  = b^(n)  + 1 , b^(n)  = b^(n) 


with  probability  b2(n)/b(n), 


But  see  Sec.  20  for  m = 3,  4 application. 





Hi'ni'o  t lio  r.r/’i'.’f  i’.;  >\t'xt  ts 

li|(u  + 1)  - |h|(n)  - l|  • t'j  {i»)/l'(ii) 

♦ |l>|  (n)  ♦ I I • I’  ,(it>/t'(ci)  . 
h , t n * 1 > ' l>  , U» ' t’  j 0> ' I'  in  ' 


1 
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and 

bj  + 2b^  + 3b^  + 4b^  = 40.00  , 
which  is  the  expected  number  of  remaining  rounds. 
hA-.  PROCRAM  NOTES 

1.  The  program  is  designed  to  show  tlie  running  liistory  of  the 
hunt  by  successive  engagements.  A PAUSE  of  I sec  displays  the  engage- 
ment time  and  tlie  remaining  expected  number  of  targets.  If  more  time 
is  needed  to  record,  key  h SE  I,  which  will  give  a 5-sec  f -x-  fiasli- 
ing  toll()wed  by  a 1-sei'  I’AUSK.  At  any  time,  stopping  tlie  program  by 
K/S  pi-rmits  an  examination  ol  the  current  status  of  the  hunt  bv 

. UKOI  A)  . htKtl  HI  . b tKfl  8)  , 

4 


t Kl 


. » , 


h I Ml  I M 


7.5  USER  INSTRUCTIONS 


7.  THE  DEER  HUNT  (m  < 4) 

A (STO  A)  B.m  (STO  B)  X STO  C p STO  D 


«p- 


wymy!;? 


II.  hi  I 
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8.  A BOMBER  PENETRATION  MODEL  (DEFENDED  BOMBERS) 


8.1.  REFERENCES 

None . 

8.2.  DISCUSSION 

This  penetration  model  is  offered  solely  as  an  example  of  how 
the  HP-67  can  assist  the  analyst  in  his  preliminary  study  of  the  fac- 
tors bearing  on  a problem  and  how  these  factors  interact.  Model- 
making is,  or  rather  should  be,  an  art  form  drawing  the  essential 
elements  from  reality  and  illuminating  their  articulation.  The  author 
holds  to  the  view  that  an  initial  model  should  be  economical  and 
transparent.  The  ornaments  come  later. 

A group  of  bombers  with  integrated  fire  control  for  mutually 
supporting  self-defense  against  interceptors  makes  a corridor  pene- 
tration to  a set  of  targets.  The  payload  of  a bomber  can  be  divided 
at  pleasure  between  defense  missiles  (AAM)  and  ground  attack  munitions 
(ASM).  This  loading  is  decided  prior  to  the  mission  by  choosing  the 
number  of  AAM  to  be  fired  against  each  interceptor  based  on  the  ex- 
pected number  to  be  encountered.  At  equal  intervals  of  time  during 
the  penetration,  a clump  of  Interceptors  comes  within  range  of  the 
bt»nib»TH'  /VAMs.  'Iw  ix'ni>  ri>  /h'lr*.  Kat'h  of  t Ik-  surviving  inti-r- 

'cplors  ol  the  I liiiap  I*. ikes  a single  pass,  .illotating  fire  iinilormlv 
.» r t-i  iimI  t lien  wllltdriws  t !»*•  battle.  As  the  battle 
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K = A + r • S,  the  payload  constant 
r = ratio  of  ASM  to  AAM  weights 
k = AAMs  expended  per  interceptor  encountered 
Q = probability  an  AAM  will  miss  an  interceptor 
1 = interceptors  per  clump 

q = probability  an  interceptor  will  miss  its  target 

N = duration  of  the  penetration  with  unit  of  time 
the  interval  between  interceptor  mass  attacks 

T = ASMS  on  target 

y**x  = y^. 

Then 

. ' B . = B • q**(I  • Q**k/B  ) (1) 

n+1  n ^ n 

because  I • Q**k  interceptors  survive  and  this  divided  by  B^  is  the 

expected  number  of  passes  per  bomber. 

The  number  of  missiles  fired  per  bomber  on  the  nth  engagement  equals 

kl/B  . Hence  tlie  number  of  AAMs  per  surviving  bomber  to  target  is 
n 

A » kl(l/B|  + 1/B,  + ...  + l/!lj^)  , li) 

' t u I liniiink:  .itmimpt  1 iui  . w«  rc  rc.i  I I /I’d  In  comb. it  . liie  luim'cr  ot 

. • 1 l I III  it  < illl  I I , ( II.  I I It  till  lit  I'llC  i . I 

s» 

» « t 


I .iij  V ‘."vt 


"^SB" 
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quickly,  and  (2)  this  exploration  is  informative  since  it  shows  the 
sensitivity  of  the  outcomes  to  k (as  well  as  to  the  other  parameters). 


Example. 


B = 10,  K = 30,  r = 3,  Q 


k 

2 

1 

3 

4 

3.5 


1/2,  I = 10,  q = 1/2,  N = 3 


A ®4  T 


7.30 

5.31 

40.19 

4.72 

1.95 

16.40 

9.87 

7.52 

50.47 

12.55 

8.73 

50.78 

11.19 

8.22 

51.54 

Note  how  flat  the  bombs-on-target  curve  is  for  k = 3 to  4.  Going  from 
k = 3.5  to  4 decreases  bombs  on  target  by  1.5  percent,  but  inaveascs 
bombers  saved  by  5 percent  of  the  original  force.  This  sr  ing  is  not 
trivial  if  bombers  are  to  be  recycled  for  follow-on  attacks.  It 
illustrates  the  insight  that  quickly  prepared  "toy"  models  can  provide. 


8.4.  PROGRAM  NOTES 


None . 


8.  A BOMBER  PENETRATION  MODEL 


OUTPUT 

OATA/UMTt 


-88- 


8.6  BOMBER  PENETRATION 


» KEY  ENTRY 


.801 

.082 

.083 

.00-f 

IB 

.805 

a 

.088 

007 

r 

0OS 

.009 

- 

.318 

011 

012 

JUi 

Oil 

015 

01o 

017 

013 

01.3 

ga 

- 

028 

"321 

822 

023 

021 

025 

02o 

027 

I 023 

‘82.- 

- 

038 

'031 

'032 

■033 

Oo-f 

■035 

■035 

"037 

■033 

!!9 



■039 

KEY  CODE 

li 
J3  il 
31 
33  K 
31 
33  «3 
31 
33  84 
31 
33  83 

31 
33  88 

i±_ 

21  U 
33  89 
38  12 
33  88 

32 
33  97 

81 
33  48 

22  13 

irrr 

18  28  48 
38  83 
38  89 
31 
38  84 
-35 
38  88 
-24 
38  85 
-41 
31 
38  88 
-35 
35  88 
52 

33-55  87 
38  88 
81  ■ 


COMMENTS 

INPUTS 

K 


STEP 

KEY  ENTRY 

KEY  CODE 

057 

CHS 

-22 

-053 

RCLl 

38  81 

059 

+ 

-55 

060 

'858 

RCL2 

38  82 

051 

-24 

052 

RCL0 

38  88 

'853 

X 

-35 

‘051 

STOC 

35  13 

■055 

RTN 

24 

1/B, 

1 hsfl 


TO  GET  yX 


■012 

♦ 

-55' 

TO  GET  N -H 

013 

RCLI 

38  48 

PEN.  COMPLETE  ? 

■011 

18-33  ■ 

■015 

CTOD 

22  14" 

LOOP 

■015 

eroc 

22  13 

• 2810 

m? 

01E 

RCLO 

38  88~ 

6 u . 

■01? 

1 X 

52' 

■050 

ST-7 

35-43  87  ■ 

031 

»Ct7 

J8  87 

852 

RCll 

JB  84" 

-33  ■ 

»«i 

M 18 

>5* 

•» 



^8N 

rii 

» II 



9.  DAMAGE  PROBABILITIES.  PVN  AND  QVN  TARGETS 


9.1,  REFERENCES 


D.  C.  Kephart,  Some  Aids  for  Estimating  Damage  Probabilities 
in  Attacks  Against  Targets  with  P and  Q Vulnerability  Numbers 
The  Rand  Corporation,  R-1168-PR,  March  1973  (For  Official 
Use  Only) . 

Physical  Vulnerability  Handbook- -Nuclear  Weapons,  Defense 
Intelligence  Agency,  1975  update. 


9.2.  DISCUSSION 


The  program  in  this  section  was  prepared  by  D.  C.  Kephart  of  The 
Rand  Corporation.  He  has  alsu  written  the  program  for  Texas  Instru- 
ments' SR-52  hand  calculator. 

The  program  gives  daimige  probabilities  for  nuclear  weapons 
applied  against  PVN  and  QVN  point  targets  at  the  optimal  airburst 
height.  For  these  two  classes  of  targets  'psi'  is  given  by: 


Dynamic  pressure,  psi  = 0.02893  x 1.44  (v  = adjusted  QVN) 


9.3.  EQUATIONS 


Notation 


VN.K  = Vulnerability  number 

V = Integer  part  of  VN.K 

K = Fractional  part  of  VN.K  = [K-factor ] /lO 
w = Warhead  yield  in  kilotons 
C = Weapon  CEP  in  feet 
A = VN  adjustment 

V = V + A adjusted  vulnerability  number 
R ■ Weapon  radius  in  feet 

P ■ SlngU'-shot  prob.'ibl  1 1 1 V of  damage  (SSPO) 
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Formulas  for  PVN 


S-f  j[f 


A = ln(S  )/ln(1.2) 


V = V + A 


R = w '-^[6383.35  x 0.8836'']  if  v ^ 20.5 


R = w '-^[1900.05  X 0.9368'']  if  v > 20.5 


P = 1 - exp  I -R‘'/2/[C^/ln(4)  + 0.04R^]| 


Formulas  for  QVN 


S satisfies  the  cubic  equation 


S = [SK  + 1 _ K]^''^ 

w 


A = ln(S^)ln(1.2^) 


V = V + A 


R = w^^^[6561  X 0.87918'']  if  v S 15.4 


R = w ' [23.42  + 2736.9  x 0.92883'']  if  v > 15.4 


t = 1 - exp  I -R^/2/[C“/ln(4)  + 0.09R^]( 


P = t if  t ^ 0.82 


Q - 2.826t  - n.94t  - 0.866  if  t ^ 0.82 


r • 0 It  0 < 1 


VN.K 

Sto 

A 

(K  = K factor,  21Q7  = 21.7,  42P6  = 42.6) 

Yield,  KT 

Sto 

B 

CEP,  ft 

Sto 

C 

PVN 

Key 

A 

^Adjusted  Vulnerability  Number  -pause- 

Calculate  < Weapon  radius  -pause- 

gvN 

Key 

E 

( Single-shot  probability  of  damage  (SSPD) 

Compute  time  10-12  sec  for  PVN,  12-22  sec  for  QVN 

Display  SSPD 

Key 

B 

Calculate  PD  for  one  more  shot 

Key  CEP,  ft 

Key 

C 

Calculate  SSPD  (After  Key  A or  Key  E;  quick  SSPDs 

using  new  CEP  values.  Cominite  time  3 sec.) 


EXAMPLE 


21.9 

Sto 

A 

1000  KT 

Sto 

B 

'000  ft 

Sto 

C 

Ad  i VN  Weapon  Radius 

SSPD 

VN  = 

21P9 

Key 

A 

12.509  13575.496 

0.97  3 

VN  = 

2iq9 

Key 

E 

17.253  7891.371 

0.732 

Disp lay 

SSPD 

Key 

B 

PD  = 0.928  (21Q9,  2 shots 

0 

Key 

6000 

Key 

C 

SSPD  = 0.627  (21Q9,  1000 

KT,  6000 

9,6  PKOGRAM  USTING 


STEP  Kcr  e/VTPV  KEV  CODE 


STEP  KEV  ENTRV  KEY  CODE 


i]J  cl  it  12  DSP  P (QVN)  (11) 

4 4 M n<M  n.-i  > \ / 


114  KCLS 

115  RTH 


-55  V (7) 
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WBl 


FOUR  DEUCKS  (PRKClSLON  A.2-LNCH  MORTAR  jiLFF) 


IQ RKFKRF.NCKS 

a.  FT  4.2-F-l,  Finug  TahLco,  Movtav,  •} .'2- >icih,  M.U),  Uepartment 
of  the  Army,  December  1954. 

b.  FM  23-92,  ‘1.2-Tnij'n  Movtav,  M60,  Department  of  tlie  Army, 
February  1961. 


10.2.  D1S(:U.SSL0N 

The  4.2-iiicli  mortar,  M30,  i.s  a rifled,  muzz  le- loaded  weapon,  known 
affectionately  to  the  Army  as  the  "four  deuce."  The  tube,  elevated  at 
anjtles  of  45“  to  60“  (800  mils  to  1065  mils)  delivers  indirect  fire  to 
almost  5.5  kilometers,  depending;  on  tlie  propellant  cliarge,  elevation, 
and  round  selected. 

Indirect  fire  units  use  a met eoro loj; i ca 1 message  (a  coded  weather 
report)  from  division  artillery  in  conjunction  with  unabridged  firing 
tables  to  prepare  firing  data.  As  of  1961  (Ref.  b)  a rather  time- 
consuming  and  largely  manual  procedure  was  used,  whicli  even  sacrificed 
almost  all  tabular  interpolation  to  save  time  and  avoid  errors.  Tliis 
section  sltows  how  tlie  procedure  of  tliat  era  would  liave  been  simplified 
liad  a programmable  hand-calculator  been  available  tlien. 

The  primary  task  is  to  reduce  the  tables,  whicli  were  based  on 
range  firings  conducted  at  the  Aberdeen  Proving  Crounds,  to  a set  of 
formulas.  This  is  accomplished  by  data  fitting,  a task  for  which  the 
HP-67  is  admirably  suited  if  one  has  at  hand  Program  3 of  the  Standard 
Pac  and  Progr.am  14  (Polynomial  Approximation)  of  Stat  Pac  1.  Mark 
that  this  data  fitting  is  purely  empirical.  ft  is  based  in  no  way  on 
the  physics  and  mathematics  of  exterior  ballistics. 

The  only  problem  posed  by  the  tailoring  of  formulas  to  number 
streams  is  the  choice  of  the  formula  type  to  be  used.  Initial  guidance 
is  provided  by  plotting  families  of  curves  and  staring  at  them.  (See 
Sec.  21.) 

The  resulting  formulas  are  simple  and  so  is  the  required  program- 
ming. You  pay  for  this  double  simplicity  but  can  exploit  it.  The 
large  number  of  constants  generated  by  fitting  exceeds  available 
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stora^o  space.  But  because  the  program  Is  short,  there  Is  space  to 
put  constants,  sometimes  rounded,  in  the  pronram  Itself.  This  is 
program/storage  tradeof f . 

Tlie  next  subsection  gives  formulas  that  correct  the  fire  for 
nonstand.ird  conditions,  in  the  order  in  which  they  will  be  programmed. 
This  Is  also  the  approximate  order  of  the  manual  calculations  In  tlie 
examples  of  Ref.  a.  The  program  significantly  modifies  the  methods 
of  that  reference  In  respect  to  automatic  interpolations,  .illowance 
of  difference  In  altitude  of  mortar  and  target,  l>alllstic  winds,  .ind 
elevat ion  correct  ions. 


Cenej^n 

(1)  This  program  is  restricted  to  the  MiO  firing  the  UK  shell 
MJ29  wltli  extension  (long-range  fire). 

(2)  Tl>e  r.inge  of  charges  is  2‘>.5  to  41  and  increments  of  l/R 
are  permitted. 

(J)  Klevat  ions  are  restricted  to  800  to  dOO  mils  (\A) , since 
the  tables  are  so  restricted  for  charges  above  32. 

(4)  Hence  charge  is  selected  so  tliat  850  ^ (47.81°)  elevation 

giv»>s  the  desired  range  ( 1 fli  = 3b0/b400°  = 

0.05b25  = 1/17.778). 

(5)  Meters  r.ulier  than  yards  are  used  (1  yd  = 0.'1144  m) . 


No tat  ion 


U^  Mtltude  of  mortar  position  in  meters 

R Chart  range  in  meters 
o 

Rj  Range  corrected  for  difference  in  altitude  of  mortar  and 
target 

R.,  R|  corrected  for  metro  and  ballistic  factors 

R(800)  Range  for  K * 800  i(i,  a function  i>f  charge  m 

A Azimuth  of  fire  (mortar  to  target)  in  mils,  t’.W  from  N 
o 


tho  initial  olovatlon 

K,,  ('orr»>ct»'d  olovatlon  (#1) 

C Tho  iitlt  lal  oharso  solootod 
o 

C.,  Corroctod  oharjte 

u)  Tho  anulo  of  fall  In  mils  (impact  an^lo) 

T Powdor  temporal uro  In  K” 

VH  Muzzle  velocity  error  in  ft/soc  fv>r  lot  used  ((1  it  not  known. 
This  can  ho  ohtalnod  onlv  t n'm  iri.il  llrin>;s.) 

r Ooviallon  in  sholl  woip.ht  t tv'm  [•[]  [•]  l-l,  0.  +11. 

(Weight  is  mo.jsurod  in  sipi.iii's  and  st  ampi'd  on  tho  sholl.) 

R|.  Altitudo  of  l argot  in  m«'t«MS 

lij^j  Altitudo  of  moteorological  d.ita  pl.ino  (NIX')  in  .''iT.' 

/Vy  Azinuith  of  ballistic  wind  in  ^ I'W  t rom  N.  (T'hls  is  the 
direction  which  the  wind  blows.  It  is  an  avor.igo  of 

winds  up  ti'  tho  maximum  t'rdlnato  of  tho  t ra  loctv'rv . ) 

Anglo  hetwoon  A^^  .»nd  .\y  - UOt)  in  dogia'os 

hensitv  ot  tho  .air  as  porcont  I'l  st.ind.ird  fi'r  tho  MOP  .lit  itiui 
<5^  llorrectod  densltv  fv'r  mortar  altitudo  rol.ilivo  t »'  MOP 
W Ballistic  wind  in 

d Hrlft  deflection  because  of  sholl  rot  it  ion  (.ilwavs  ti'  tho 
r I gilt ) 

I)  Deflection  correct  ii'ii,  fin.il  (#1) 

Pj  Correction  for  range  wind  (tall  or  head),  motors 

O,  Correction  for  round  weight,  motors 

Pj  Correction  for  powdor  tomperaturo,  motors 

P,  Ciirroction  for  actual  air  donsitv,  motors 

■4 

P^  Correction  for  VC.  meters 
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P " - (,Pj  + P»  + Pj  *-’5)  • 

P Clian>;o  la  muzzle  velocity  due  to  change  iu  pi'wdt'r  temper- 
ature (,tt/sec) 

Mey_o 

/V^,  W,  are  known  when  the  Metro  Messa>{e  is  "si'lved."  The 
message  also  >tives  air  temperature,  which  is  not  r«‘ levant  tor  mortar 
calculations.  The  message  has  12  lines  in  addition  to  the  lu-ading. 

The  initial  digit  is  the  standard  altitude  number.  The  correspondence 


Alt.  No. 

0 

1 

) 

i 

4 

S 

Height  (ft) 

0 

600 

1 soo 

loot) 

4 SOO 

6000 

Alt.  No. 

6 

7 

8 

9 

0 

1 

Height  (ft) 

9000 

12000 

1 SOOO 

18000 

24000 

lOOOO 

Tlu'  lines  of  the  message  give  .1^,  W,  5 f.ind  air  temperature) 
appropriate  for  the  maximum  ordinate  of  the  tralectory  resulting  fri>m 
any  particular  combiiiatlon  of  range,  charge,  and  eli'v.it  iiui. 

For  this  progr.tm  (800  ^ K £ 900  rfi,  R 2:  12S0  m)  , alw.tys  use  line 

4 

C - 2S.‘>  l.lne  J 

C » 41  l.iite  '> 

C S lb  and  K S 8‘>0  l.lne  “>  . 

Formulas 


F’  " 8. SO  lii 
o 


C - 1 1.9 12  exp  iO. 00021 1 R ) 
o ' o 


U)  - (867  + 4C  )(  160/6400) 

• o 


R , - R + (IL.  - H )/tan  u) 

1 o To 

(Rj  > R^^  if  target  is  above  mortar) 
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R_,  ■ Rj  + P 


This  Is  tlu'  corroitoii  final  ran>;o. 


11.4)2  I'xp  (0.0002)1  Rt> 

(Tliis  is  tiu'  cor roi' t Oil  iTi.ir^c  for  K » )!>0  rfi. 
It  is  r.'  hi'  to  the  no.ircst  1/2 

ch.ir^e.  ) 


K(800)  = 52.11  + 81.7)  0, 


+ 2.4(11'.-,  - tl.Oal  (■  , 


(Fitted  fri'm  bottom  line  J.it.t  of  t.ibles) 


K,  - 800  + 1(1. 41b  - 0.04)  ('.,1  1R(.800)  - R.,1 


This  is  the  fln.il  elevation  at  which  to  lay  the  tube.  Note  tlie 

eorrection  of  t'.  usinv  K.  = 850  tfi  .ind  tlte  tinal  R.,.  But  then  K = 850  i/i 
o ~ 

is  foiRotlen  and  the  procedure  is  to  ro  to  tormulas  to  Ret  K.^  wif^h 

respect  to  the  b.tseline  eli'vation  K = 800  rfi.  This  procedure  yields  a 
small  correction  to  the  initial  850  t(i. 

Note:  Dec.tuse  of  lack  of  proRram/storaRc  space,  formula  (11)  is 

rewritten  .is 


r = -(0.02  AT  + 0.005  Al'-)  AT  5 0 


0.05  AT  + 0.00)8  .AT' 


T ^ 0 . 


.ind  the  user  enters  AT  = T - 70°,  the  deviation  from  the  standard 
powder  temperature  with  the  proper  sIru. 


Checks 

Til  provide  priiRf.im  verification,  .i  st  ep-bv-step  e.x.imple  is  cal- 
cul.ited  m.inu.illv.  The.se  re.sults  are  then  compared  with  the  Army's 
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flold  proci'duro , which  docs  not  require  Interpolation  and  rounds  oil 
variiuis  values.  instead  ot  meters  are  employed  in  these  checks. 


r‘-  • SUO  (nearest  10  vd)  , H - 1505  it  . \L  - 1210  it  . 

q O r 

A -<i825  t . w>;t  - Q (r  - -1)  . T » 55”F  . 
u , 

\ 

VK  « - 12  it /'sec,. 


Solvin>;  the  Metro  Mess.i>;e  (heading  and  line  5), 

M 1 F 1 2 0 8 10  1 

- s 2 h 2 5 d 5 7 8 5 . 

we  >;et  : messa>;e  irom  Station  IF,  MOF  = 1200  it,  as  oi  0810  hr;  this 
is  ms>;  tvpe  1 (ior  mortars);  standard  altitude  ior  line  5 is  nOOO  it; 
ballistic  wind  blows  irom  2b00  rfi,  strenjtth  25  mplt;  air  densitv  is 
d5.7  percent,  .uui  air  temperature  is  SS'F. 


The  .»bove  conditions  were  taken  irom  Ref.  a.  Thev  could  be 
realised,  ior  example,  at  Hunter  l,i>;j;ett  Military  Reserv.uion  in 
Caliiornia  on  a dav  in  November  witli  Santa  Ana  winds  blowing.  Use 
M.ip  Series  V8d5S,  Slteet  1755  1 NW,  1:25000.  Put  the  mortar  position 
on  Hill  15lh  (tibdhO  7 IdOO)  and  the  target  close  to  the  iunct  iori  oi 
two  dirt  ro.ids  .uui  almost  in  the  bed  of  Fria  Creek  (b20dS  74000). 


.3048  STO  7,  9/160  h ST  I 

f P <->  S 

11.932  STO  0,  .000231  STO  1, 

4.682  STO  2,  17.879  STO  3, 

72.914  STO  4,  .0038  STO  5, 

11.037  STO  6,  .7293  STO  7, 

52.51  STO  8,  81.73  STO  9. 

f P <->  S 

RUN.  f W/DATA.  RECORD  BOTH  SIDES. 


10.5  USER  INSTRUCTIONS 


10.  FOUR  DEUCES 


INSTRUCTIONS 


INPUT 

OATA/UNITS 


OUTPUT 

DATA/UNITS 


LOAD 

DATA  AND  PROGRAM  CARDS 

EXAMP 

Ho(m) 

STO  A 

459 

Ro(m) 

STO  B 

4864. E 

Ao()<<) 

STO  C 

4825 

Hj-(m) 

STO  D 

369 

r 

STO  0 

-1 

0 

1- 

STO  1 

-15 

VE  (fp*; 

I STO  2 

-12 

HmW 

STO  3 

1200 

AwO^i) 

STO  4 

2600 

W (mph)  STO  5 

25 

STO  6 

95.7 

13  PRESS  A/7  SECS,  SEE  CO^ECTED  AZ 

14  PRESS  B/7  SECS,j££  CORRECTED  CJ^/^GE' 

15  KEY  IN,  ROUNDED  TO  NEAREST  1/2 

16  PRESS  R/S  {5  SECS),  SEE  CORRECTED  EL 


WARNING.  IF  NEW  RUN  IS  TO  BE  MADE 

WITHOUT  RELOADING  THEN  0 sfo  8. 

ALSO  REGISTERS  B AND  6 MUST  BE  CHECKED. 
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1*11 1 


LL-  A USKR  jV^iOAT  1 ON 


n . I . KKl'KRraCE 

a.  1,.  N.  Pockliam  aiul  R.  W.  Davis,  ,1  Stini'l  ij’it'ii  rvojhifjat  i <.»i  Modt'l 
j\'>‘  I,lUW>‘  Air  Korre  Weapons  Laboratory,  Teeli- 

nleal  Report  AKWl,-TR-7 2-95 , Au^;ust  1972  , ASTIA  No.  A1)902736L. 

..J’ ' 

Tlu'  laser  e(iuatloii  programmed  In  this  section  was  hrounlit  to  my 
attention  by  Lieutenant  Oolonel  R.  S.  Dei.aney,  USAF.  The  ecpiat  um 
applies  to  prcipa>;at  Ion  In  the  atmosphere.  A listing  of  the  variables 
and  parameters  used  In  the  equation  shows  the  factors  consldt'red  In  it. 


S^mbo  1 

Meaji  Ing 

U nits 

F 

power 

wa  1 1 s 

R 

ranne 

km 

1 

avera>;e  intensity 

w.it  t s/cnr 

h 

blockaite  factor 

-- 

K 

thermal  blooming  factor 

— 

a 

atmospheric  ext Inct Ion 

1 / km 

pi'wer  reduction  factor 

— 

'^2 

beamspread  factor 

— 

X 

wave  1 engt  h 

microns 

1) 

diameter  of  primary 
output  mirror 

meters 

^TR 

one  s i>;ma  j liter/ 1 racker 

mlcrorad ians 

«PL 

one  sl>>ma  j liter /plat  form 

micro rad  Ians 

^8L 

one  si^ma  ) 1 tier /boundary 
layer 

mlcrorad Ians 

^vr 

one  slRnia  ) 1 1 ter/atnu'sphere 

mlcrorad ians 

0 

an^le  between  beam  and 
target  normal 

deg 
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U.  >.  KQUATIONS 


100  b • K • ji  _•  oxp  •^_o_os  J9  

kj  • itR-’i(0.9  k.,  X/0)"  + 4 Oj-.j  4-  4 o‘^.)l 


Till'  mimbi'r  100  Is  noi'ili'il  to  tlio  intoiisltv  at  tho  tar^ot  In 

> 

watls/i'Mi"  will'll  R is  ill  k i loiiii't ors . 

Wi'  will  |irogr>im  tliroo  proli loms : 


• Civi'ii  r anil  R,  fiiul  I 

• Oivi'ii  I aiul  R,  fiiul  I' 

• lUvi'ti  1 ami  1’,  fiiul  R. 


ill  of  fort  , till'  proftrani  is  a 0 l>t  i t i zoil  luimnuram. 

Wi'  will  nso  Ni'wti'ii's  nii'llioil  to  not  R nivi'ii  1 aiui  P.  Wo  luivo 


who  I'l' 


1 - l.Po“"'^/R'  , 


1 tHl  ♦ li  • K QOS  0 

k • III  (0.9  K.,  \/0)“  4 .',o‘ 


“ ^R  "pi,  ^’at  • 


Wo  want  till'  root  of 


—1^  L?  ■> 

- M'o  /\C  - \ - 


f IR)  - -(a  4 2/R)  l,Po  “"/ir  , 


f 


lUJi V J I! I ipilfpipp^ 


IPBMIWWP 


-no- 


t 


1 


I 

j 

y 

\ 

j 

4 

i 

I 

I 


i 


The  only  point  of  Interest  is  tlie  determination  of  a good  start - 

-cxR  2 

Inn  value  The  function  e /R  is  ci'iu'ave  upward  and  extremely 

flat  for  even  moderately  larne  R.  If  an  Initial  R^^  Is  picked  that  Is 
greater  than  R,  it  may  well  happen  that  the  flat  tnnnent  projected 
backward  will  generate  negative  values  for  the  successive  R^. 

For  the  root  of  (4), 

R^  - Ae”"*^  . A - hi'/ 1 . 

A first  approximation  is  /a,  hut  this  Is  clearly  too  large.  Uut 
taking 


a A 

R =•  /a  e ^ 
o 


(7) 


a value  smaller  than  the  actual  root  is  found. 

The  program  also  computes  the  Illuminated  area  by 

O.Ol  itR"((0.9  K.,  X/nl"  + 4a"I  sec  0 . (8) 


(1)  The  program  Is  an  example  of  the  use  of  the  flag  K1  to  find 
the  solution  of  any  one  of  three  problems.  The  program  structure  is: 


P 

*fhm.A 

hF?3 

t'.TOl 

Calculate  P 
*f  1.111, 1 
.STt)  A 
hRTN 


R 

*fhllhll 

hF?} 

trro> 

Ca  leu  I.)  to  R 
*fhlll,2 
STC  11 
hR'l'N 


I 

>^11.111.0 
hF?  1 
CTOl 

Calculate  1 
*ri.ni.i 
STC  C 
hRTN 


Flag  F3  Is  sot  by  data  ontry,  and  clearod  by  tost.  If  you  key  P and 
press  A,  F3  is  sot,  P is  stored  in  A,  and  F3  is  cleared.  If  you  then 
key  1 and  press  C,  I is  stored  in  C.  Now  press  B.  Since  F3  is  not 
set,  tlie  step  "GTO  2"  is  skipped  and  R is  calculated. 

(2)  For  the  preceding  problem,  an  "h  PAUSE"  sliows  the  successive 
steps  in  the  convergence  to  two  decimal  places.  When  P is  large 
and/or  I is  small,  tlie  convergence  is  slow  because  the  value  of 
from  (7)  is  small.  You  can  speed  up  convergence  by  choosing  an  R^ 
that  is  large  but  still  smaller  than  the  expected  final  root.  Then 
execute  the  sequence:  Key  R^,  STO  B,  GTO  0,  R/S. 


Example.  Let  the  parameters  be: 


b = 0.9,  K = I,  a = O.l,  k^  = 1.5,  k.,  = 20/9, 


X = 3.6,  0 = 1,  each  of  the  four  sigmas  = 4 (so  tliat 


O Is  64),  0 = 60°.  Then  PRESS  E to  initialise. 


(1)  P = 10  (EEX  6),  PRESS  A,  R = 4,  PRESS  B. 


Now  PRESS  C to  got  1 = 1299.60  watts/cm  . 


PRESS  D to  get  309.47  cm  for  ttie  illuminated  area  at  this  range 


(2)  4,  PRESS  B.  1299.60,  PRESS  C.  PRESS  A to  got 


P = 1 000  002.64  watts  (because  of  roundoff  in  1). 


(3)  10  , PRESS  A,  1299.60,  PRESS  C.  PRESS  B to  get 


R = 4 after  two  iterations  (3.99,  4.00). 


11.6  PROGRAM  USTINC 


A in^E  AA  AA\M]iLK) 


IJ-A .._  kkkerknce 

a.  E.  W.  Paxson,  l\wt  ujl  lu  P:  i\‘i‘irit  natoru  Mortar  yiro,  Tho 
Rand  Corporation,  P-58b7,  February  1977. 

DISCUSSION 

In  manual  war  gaming,  people  instead  of  computers  make  tlie  tac- 
tical decisions  depending  on  tl»e  situation.  Yet  there  are  recurring 
events,  such  as  firefights  or  air  intercepts,  whose  outcimies  must 
be  assessed  systematically  on  tho  basis  of  agreed-upon  rules  and 
planning  factors.  The  assessment  consists  of  taking  a sample  from 
a probability  distribution  function  for  outcomes,  since  the  game  cannot 
move  on  without  a definite  result.  War  gamers  call  this  "shaking  the 
dice."  It  will  not  do  to  say,  "The  probability  is  0.4  that  at  least 
three  tanks  of  a company  of  10  will  be  destroyed."  The  game  demands 
a statement  such  as,  "In  this  firefiglit  3 tanks  were  immobilized  and 
one  w.is  set  afire." 

* In  any  game,  there  are  usually  large  numbers  of  recurring  events 
of  the  same  type.  One  expects  tliat  the  results  for  this  set  of  events 
will  average  out — giving  the  mean  beliavior  of  the  model  underlying  the 
event  type.  To  clarify  this  statement,  t l»e  simple  f und.imenta  1 principl 
of  the  Monte  Carlo  method  (sampling  from  a probability  distribution) 
is  invoked.  For  example,  suppose  an  event  can  iiave  four  outcomes 
0^,  0.,,  0^,  0^  with  respective  probabilities  Pj,  p^,  p^,  p^,  where 

p^  + P,  + P3  + P4  = 1.  Put  Pj  = p^,  P,  = p^  + p^,  P3  = Pj  -F  p,  + P3, 

P^  = p^  + p,,  + P3  P_^  “ Take  a large  number  M of  random  numbers 

uniformly  distributed  over  the  interval  (0,  1).  By  the  Law  of  Large 
Numbers,  one  expects,  then,  that  of  tlie  M random  numbers, 

PjM  will  be  in  the  interval  0 to  Pj, 
p,M  will  be  in  the  interval  P^  to  P., , 

PjM  will  be  in  the  interval  P.,  to  P^, 


P.M  will  be  in  the  Interval  P,  to  1. 


This  section  shows  by  an  example  lu'w  liaiul  calculator  programs 
can  provide  assessments  of  this  nature  to  speed  np  the  plav  ot  manual 
war  games. 

A g.ime  for  series  of  games)  is  set  np  to  test  the  beh.ivior  in 
a full  tactical  environment  of  an  innovative  weapon  svsti'm,  which 
is  to  destroy  enemy  armor  as  part  of  a ci’mbined  arms  ti'rce. 

I'he  proposed  svstem  envisages  a new  tvpe  of  mortar  nnind  that 
has  a heat-seeking  sensor  head  controlling  the  maneuver  of  the  round 
to  a target  during  the  steep  terminal  phase  of  the  trajectorv.  These 
rounds  are  ripple-fired  at  a set  of  armored  targets,  picking  targets 
at  random.  The  sensor  head  will  reject  targets  previously  set  afire 
or  exploded  iK-kills)  to  avoid  the  mot h-and-f 1 ame  effect.  But  pre- 
viously immobilised  vehicles  (.M-kills)  m.iy  wastefully  absorb  rounds 
that  home  on  the  still-warm  engines.  Such  overkill  is  a common  battle- 
field event. 

If  at  .iny  time  t during  the  eng.igement  there  are  j vehicles  still 
moving  and  k vehicles  either  still  moving  or  immobilized  by  previous 
fire,  then  the  probability  Pfj,  k,  l)  of  the  state  fj.  kl  at  time  t 
can  be  determined  analytically  (see  Ref.  .il.  But  the  .inalvsis  is 
complex,  as  are  the  resulting  formulas  for  l'(i.  k.  1 1 . For  use  in 
war  gaming,  it  would  still  be  necessary  to  s.imple  with  respect  to 
P(|,  k,  t).  It  is  much  more  economical  to  adopt  the  pri'cedure  of  the 
following  subsection. 

12.3.  FQFATIONS 

Let  r be  the  probability  that  a round  gets  an  M-kill,  immobiliz- 
ing the  target.  Let  s be  the  probab^ity  of  a R-kill,  with  the  tar- 
get exploded  or  set  afire.  Initially,  there  are  moving  targets, 
against  which  N rounds  can  be  fired  during  the  target  exposure  interval. 

Then  the  state  changes  per  round,  with  their  assoc  i.tted  proba- 
bilities, are  shown  in  the  following  table,  rt'member ing  that  of  k 
targets  one  is  picked  at  random: 
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State  Change  I’robabilitv  Reason 


The  corresponding  bounds  in  the  interval  (0,  1)  are 
Pj  = 1 - s - rj/k,  P.,  = 1 - s,  P'3  = 1 - s + s,i/k.  P4  = 1- 

Pseudorandom  numbers  over  the  interval  (0,  1)  will  be  generiited 
bv  the  multiplicative  Linear  congruent  ial  method.  Let  R^^  be  the 
initial  "seed,"  a seven-digit  decimal  fraction.  Let  m be  the  multi- 
plier, and  let  R^^  be  the  ith  pseudorandom  number  generated.  Then 

= fractional  part  of  (mR^)  . 

In  the  Hewlett-Packard  Stat  Pac  1 (section  04),  the  values 

R = 0.5284163  and  m = 997  are  used.  Witli  these  choices,  the  number 
o 

of  different  random  numbers  before  returning  to  R^^  (the  period  of  the 

sequence)  is  500,000.  The  sequence  passes  ttie  standard  tests  for 
* 

r.indomness . 

12.4.  PROGRAM  NOTES 

Figure  12.1  provides  the  logical  flow  of  the  programming. 

It  is  essential  to  save  the  last  random  number,  to  be  used  as 
the  "seed"  for  the  next  evaluation. 

*0ther  "good"  pairs  are:  R = O.lllllll,  m = 291;  R = 0.7742713, 
m = 997. 


■I 

A 


12.6  USIR  INSTRUCTIONS 


SHAKING  THL  DICE 


INSTNUCnONS 


LOAD  PROGRAM 

STO  .5284163  OR  OTHER  SEED  IN  Rq 
STO  9V7  OR  OTHER  MULTtPL'ER  IN  R, 
STO  r IN  R^ 

STO  s IN  R3 

INITIAL  TARGETS  A,  ENTER 

ROUNDS  FIRED  N,  KEY,  PRESS  A 
(NO  ROUNDS  LEFT  .1 
2 MOVING  TARGETS, 

1  M-KILL,  2 K- KILLS 

RECORD  OR  SAVt  NEW  SEED, 

(HERE  NEW  SEED  IS  IN  R,,.! 

NEXT  RUN  WITH  SAME  VALUES 

0 MOVING  TARGETS 

2 M- KILLS,  3 K- KILLS 
RECORD  OR  SAVE. 

NEXT  RUN,  WITH  NEW  SEED  IN  R^ 

3 MOVING  TARGETS 

NO  M- KILLS,  2 K- KILLS. 
RECORD. 


(NPUT 

OATA,UNITS 

KKVS 

OUTPUT 
DATA  UNITS 

EXAMPLES 

1 M 1 

1 II  1 

1 II  1 

1 II  1 

1 II  1 

1 II  1 

.2 

ISTOII  2 1 

0.20 

1 II  1 

.3 

ISTOll  3 1 

0.30 

1 II  1 

5 

1 1 1 r NT  1 

5. IX) 

1 1 1 1 

6 

1 II  A ! 

0. 

1 II  1 

! II  1 

2. 

1 II  1 

1 II  1 

3. 

1 II  1 

.7-‘)OOEI27 

5 

1 IIENTI 

0 

1 II  1 

6 

1 1 1 A 1 

0 

1 II  1 

1 II  1 

1 II  1 

1 II  1 

,7336tf8o 

1 II  1 

5 

1 IIENTl 

0 

1 II  1 

6 

1 1 1 A I 

3 

1 N 1 

1 II  1 

1 II  1 

3 

1 II  1 

8173707 

12.6  SHAKING  THE  DICE 


KtY  tWTWY  KtV  COOf 
Oi.  ttf'..') 

ee:  srci  ;5  ,c. 

.s»i 

6C-<  tliih  ;i 

oo;  src.~ 


COMMfNTS  STFP 


t^BiS 

^ ' 

kClB 

.c* 

••i.' 

br:5 

I’f 

i 

T • 

C 4 

kCL*i 

, « 

Tc 

■t- 

j.: 

sroi 

■'r 

• J 

A.'...'' 

\ 

. 7^ 

- 

r 

•»  « 

» C 

C > 

- 

. r 
1% 

$104 

IH 

1 

T t 

kii! 

t 

- 

•#  J 

STi;? 

■'r 

jr 

A.  k. 

7c 

1 j 

kCL3 

v'  C 

1 0 

, **r 

* 

rr 

s:0B 

■»r 

jt 

l^:.4 

or 

icCl 

J.* 

.\»r 

»6- 

, •jr 

bUV 

C w 

Rc:r 

tJ 

kCii 

.V>V’ 

c • 

. -jc 

il03 

JJ" 

rj 

RCic 

ROlO 

Of 

\iV' 

• "?*■ 

0104 

i-i 

RUB 

Jo 

i » 

1 

i'i 

- 

-J5 

$10$ 

J' 

,*  ^ 

B$S1 

;t  :r 

•*0 

6 10$ 

t J 

iio:^ 

05 

iLti.r 

o: 

RCLO 

Jo 

lie 

RCLl 

Jo 

r; 

V 

- 7*^ 

FSi 

SlOO 

■^r 

Of 

NEXT  RANDOM  NR 
TESTS  FOR 
STATE  CHANGES 


STro 

KkY  tNTHY 

KtY  coni 

COMMf  N1S 

es: 

ps:j 

-c'  J5  'ic 

KJO  TMAMrXP 

ess 

erne 

::  i: 

es.'S 

eros 

JJ  0! 

iM'i 

PSfl 

4i$i.$ 

:;  os 

Pfj 

RCLO 

Jt  1 4 

eo: 

J 

0! 

6t$ 

4f 

Ot>4 

STSS 

“t  . « 
i * 

j-’ 

ess 

DSSI 

io  i; 

Ooo 

ST  OS 

j; 

es.' 

bTi'S 

^ * *r 

ess 

»lF4.  ■» 

- . C-t 

et.? 

Jc  . I 

o’r 

ere 

1 

0* 

e'l 

_ . r 

•1  • 

ers 

STSh 

7 ^ 

j-' 

arj 

R.'se 

J; 

074 

1 

01 

ors 

- 

■4$ 

07$ 

s;jt 

^ i *. 

k-i 

077 

$$:i 

-i  •ic 

07$ 

$10$ 

vJ 

079 

*LB1.S 

I'S 

C>AO 

ese 

PSFO 

-c'J  CO 

esi 

Rl'l.  1 

36  4t 

ess 
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1 ) . OPT  IMIW  OJK  RESqilRCKS 


t 


D.l.  RKKERENCKS 


a.  Italian  Luss  and  Shiv  K.  (aipta,  "Allocation  of  Eftort  Resources 

/\inonn  Competing  Activities,"  lu'th'tnh’li,  Vol.  23, 

No.  2,  M.irch-Apr  1 1 l't75. 

b.  A.  Charnes  and  W.  W.  Cooper,  "The  Theory  of  Search:  Optimum 
Dlstr  Ihut  ions  of  Searcli  Effort,"  Ma>Uhji’.rr,,'ni  t’liot' , Vol.  S, 


I 3.2.  DISCIIS.SION 

Tliis  is  a topic  in  nonlinear  (convex)  programming.  3'he  bibliog- 
raphy of  Ref.  b indicates  optimum  search  as  one  milltarv  motlv.ation 
for  the  study  of  programming  of  this  nature.  Another  milltarv  appli- 
cation asks  for  the  optimum  allocation  of  weapons  to  a target  system 
organized  into  classes  of  targets  of  given  number  and  value  (or  wliich 
the  weapons'  kill  probabilities  difi'er.  tMvilian  appl  Icat  iiuis  arise 
in  al  loc.it  ing  advertising  budgets,  portfolio  selection,  .ind  budgeting 
(Kef.  a). 

The  problem  is  formulated  as  follows.  Let  11  be  t lu’  .iv.i  i l.ible 

resources  to  be  allocated  to  N activities  in  the  amounts  x,,  x.,  ...,  x„ , 

1 Z N 

where 


xj  a^  S:  0 .*  (1) 

If  Xj  is  applied  to  activity  I,  the  return  on  the  investment  is 
Qj(Xj),  which  is  a differentiable  and  strictly  concave  Increasing 
function.  I'or  example,  in  the  weapon  al.'locat  ion  appl  ic.it  ion,  if  there 
iire  Tj  targets  in  class  1 all  of  value  Vj,  .ind  if  the  SSl’K  is  1’^,  then 


+ X.,  + 


+ X, 


r ■'‘l^'l 

Q,(Xj)  - TjVj  1 - (1  - P^)  , 


* 

Higher  ;uithorlty  or  other  considerations  m.iv  dlct.ite  that  s^^me 
act Ivlt-ies  he  assigned  minimum  (nonzero)  resources. 


! 


-123- 


stnot'  on  the  averano  weapons  are  applied  to  each  target  and  the 

square  braekets  contain  the  kill  probability  per  target.  Figure  13.1 
shows  [1  - (1  - 0.!))  ).  Note  the  rapid  decrease  In  marginal  re- 

turns tor  the  larger  values  of  x.  One  "standard"  form  for  Qj(Xj)  Is 

Qj(Xj)  - .S^ll  - exp  (-m^x^)l  , (3) 


which  yields  (2)  If  and  m^  - (I  - 1’^). 

We  now  want  to  m.ixlmlze 


R » Q,(Xj)  + 


subject  to  (1). 

11.1.  KQIIATIONS 
From  ( 1) , 


9Q^(Xj)/3xj  - S^m^  exp  (-m^x^' 


Relndex  tl\e  activities  so  that  S.m,  S,  . , m,.,.  wlilch  are  the  marginal 

I 1 IT  1 1 + 1 

returns  at  the  origin.  Introduce  the  l.agranglan  multipliers  and 
max imlze 


I'li'N'  - "t  • ( g - »]| 


To  tills  eiul 


9Q,(x  ) 9Q  (x.) 

• a M m — -- 

3x.  t 3x. 


t,  1 - 1 i 


express  the  marginal  equilibrium.  We  find  that  the  successive 
.ire  ci'iinected  hv 


?S?)5 
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t+1 

**  E »/'"(  “ («e+,  "V+i 


• Mj  1/m.  . U>) 


wlu'ro  "**"  mi'.iiis  "to  t ho  pi'wor"  .inii  ” I'l”']  (-hm^). 

Tho  solvit  lou  alni'rlllim  is  t vv  sii'p  at  tlial  d tor  v»'liloli 


^e+i  "V-H 


'd+l 


t71 


riu'ii 


1 

m , 


S,m, 


i - I 


aiul  tho  maximum  rot  urn  is 


C 

E 

i-i 


■"i'’'i  - ."d 
m . 

i 


I'M 


rir.uri'  11.2  slu'ws  tho  phvsiral  roa  I i .*..il  i vMi  v'l  this  a 1 r.or  i t lim . 

Tho  (1  th.it  p.ivi's  t lu'  m.iximum  roturii  is  nsoil  I irst,  thou  tho  v,'  th.it 

A 

vloUis  tho  iioxl-lil>>liost  rotiirii,  auvl  sv'  on.  Kiiui  tho  x.  whoro  tlio 

1 

slopi's  of  tho  v^'-oiirvos  afv-  all  ovpial,  siiu'o  L.iklun  ,i  unit  of  n'Siviiroos 
Irv’iii  i>iu'  pooki't  .iiul  putting;  it  iutv'  ain't  hor  Um'vinp,  away  I rom  tho 
mar>;iual  v'qii  i I ihr  ium)  ilov'roasos  ovor.ill  roturus.  tlU'.irlv.  slop  wluni 
ovou  till’  m.ir>;iii.il  rv’turn  .it  Llii'  oripiii  ot  .in  aoiiviiv  o.iiiiiot  ov'Ut  r ihiito 
.IS  iiuioli  .IS  its  moro  Uu'r.ilivo  foili'ws. 

,1 1 is  oharao  t or  i St  i V'  of  siiv'h  solutii'iis  that  soiiu'  aotivitii'S  tv'- 

I 

ooivi'.no  ol  fort  . for  o.xamplv',  in  sv'mo  air  vlofv'iiso  prohloms,  i.ir^tots 
ol  vi'rv  low  valiio  ^ot  no  liofoiiso  a 1 li'o.il  ion . 

A oi’tis  ivlorovi  I’l'  oxporionooil  giioss  o.iii  f ri'vnioiil  ly  oi'tiio  within  h 
pt'-roont  v'f  till'  v-aloiilatoil  optimum.  This  is  ap.aiu  typioal  of  t hoso 
prv'hloms.  Hut  (11  you  iii'V’v'r  kiu'w  how  olosv'  yiui  art'  to  t hi'  o|'t  imum, 

(2)  tho  ImJ'rovomoul  l'oin>’  iulorost  ou  an  luvostmout,  iiiakliu;  many  iu- 
vostmoiils  /ulits  up,  auil  (11  if,  for  roasous  oxtorual  to  tho  moilol  . 
oortalu  a'ooalions  to  oi'ri.iin  oat  ogor  ios  .iro  spootftoil,  you  will 
know  what  poually  is  palii. 
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Kxamp 1 e . Basod  on  ostlmatod  initial  position  and  drift,  a damaged 
ship  may  be  In  one  of  fi>ur  search  sectors  with  estimated  probabilities 
0.5,  0.3,  0.15,  0.05.  Because  of  weather  and  sea  conditions,  the  con- 
ditional probabilities  that  one  search  sortie  will  detect  the  target 
in  these  sectors  are  0.1,  0.1,  0,4,  0.2,  respectively.  Find  tlie  best 
allocation  of  20  sorties  to  maximize  the  detection  probability. 

We  have 


Sec  tor 

s 

r. 

m 

sm 

Re  index 

1 

0.5 

0.1 

0.  105 

0.053 

3 

•> 

0.3 

0.3 

0.357 

0.107 

1 

1 

0.15 

0.4 

0.511 

0.077 

2 

4 

0.05 

0,2 

0,223 

0.011 

4 

* 

* 

is 

A 

1 yields 

xi  = 5. 

.4,  x^ 

= 3.1, 

II 

,5,  x^  = 0 

is  20  as  it  must  be.  The  detection  probability  is  0.73.  Figure  13.3 
sliows  the  Q functions  and  tlie  solution  points  and  illustrates  tlie 
algorithm.  But  it  indicates  also  that  the  schematic  of  Fig.  13.2  is 


not  the  general  case.  That  is,  it  is  not  necessarily  Irvie  iltat 
* * * 

X,  i X , S X.  . . . . 


3 • • • • 


I 3 . 4 . _ PROtllOUl  N0Ti;S 

Although  sliort  and  simple  in  structure,  the  progr.am  has  several 
elements  of  interest. 


(1)  I.Bl.  A loads  l/m^  using  indirect  storing.  l.BL  B lo.ids  S.m^ 


into  protected  secondary  storage,  also  using  indirect  storing.  But 
c.ire  must  be  taken  to  indicate  wliether  f P -*-*■  S puts  the  memory  in 
the  primary  or  secondary  storage  mode.  This  is  particularly  important 


when  exiting  from  .i  loop  in  the  middle  of  a label. 

* 


(2)  I.Bl.  I determines  the  successive  x^  (Kq.  (8))  by  f 1S2  up 


to  1. 

* 

(1)  I.BL  2 llu'ii  detetmlnes  R by  forming  the  sum  (F.q.  (6>)  In 
reverse  order  from  i down  to  1.  f DSZ  provides  the  I'xit  by  skipping 
on  zero. 


» tt£L.  ' - « ' itf.li— » ' 
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M,=  .0158 

Q // 

0.2  - / / 


/■Ms 


2 4 6 8 10  12  14 

Number  of  sorties 

Fig.  13.3 — A search  problem 
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13.5  USER  INSTRUCTIONS 


13.  OPTIMUM  ALLOCATION  OF  RESOURCES 


quality  HUCnCAfilii 

COi  Y i UKMSHSD  TO  DDC 


13.6  OPTIMUM  ALLOCATION  OF  RESOURCES 


KCV  INTfiV 

CiJI  *i,6Lh 

oc:  : X 


ZllU I NDIRECT  STORAGE 
. 5i  OF  1 m; 

STO  (i) 


SINCE  10  IN  Rj, 
INDIRECT  STORAGE 
OF  S;m;  IN  11,12,... 


h ST  I (RE-INDEX) 
1/m, 


KEY  coot 
JS  13 
51 . 


M , (6) 

Sj  + j t 1 - M^  t ,(7) 

PRI 

i 

I 

1 

NEW  PARTIAL  SUM 


i i + 1 
I 

1 1 

( ) , m i + 1 ) 
SEC 


M i M (6)r 
SWITCH  ' 
LOOP 


Cc':' 

-24 

I’oT 

LN 

3Z 

Pis 

16-51 

Os? 

pql; 

36  45 

b^ 

ore 

X 

-35 

P s 

51 

c 

p:s 

16-51 

f'r 

f:li 

36  46 

07- 

PCLC 

36  13 

L’'r 

16-33 

GTOf 

22  »2 

c' 

or  01 

22  Bi 

t . . 

iSU 

21  92 

PCL ! 

36  45 

oai’ 

~ ‘»r  - 

w 

RCL  A 

36  11 

e;': 

- 

-45 

p;rs 

16-51 

RCL ! 

36  45 

• 

i L ■ 

X 

-35 

' cjr 

ST*0 

35 

-55  66 

li’:' 

16-51 

ft” 

cs:i 

16 

25  46 

iv: 

C-T02 

22  92 

' Of.' 

o^S 

16-51 

bi^ 

PCL6 

36  66 

? 3 * 

RTN 

24 

r - - -- 

LEFT  030  IN  SEC 
Sim; 


PRI 

RCL ( i ) , 1 m j 
Xi*  (8) 

SEC 


FINISHED  ? 

LOOP 

LEFT  076  IN  SEC 
RCL(i)  , Sj  m i 


LOOP  FOR  SUM 
PRI 


Sfmj  Km, 


LABELS 

‘’(6),  TEST K 


- T - 

B . 


!" 

" 1 my 

Sg  m 
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EXAMPLE.  6 = 200  w«apon$  are  fo  be  allocated  to  150  target! 


in  5 

value  classes 

to  get 

maximum 

return . 

Ti 

''i 

Pi 

Gloss 

Number 

Value 

SSPK 

S 

m 

Sm 

Re-index 

1 

10 

5 

.05 

50 

.00513 

.25647 

5 

2 

20 

4 

.2 

80 

.01116 

.89257 

3 

3 

30 

3 

.3 

90 

.01189 

1.07002 

1 

4 

40 

2 

.4 

80 

.01277 

1.02165 

2 

5 

50 

1 

.5 

50 

.01386 

.69315 

4 

The  formulas  are;  Q;  = 5;  [l  — (l-Pj)**'  , 


Si=Ti  vj,  mi=-in(l-Pi),Ti. 


3 

KEY  IN  m,  PRESS  A 

.01189 

1 1 1 A 1 

84.10 

.01277 

1 1 1 A 1 

78.31 

.01116 

1 1 1 A 1 

89.61 

.01386 

1 11*  A 1 

72.15 

.00513 

1 1 1 A 1 

194.93 

4 

10  h STI 

10 

1 h 1 1 ST  1 1 

10.00 

5 

KEY  IN  S,  PRESS  B 

1.07002 

1 1 1 B j 

1.07 

1.02165 

1 11  B 1 

1.02 

.89257 

1 1 1 B 1 

0.89 

.69315 

1 1 1 6 1 

0.69 

.25647 

1 11  B 1 

0.26 

6 

200  STO  B 

200 

ISTOll  B 1 

200.00 

DSP  0 

1 II  1 

200. 

7 

PRESS  C 

1 11  c 1 

4. 

(NO  WEAPONS  APPLIED  TO  CLASS  1) 

1 II  1 

8 

R,S 

1 1 1 R..^  1 

65. 

R,5 

1 1 1 R/s  1 

57. 

R/S 

1 1 1 R.-S  1 

53. 

. 

1 l(R.^I 

25 

(SUM  IS  200) 

1 II  1 

9 

VS  . 

1 1 1 R,  is  1 

140. 

(THE  VALUE  ACHIEVED  IS  140) 

1 II  1 

1 II  1 

1 It  1 

1 II  1 

1 11  1 

1 II  1 

1 II  1 

1 II  1 

1 U 1 

j 

i 


i 

♦i' 


j 

h 

.1 
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lA.  LOC-LINEAR  CUMULATIVE  AVERAGE  AND  UNIT  COSTING 


14.1.  REFERENCES 


H.  E.  Boren,  Jr.  and  H.  G.  Campbell,  Learning  Curve  Tables, 
The  Rand  Corporation,  RM-6191-PR,  April  1970  (3  vols.). 

R.  W.  Hamming,  Numerical  Methods  for  Scientists  and  Engineer. 
McG:raw-Hill,  New  York,  1962. 


14.2;  DISCUSSION 


Learning  curve  theory  assumes  that  each  time  the  total  quantity 
of  items  produced  doubles,  the  cost  per  item  is  reduced  to  a constant 
percentage  of  the  previous  cost.  The  relationship  is  given  by  the 
power  or  log-linear  relation 


where  x is  the  cumulative  pruduction  quantity.  If  y is  the  average 
cost  of  the  first  x units,  we  have  the  cumulative  average  learning 
curve.  If  y is  the  cost  of  the  Kth  unit,  we  have  the  unit  learning 


For  plotting  purposes,  the  midpoint  x^ 
lot  average  cost,  is  ti  be  determined. 


14.3.  EQUATIONS 


If  S is  the  fraction  to  which  cost  decreases  (the  learning  curve 
percentage)  when  the  quantity  is  doubled,  the  slope  of  the  learning 


Then,  with  a the  cost  of  the  first  unit 


t 


is  till'  aviM’ano  I'osl  of  tiu'  f irst  x units,  ami  llio  total  cost  lor  tlu‘ 
I Irst  X units  is 


XV 


Ivt  1 


( I 


Thi'  unit  oost  at  tlu'  \th  unit  is 

r l>+  I , , , l'+  I 1 

Vu  ■ - (X  - I)  J . 


(4 


l.ot  (x  . V ) lu'  t lu'  miiipoinl  tor  t ho  first  ItU  of  n units.  'I’lu' 
"’h 

V ■ V • .'in  , and  iisiiu;  (•(),  x is  I In'  si'lullon  I'f  Iho  o(|uat  ion 

l\l  O III 


ht  I , I' 

t - (x  - 1)  - n 

III  III 


(') 


Tlio  .iliovo  oi|ual  ions  apiilv  to  t lio  lo^^-lliUMf  oiimul.'K  Ivo  o.iso. 

for  tlio  lop.-lliu'ar  unit  ourvo,  I lio  unit  oust  at  t ho  \tli  unit  is 


h 

V ax 
u 

.uul  t ho  oumulatlvi'  avorap.o  oust  for  tho  first  n units  is 


V 

o 


li 

X 


((< 


(/ 


slnoo  tho  total  oust  is 


T - a 


n 


h 

X 


(« 


TIu’  nildnolnt  (x  , v ) Is  ilot  c'rnilnod  liv 
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Tho  only  programming  probloms  arc>  posod  by  Kqs.  (S)  and  (8). 

In  applying  Newton's  method  to  (5),  a good  first  estimate  is 

needed  for  x . We  have 
m 


f(x) 


b+1 

X 


(x 


1) 


b+1 


b 

n 


b+1 

= X 


b+1 

X 


(1  - l/x)*”^^ 


b 

n 


and  using  the  first  two  terms  In  the  expansion  of  (1  - 1/x)  , the 

first  estimate  for  x is 


Xjj  = n/(b  + 1)'^’’  . (10) 

which  is  very  close  for  large  n. 

To  get  an  excellent  approximation  for  the  sum  In  Kqs.  (7)  and 
(8),  use  the  Cregory  formula  (Ref.  a,  p.  158): 

n ^ 

f f(x)  dx  = .,  ' i ^ ^^*0  ~ ^*n-l^ 

0 

-■h  * ■■■  ■ 

Because  x*’  is  integrable,  the  formula  can  be  applied  in  tlie  backward 
direction  to  get  the  sum.  Since  x^’  is  steep  for  small  x,  we  start 
the  sum  at  x = 4.  (hamming  calls  this  "low  cunning.")  We  have 

n+1 

X:  X*’  = I + 2*’  + 1^’  + f x‘’’dx  - i [l*’  + (n  + 1)'’] 
x=l 

1 f.b  ,b  . i\8  , hi 

- j-y  14  - 1 - (n  + 1)  + 

+ 1^5*’  -2*4*’+  l’’  + (n  + 1)'’  - 2 • n’  + (n  - o’’]  , 
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yioldinf;  tlio  l iitlior  liio  U'n.inl  ri'sult 


x*’  - I + 2*’  + 

x«l 


3 2— 

8 " h + 1 


>*’  - ' s'’ 

(1  2'» 


(11) 


+ — , (n  - 1)  - , « + 

2*4  (i 


n + 1 3 . I \ 1’ 

ir-TT'H  • 


PRDCRvlM  NOTKS 


1.  In  I'iiulinn  a flrsl  aiiprox  imal  ion  to  tho  root  of  lai . (S)  hv 
oxprossion  (111),  a valno  loss  ll\an  1 arisos  tor  small  n ('•'•>),  I'liis 
wvniUl  load  to  an  llrror  sijtnal  lu-oanso  of  (x  - 1)'’^'.  I.IU.  2 iirovidos 
a start  inj;  valno  of  l.Ol  to  avoid  tliis. 

2.  I'lir  small  valnos  of  n (-'S'),  j;i-t  ) x'’  bv  mamial  oalonlation 
(not  proy,rammod ) . 

3.  Tbo  |>roy,ramm  i nj;  ilono  lion'  involvod  oxpor  imont  iny  with  tlu'  rola- 
tivoly  nnoommon  soooiul-ordor  Nowton  mothod  for  tho  root  ol  f(x)  “ 0. 

A 

I'ho  formnla  is 


n+1 


t (x,^) 

’‘n  ) 

u 


\ + 


r(x  ) r^x  ) 

n \\ 

2(f'(x,))' 


I'omparod  witii  tlio  first-ordor  nu'l  hod's 


’'n+1  ’'n  f^(x  ) 

n 


la'nvory.i’noi'  to  a y.ivon  .loi-nraoy  is  somowiiat  fast  or  (20  poroi'iit  in 
rnnniny,  t imo)  hnt  at  additional  ptM>;ramm  inp,  v-ost  to  pot  l'(x).  TIu' 
t i rst -iirdi'f  n.othoil  is  nsod,  omplovinp  OSi’  2,  I KNO,  to  pot  (wo-dooimal- 
pl.u'i'  ai'i'iirai'v. 

•4.  Ohvionslv,  a «•  1 in  propramminp,  sinoo  a is  simpiv  a multiplii'r 


Thoro  is  a t vpoprai'li  i oa  I orror  in  this  tinmnla  as  pivi'n  on  p.  82 
ot  Kof.  a. 


14.5  USER  INSTRUCTIONS 


^114.  LOG- LINEAR  CUMULATIVE  AVERAGE  AND  UNIT  COSTINGS^ 


INSTRUCTIONS 

INPUT 

OATA/UMTS 

KEY  IN  % AS  .XX,  PRESS  A 

.69 

KEY  IN  n,  PRESS  R.'S 

81 

FOR  CUM.  AVG.,  PRESS  B. 

SEE  CUM.  AVG.  COST 

PRESS  R S.  SEE  CUM.  TOT.  COST. 

PRESS  R 'S.  SEE  COST  OF  nTH  UNIT 

PRESS  C.  SEE  MIDPOINT  Xn, 

FOR  UNIT  COSTING,  REPEAT  1 AND  2 

ABOVE,  DSP  3. 

.69 

81 

PRESS  D.  SEE  COST  OF  nTH  UNIT. 

PRESS  R/S.  SEE  TOTAL  COST. 

PRESS  R S.  SEE  AVG  COST. 

PRESS  R 'S.  SEE  MIDPOINT  X.„. 

OUTPUT 

OATAUNITS 


7.706 

0.044 

19.85 


0.690 

0.465 

0.095 

15.021 

0.185 

23.281 


4 


h 


?4.6  PROGRAM  USTING 


15.  TIMK-PHASED  PROCUREMENT  COSTING 


15.1.  REFERENCES 


None . 


15.2.  DISCUSSION 

The  costing  approach  of  this  section  is  a hitherto  undocumented 
model  owing  to  H.'G.  Massey  of  The  Rand  Corporation.  For  definite- 
ness, the  model  Will  be  explained  in  terms  of  aircraft  procurement, 
although  it  applies  equally  to'  the  procurement  of  any  system  with 
any  number  of  components.  The  initial  batcli  of  test  articles  is  not 
priced,  it  being  assumed  that  they  are  funded  from  another  account. 
But  the  average  prices  of  the  batch's  components  provide  the  costs 
of  initial  articles  for  learning  curve  calculations.  Overruns  and 
inflation  are  not  included.  The  model  determines,  by  fiscal  year. 
New  Obi igat ional 'Authority  (NOA)  required  to  meet  the  production/ 
delivery  schedule. 

15.3.  EQUATION.S 

Specify  that  two  required  sequences  of  possessed  aircraft 

n = I,  ...,  N,  are  to  be  in  the  fleet  at  the  end  of  year  n,  wliere 

S are  squadron  or  UE  aircraft  and  T are  training  aircraft, 
n n 

If  A and  B arc  flying  hours  per  year  per  aircraft  (Fll/Y),  for 
squadron  and  training  aircraft,  the  iVirnilati lU'  fleet  flying  hours 
tlirough  year  n are  approximately 


H = V [AS,  + BT. 

n 4-'  i i 
1=1 


The  ■.‘unulatLvc  fleet  attrition  Is  given  by 


a = C • H 
n n 
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A fraotlon  A ot  tho  fleet  Is  assluned  to  eoniniand  support  (pipeline). 
With  allowance  for  these  twi'  factors,  the  ^w:ui'<iti th^  number  of  air- 
craft to  be  p-ocured  tlirougb.  year  n is 

Q - (I  + A)  P + a . (3) 
'n  n n 


If  the  fleet  Is  to  be  kept  in  a steadv-state  condition  (S  , T ) 

N N 

fi>r  M more  years,  aircraft  will  be  delivered  in  year  N + 1 to  meet  the 

attrition  requirements  of  these  M years.  Hence  H„ . , and  a^, . , will  be 

N+ I N+ 1 

calculated  to  determine  this  final  buy. 

An  aircraft  lias  three  major  components:  engines,  airframe,  and 
avionics.  There  is  a cumulative  buy  program  f.'r  each  component,  allow- 
ing for  lead  times  and  learning  curve  effects,  (aimulative  average 
costing  is  most  convenient  (see  Sec.  14). 

Take  engines  first.  The  procurement  parameters  can  be  written 
as  the  string  of  numbers 


I'lj’  ^11’  *^12’  ^12’  ^'l3’  *'!’  ^1^  ’ 

where  the  first  subscript  refers  to  component  I (engines)  .ind  a^  is 
the  cost  of  the  first  article,  p^^j  Is  the  learning  curve  percentage 
(or  rather  fraction)  up  to  article  number  p^  is  the  learning 

curve  percentage  for  articles  Xp  + I up  to  x^.,,  and  p^  ^ is  the  sub- 
sequent percentage.  Of  course,  a learning  curve  may  have  only  one 
or  two  segments.  The  lead  time  In  months  Is  tj,  so  that  an  engine 
takes  tj  months  from  start  of  fabrication  to  delivery  to  the  fleet 
as  part  of  a complete  aircraft.  (Mating  of  engines  and  avionics  with 
the  airfriune  Is  taken  to  he  part  of  the  airframe  lead  time  and  cost.) 
The  number  of  engines  required  per  aircraft  is  Up  allowing  for  multi- 
engine  models  and  spares.  With  these  production  parameters: 

I'p  to  Xp  the  cumulative  average  total  cost  for  x articles  is 


(A) 
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whero  - fi  -• 

From  Xjj  + 1 to  the  cumulative  total  cost  is 


ai  • Xii  • X 


and  above  ^ 


'*l  ■ ’'ll 


^2-^3  ^’13-"^ 


These  expressions  are  verified  by  putting  x «=  Xj^^  in  (3)  and  x = x^,^ 
in  (b) . These  expressions  define  a c’tiirtiilatLt^c  cost  function  Cj(x) 
for  all  articles  numbered  1 through  x. 

Next  express  t in  years  and  write 


t|  - INT(t^)  + FRAC(tj)  . 


Thus  a 2b-month  lead  time  is  2.167  = 2 + 0.167.  Then  fabrication 
starts  at  -tj^,  but  during  the  fiscal  year  -lNT(tj),  only  pQj  FRAC(tj) 
articles  will  be  started,  requiring  NOA  of 


CjCuQj  FRAC(t^))  . 


To  the  end  of  tlie  next  FY,  -INT(tj^)  + 1,  the  i-tirmlutive  number  of 
engines  started  is 


+ m(Q2  - Q^)  FRAC(tj)) 


and  Nl)A  for  that  FY  is 


Cj(hQj  + h(Q2  - Qj)  FRAC(t^)) 


By  example,  if  tj^  • 26  months,  - 16.11,  ” 35.17, 
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.. 


PROD  FOP 
Ql  a/c 


PROD  FOR 
% - a/c 


= 16.91  = 35.17 

(7  and  50  are  now  rounded  value.s) 


Note  that  deliveries  to  Che  fleet  are  made  at  a uniform  rate  during 
a year.  So  that  formula ’(7)  will  apply  for  the  first  year  and  the 
last  year,  store  = 0 and  store  as  Then,  for  the  last 

FY,  we  will  have  the  correct  cumulative  total 

The  above  procedure  applies  equally  to  the  remaining  two  com- 
ponents. The  program  then  generates  the  c'unu iatith^  NOA  component  by 
component.  The  user  then  adds  vertically  by  FY,  and  takes  differences 
horizontal ly' to  get  the  final  NOA  by  FY  that  the  procurement  program 
would  demand. 

The  model  well  illustrates  the  simplicity  resulting  when  cumula- 
tive average  costing  is  used  rather  than  unit  costing. 


15.4.  PROORiVM  NOTES 

1.  l.BL  A uses  indirect  addressing  and  simple  looping  to  produce 
the  cumulative  flying  hours  H.. 

2.  l.Bh  B computes  required  cumulative  aircraft  and  puts  these 
in  the  primary  storage  originally  occupied  by  the  squadron  and  train- 
ing aircraft  of  the  original  schedule. 

3.  The  CTO  E of  line  074  will  lead  to  the  storage  of  I'l 

also  (see  above  for  reason). 

4.  LBl.  C calculates  the  cobf f Iclents  required  for  a segmented 
learning  curve. 

5.  LBl,  6 calculates  the  successive  cumulative  number  of  articles 
produced  and  NOA  required  as  tlie  fiscal  years  are  incremented. 

6.  "Packed"  storage  is  used  extensively.  Thus,  S.,  = 24  and 
T2  “ 6 are  stored  in  Register  2 as  24.06.  Tlien  S.,  is  retrieved  by 
'f  INT’  and  T^  by  'g  FR/\C,  EE.X  2,  X'.  This  storage  device  is  useful 
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who»  storago  space  is  at  a promiiim,  but  care  must  bo  taken  In  fixing 
the  number  of  places  alter  the  decLmal  point.  Thus  storing  24. b in- 
stead of  24.0b  would  cause  the  program  to  produce  bO  for  T.,  instead 
of  b.  You  also  pay  for  "packing"  in  the  coin  of  program  steps. 

?.  To  hold  the  ]>rogram  to  224  steps,  the  user  is  asked  to  do 
considerable  I'ut  simple  inputting,  such  as:  Oiven  a lead  time  t of 
2b  months,  key  in  2b  and  then  'KNTl.K,  12,  t,  STO  b'. 

tS.  For  lack  of  storage  space,  the  user  records  output  as  pro- 
duced and  does  some  final  additions  and  subtractions  to  get  the  NOA 
by  fiscal  year. 

Kxamjile.  I'he  required  buildup  schedule  is: 


Fnd  of  FY 

1 

•> 

:i 

4 

5 

b 

Squadron  aeft 

12 

24 

48 

84 

108 

108 

fraining  aeft 

1 

b 

12 

2\ 

27 

27 

The  steady-state  fleet  is  ti>  be  re.aclu'd  at  the  end  of  vear  5.  This 
will  be  kept  constant  through  year  10.  KH/Y  = 240  for  squadron  aeft 
and  720  for  training  aeft.  The  cumulative  attrition  coefficient  is 
O.OOOIS  and  the  exponent  is  1.05.  The  command  support  factor  is  5 
percent.  The  engine  learning  curve  has  three  segments,  and  the  param- 
eters for  this  component  are 

'*1  " '’ll  ° ’'ll  " '’l2  " ’'12  ° 

'’ll  ° *'1  “ '‘l  " • 

For  tiu'  airframe  (two  segments): 

■'2  “ '’21  ' ^-.1  = P2->  ° ^’-7’,  t.,  = 20.  u,,  = 1 . 

For  avionics  (one  segment): 

a,  = 10,  p^j  = 0.75,  t,  = 14,  U,  = 1 . 


3 


1 


4 


FH  YR  1 


FH/Y 


CUM  FH 
N ? 


CUM  FH 
M MORE  YRS 


REQ  A/C 
YR  1 


CUM  A/C 
REQ  N + M? 


STORED  IN 
PRl  I TO  N+1 


ATTRITION 
CMD  SUPP 


(Sec.  storage) 

(Overwrite  sec  storage  with  learning  curve  constants) 


1,  2,  3 
SEGMENTS  ? 


CHANGE  LEARNING 
CURVE  %'$  TO  b's 


FY 

STO  B 


CUM  COST 
1 SEGMENT? 


CUM  COST 
2 SEGMENTS  ? 


CUM  COST 
3 SEGMENTS 


/ 

Y 

r‘AI  r*l  11  ATP 

8 

Y 

r* A 1 II  ATP 

DSP  FY 

INCREMENT  FY 


(Note  multiple  use  of  registers  in  this  program.) 

Fig.  15.1  — Time  phased  procurement  flowcSort 


15.5  USER  INSTRUCTIONS 


15.  TIME  PHASED  PROCUREMENT  COSTING 


INSTP0CT1ONS  EXAMPLE  OF  14.4 

OATA/UMT^ 

KEVS 

STO  N.M  IN  A.  N YEARS  TO  FINAL 

1 i r 

BUILDUP,  M ADDITIONAL  YEARS  OF 

i II  1 

STEADY  STATE  FLEET  OPERATION 

i_..:)i:  1 

(N<  7,  NOTE  5.05  AND  NOT  5.5): 

5.05 

ISTOIIA  1 

STO  1 + X IN  B,  WHERE  X IS  THE 

1 if  ! 

COMMAND  SUPPORT  FACTOR 

1.05 

1 STO  II  B ) 

STO  C IN  C AND  d IND. 

.00015 

1 STO  i Ic  ! 

ATTRITION  FACTORS.  EQ  (2). 

1.05 

(STOl  1 D 1 

STO  A.  B IN  0.  A IS  FHA  PER 

240.720 

1 STO  1 1 0 1 

SO.  A/C,  B IS  FH/Y  PER  TRG  A/C. 

1 II  1 

(B  HAS  3 DIGITS.  INITIAL  ZERO  IF  NEEDED) 

1 II  1 

STO  S„  . Tn  IN  Rn  , n = 1,  ....  n. 

12.03 

1 STO  1 1 1 1 

AND  STO  Sn.  Tn  IN  Rn+,ALSO. 

24.06 

1 STO  1 1 2 1 

S„  IS  REQD  SO  A/C  AT  END  OF  YR  n . 

48.12 

ISTOll  3 1 

Tn  IS  REQD  TRG  A/C  AT  END  OF  YR  n- 

84.21 

1 STO  II  4 1 

108.27 

1 STO  1 1 5 1 

ISTOlI  6 1 

PRESS  A. 

1 1 1 A 1 

TO  SEE  CUM  FH 

1 f IIP— s| 

1 RCL 1 1 1 1 

IRCLII  2 1 

iRCLll  3 1 

1 RCL  1 ( 4'  1 

IRCLll  5 1 

1 RCL 1 1 6 1 

f P^S(IMPORTANT) 

1 f IIP— Si 

PRESS  B 

1 1 1 1 

TO  SEE  CUM.  A/C  REQ  (UNROUNDED) 

IRCLll  1.  1 

IRCLll  2 1 

lRCLir.3  1 

IRCLll  4 J 

1 RCLl 1 5 1 

1 RCL] 1 6 1 

f p^S  (IMPORTANT).  FOR  ENGINES: 

1 f 1 ip— si 

at  STO  0 

10 

isToi  1 0 i 

P,,  STO  1 

.9 

ISTOll  L 1 

X,,  STO  2 

60 

ISTOll  2 1 

P,2  STO  3 

.8 

1 STO  1 1 3 1 

1 II  1 

OUTPUT 

OATA/UMTS 


J^^05 

1.5  -04 

’-05 
240.72 


12.03 
24.06 
48.12 
84.21 
108.27 
108.27 
45360 
45360 
50^ 
15120 
35280 
70560 
115920 
342720 
342 TTO 

238.97 
16.91 


35.17 

71.93 

128.75 

172.W 

238.97 

238.97 

10.00 

M.OO 

o780 


15.5  USER  INSTRUCTIONS 


88§88g^S8 


STEP  KEY  ENTRY  KEY  CODE 


BBl  tLBLA 


15.6  TIME  PHASED  PROCUREMENT 

COMMENTS  STEP  KEY  ENTRY  KEY  CODE 


£ 


COMMENTS 


1 


OR  GTO  1 


FHA 

INDIRECT 

Si 


ASj  + BTj 


ISZI  16  26  46 


DSZI  16  ZS  46 


ISZI  16  26  46 


PARTIAL  SUM 
PR  I 


FINISHED  FOR  N? 
LOOP 


ISZI  16  26  46 


M.  (ASn+BTn) 
SEC 

Hn  (1) 

TOTAL  FH 


ISZI  16  26  46 


N + 1 

FINISHED? 

LOOP 


(1+X)  (Sn+Tn) 

SEC 

Hn 


an  (2) 

Qn  (3) 

PRI 

RTN  TO  063  / 067 


SEC  (USER  INST  9) 


1 SEGMENT? 


' S,,T, 

^ S2 

la 

3 

4 

5 

6 

“ X 

11 

S3  - 

Pl2 

Xi? 

S5  _ 

Pl3 

t,/12 

S7 

® Sn.  Tn 

9 

se 

-ni- 


ls.6 PROGRAM  USTING 


Nf 

Mr 

M' 

M.- 

irr 

is: 

js: 

Im 

ISS 

.'Sr 

155 


« 

5 TOC 
$roE 
5 TO  I 

p:l6 

TNT 
CHS 
STOP 
IfCLC 

me 

5 TOC 

P**S 

e 

STOP 

7TT“inFnr 


•0 

35  J4 
35  15 

35  46 

36  $6 
l>i  34 

35 

36  $6 
16  44 

35  13 
16-51 

M 

35  M 


CLEAR 

1/12 

FY 

PRI 


Q n 

(DELETE  161) 
Q n ♦ 1 

FRAC 


Trr 

T 

il9' 

'129 

16-51 

'll-' 

RCL4 

36  04 

:rc 

RCL9 

36  09 

l.'S 

VO’" 

16-33 

2 SEGMENTS  ? 

1 - • 

V 

-35 

I 

CTOC 

22  14 

ir: 

PSF0 

-63  00 

n: 

PPL  3 

36  03 

1 *'7 

PNC 

16  24 

:i5 

PCI  5 

36  05 

r4 

STOP 

35  14 

(8) 

1!.' 

- 

-45 

r?P7 

-63  02 

l.v 

;:o 

V» 

31 

t Tt* 

RCL2 

36  02 

X 1 1 

151 

PCM’ 

36  0.' 

♦ 

.MV 

-41 

♦ fc-' 

\ 

-35 

» **c 

.VtV’ 

16-35 

5 TOO 

35  00 

ir? 

CTCO 

22  00 

1st  SEGMENT 

* 2' 

CTOC 

22  14 

'l»0 

ISC 

PCI  4 

36  04 

"•IS!.! 

il  05 

- 

is; 

.\H' 

-41 

i 

isri 

16  26  46 

' I C" 

VtV' 

16-35 

i;' 

PCLI 

36  46 

I 

CTOS 

22  00  ' 

2nd  SEGMENT 

. 

00 

• ■ " 

« r * 

PC:  5 

36  05  ' 

l.'S 

.\=V’ 

16-33 

ALL  CONSTANTS? 

:35 

I 

01 

V 

tvU' 

;7i’ 

PTN 

24  ■ 

If  5 

•f 

-55  ■ 

• 7 ♦ 

RCL  i 

36  45’ 

Pii 

V' 

31  1 

4 « - 

LH 

32  ■ 

;sr 

PCS 

36  00  ' 

155 

02  ■ 

;s= 

\ 

-35  ■ 

3rd  SEGMENT 

I'l' 

IN 

32 

190 

■ 4 0 
. ri 

STOE 

35  15 

f Jf 

-24 

♦ r • 

CT09 

22  09 

;7;' 

co; 

35  45 

b|i 

is: 

• LflS 

21  09 

;5' 

:s:i 

16  26  46 

RCLB 

3i> 

:.r 

RCL  i 

36  45 

><li 

* ^ « 

F $ 

51  ■ 

DISPLAY  FY{152) 

12S 

EE.V 

-23" 

SEE  USER  INST.  11 

1 3r 

RCie 

3i>  12' 

Ul> 

MO 

5» 

09" 

" • C.' 

1 

01  ■ 

14; 

.\=V’ 

16-33 

. c*- 

♦ 

-55 

:4: 

PTN 

24 

ISf 

STOP 

35  12" 

: 45 

CTOS 

22  05  “ 

LOOP  FOR  SEGMENTS 

l.'S 

MS 

16-51  ' 

.Cl 


CO-' 


.ll. 

rcr 

TTH- 


fT06 


KCil 


S'> 

RCIO 

\ 

STOE 

CTO? 

'TUU 


22  96 
IT  87 
36  01 
01 
-55 
31 
36  00 
-35 
35  15 
22  09 
"TTW 


^^2  s 

H|V  B' 

B 

|9 

B 

E 

|9 

S 

I M 

^^9  S 

LABELS 


[p 


"TTWr 

Ki 

Ri 

:s:i 

STOl 

gTW 


FLAPS 


-31 
-31 
10  26  46 
35  45 



1st  SEGMENT 


2nd  SEGMENT 


TO  GET 


Qjj 


11. 


SET  STATUS 


FLAPS 


TRIP 


DISP 


ON  Off 

0 n (1 
<1111 
S 11  11 
3 1.1  II 


DEG  1 1 
GRAD  11 
RAD  1 1 


FIX  11 
SCI  1 1 
ENG  1 I 
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L^_-  _ c:0S'17»KNKKIT  STRKA>tS 


1 6 .J  ^ 

None. 

1().2.  DKSg'SSlON 

The  model  of  tills  section  is  also  due  to  H.  C. . Massey  of  Tlie 
Rand  Corporat  Isin . It  de.ils  with  decisions  to  spend  money  now  as 
opposed  to  later  during  tlie  life  cycle  of  a weapon  system.  Tliat  is, 
should  engineering  development  monies  he  spent  now  with  the  expecta- 
tion that  future  operating  and  support  costs  will  be  lower?  The 
planner  must  decide,  for  example,  wliether  to  install  engine  diagnostic 
equipment  now,  assuming  that  future  maintenance  will  otherwise  be  less 
efficient  and  therefore  more  costly.  Tlie  model  quantifies  such  de- 
cision problems,  using  as  a yardstick  the  present  value  of  a dis- 
counted stream  of  expenditures  and  benefits,  both  of  which  are 
e.xpressed  in  constant  dollars  (no  allowance  for  inflation). 

A simple  example  will  illustrate.  Suppose  we  assume  that  $10M 
spent  now  will  lead  to  operating  .ind  support  (O&S)  costs  of  $2riM  8 
years  from  now;  if  no  money  is  spent  now,  tliese  future  costs  will  be 
$40M.  Let  the  discount  rate  be  10  percent,  as  lairrently  mandated  by 
the  Oepartment  of  Oefense.  If  the  SIOM  were  invested  at  10  percent 
compounded  interest  for  8 years,  it  would  yield  $21. AM.  Oonsequent  ly , 
we  would  save  $l.4M  by  not  improving  the  system  now.  Hut  if  the  rate 
were  5 percent,  we  would  l-Ot'C  $5.2M  by  lu't  improving  the  system  now. 

It  follows  that  tlie  rate  mandated  or  assumed  has  a controlling  impact 
on  the  decision. 

With  respect  to  estimating  the  future  ciists  ($20M  and  $40M  above), 
one  other  crucial  point  must  be  maile.  .Suppose  the  system  in  question 
is  a fleet  of  aircraft  of  a given  tvpe.  Then  we  must  keep  the  oper- 
ational capability  constant  in  the  two  cases.  Th.it  is,  tlie  costs 

must  be  assessed  in  both  cases  to  produce  the  s.inie  in-commission  rate 
or  other  measure  of  operational  capability  for  the  fleet. 


r 


Set  a time  horizon  of  N future  fiscal  years,  tlu*  expected  life 
of  the  weapon  system.  Let  Cj  be  the  costs  of  the  proposed  near-term 
improvements  for  years  1,  ....  N.  t.et  be  the  assessed  Cn^'ri'trunitt'd 
O&S  ffaiHni/i'’  given  tlie  stream  C^.  Tlien  tlie  discounted  present  value 
of  the  benefit/cost  stream  at  a rate  r is 


l’(r)  = S (B,  - t^)(l  + r)‘ 
i=l 


since  by  convention  present  values  are  staled  for  fiscal  year  1. 
The  "breakeven"  year  j is  defined  as  tlie  first  i for  which 


(B  - C.)(l  + r)' 
i=l  ^ ‘ 


if  r(r)  > 0. 

The  "internal  rate  of  return"  is  that  r*  for  which 

P(r*)  = 0 . (I 

If  the  actual  rate  r is  greater  than  r*,  then  l’(r)  ^ 0,  because 
(I  + - (1  + In  this  case,  tlie  near-term  iniprovei'ient 

investment  would  not  be  madi>. 

The  value  r*  is  found  by  Newton's  method.  We  have,  using 

l)l  = B.  - C^, 


P'  (r)  * (I  - 1)  • 1)  • (1  + !■)■ 

i = l ‘ 


M 


P(r)  - i • I)  • (1  + r)' 
i = l 


f 


which  Is  wrltttMi  In  this  form  to  simplify  the' programming.  Finally 


1 + r = 1 + r , 
n n-1 


j|  * Vi>  • ''Vj.' 

N 

- E I • 0,  • (1  + r) 

i«l  ^ 


(5) 


•i+i 


16.4.  PRCXIRAM  notes 

(1)  Tills  program  uses  Indirect  addressing  in  a natural  and 
straightforward  fashion.  The  successive  O.'s  are  stored  in  registers 
I,  ....  i9  (19  years  is  the  maximum  for  this  program),  hooping  by 

f IS/,  is  then  simple. 

(2)  To  get  the  breakeven  year,  we  want  to  test  for  a change  of 
sign  in 

.1 

P(J.  r)  = 5^  I)  (1  + 
i=l 

as  this  partial  sum  crosses  the  time  a.xis.  A simple  procedure  is  to 
test  the  ratios  l’(,i  + 1,  r)/P(j,  r)  to  see  wlien,  if  at  all,  the  ratio 
Is  negativt^. 

(3)  In  using  Newton's  metliod,  a simple  first-guess  at  tjie  root 
is  r = 0.  But  I + r or  1 is  then  stored  in  B.  If  any  Iteration  pro- 
duces 1 + r < 1,  stop.  In  this  case  there  is  no  internal  rate  of 


return. 


INSTRUCTKMIS 


STO  N \t±A  

STO  Di  = Bi  -c,  IN_r ^ 

Dj. : 

. 03 

D4  _ 

05  _ 

IF  N > 9,  f P--S 

STO  D,o  IN  b _ 

STO  D,,  IN  1,  UP  TO  b 

f P— S (IMPORf ANT) 
STO  I + r IN  B (r  =J0%) 

PRESS  A. _ ^ BRjAKEVtN  YEAR 
PRESS  R/S.'  SEE“P(r) 


INPUT 

DATA/UNITS 


TO  FIND  r*,  THE  INTERNAL  RATE  OF  RETURN, 
STORE  A FIRST  GUESS  1 + r IN  B 
PRESS  B.  SEE  f -X-.  _ _ 

PRESS  R/S  WHEN  DISPLAY  REPEATS. 

(SHOWS  CONVERGENCE.)  ~ 


(THE  INTERNAL  RATE  OF  RETURN  IS  13.9%, 
WHICH  IS  GREATER  THAN  THE  MANDATED  10%. 
MAKE  THE  IMPROVEMENT.) 


KEYS 

OUTPUT 

0ATA/UNIT8 

|STO|  1 A 1 

5.00 

ISTOll  1 1 

- 5.bo'  ‘ 

|StO|  1 2 1 

-3.00 

ISTOI  ( '3  1 

1.00 

isTOli  4 1 

4.00 

isTOli  5 1 

7.00 

1 f 1 IP— Si 

1 II  1 

1 II  1 

1 f IIP— SI 

1 STO  1 1 B 1 

1.10 

1 ' 1 1 A 1 

‘5.000 

1 1 1 R/S  1 

0.886 

1 II  1 

1 1!  1 

ISTOll  B 1 

1.100 

1 1 1 B 1 

1.136 

1 II  1 

1.139 

1 II  1 

”Ti39 

1 II  1 

1 II  1 

1 II  1 

1 II  1 

1 II  I 

1 II  1 

1;  II  /I 

r II  ! 

'1  w 

1 l^^l 

k.1'  - II  1 

1 II  1 

\ 

1 II  1 

t 

r'  * 

1 II  1 

1 il  1 

1 II  1 

1 II  1 
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17.  THE  NORMAL  FUNCTION  AND  ITS  INVERSE 


17,1,  REFERENCES 


None , 


17.2,  DISCUSSION 

This  program  is  frequently  used  in  conjunction  with  others,  such 
as  the  Q function  (Program  18) , Extensive  application  of  the  error 
function  is  made  in  Hewlett  Packard  HP-65  Programs  for  Evaluating 
Effectiveness  of  Field  Artillery  Weapons,  prepared  for  the  Joint 
Technical  Coordinating  Group  for  Munitions  Effectiveness  (Surface- 
to-Surface)  by  Booz-Allen  Applied  Research,  Shalimar,  Florida, 
September  1975. 

17.3,  EQUATIONS 

The  normal  distribution,  or  function,  is 


F(x)  = 


exp  (-t  /2)  dt  , 


where  x is  in  units  of  a.  An  early  approximation,  due  to  J.  D.  Williams, 


F„(x)  - 1 t ^ 0 , 


which  has  a maximum  error  of  about  0,002,  R,  N,  Snow  replaces  the 
curly  brackets  in  (2)  by 


and  obtains  an  accuracy  better  than  0,0001, 


r 


ipBJlIJil  . «|J  Jfl 


r:=j:7-; 


I 

I 

i 

i 


i 


'1 

1 
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A first  approximation  to  the  inverse  is  obtained  by  solving  (2) 
for  X, 


^0 


TT  Ct!  {1  - (2Fq  - 1)^}'' 


1/2 


(3) 


For  greater  accuracy,  x^  is  used  in  solving  the  equation  F(x)  - y = 0 
by  Newton's  method,  which  is  appropriate  since  the  derivative  is 
simple , 


F^ (x)  = — ^ exp  (-x^/2)  . 


The  error  function  is 


erf (x) 


and 


(4) 


erf(x)  = 2 • F(/2  • x)  - 1 . (5) 

Given  a value  E to  find  the  corresponding  x,  use  this  program  with 
F = 1/2  [1  + E]  and  divide  the  resulting  value  by  /2. 

Example  1.  Find  erf (0.5). 

Key  in  1//2,  PRESS  A,  then  2,  x,  1 - . 

Answer  is  0.5206. 


Example  2.  Find  x for  erf(x)  = 0.5206. 

Key  in  F = (1  + 0.5206)/2,  PRESS  B,  PRESS  C,  /2,  f . 
Answer  is  0.50. 


17.4.  PROGRAM  NOTES 


None. 


17.5  USER  INSTRUCTIONS 


17.  THE  NORMAL  FUNCTION  AND  ITS  INVERSE 


INSTRUCTIONS 


INPUT 

DATA/UNITS 


OUTPUT 

OATA/UNIT8 


KEY  X 


PRESS  A 


OUTPUT  IS  F(x) 


KEY  F 


PRESS  B 


OUTPUT  IS  X 


NOW  PRESS  ^FO^i^RE^INED  VALUE  OF  x, 
AND  CONTINUE  IF  DESIRED 


(PRESS  A AGAIN  TO  RETRIEVE  F) 
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18.  THE  Q FUNCTION  (OFFSET  COVERAGE  FUNCTION) 


18.1.  REFERENCES 

J.  I.  Marcum  and  P.  Swerling,  "Studies  of  Target  Detection 
by  Pulsed  Radar,"  IRE  Transactions  on  Information  Theory, 

Voi.  IT-6,  No.  2,  April  1960.  (Reprints  of  Rand  RM-754, 
December  1947,  and  RM-1217,  March  1954.) 

b.  J.  I.  Marcum,  Tables  of  Q Functions,  The  Rand  Corporation, 
RM-339  (ASTIA  No.  AD  116551),  January  1950. 

c.  D.  P.  Meyer  and  H.  A.  Mayer,  Radar  Detection,  Academic  Press, 
New  York  and  London,  1973. 

d.  L.  A.  Wainstein  and  V.  D.  Zubakov,  Extraction  of  Signals 
from  Noise,  Prentice-Hall,  Englewood  Cliffs,  N.J.,  1962. 

e.  L.  E.  Brennan  and  I.  S.  Reed,  "An  Iterative  Method  of  Com- 
puting the  Q Function,"  IRE  Transactions  on  Information 
Theory,  Vol.  IT-11,  No.  2,  April  1965. 


18.2.  DISCUSSION 
The  Q function 


g(r,R)  = 


(1) 


is  basic  in  radar  detection  theory.  It  is  expressible  in  Lommel 
functions  of  the  first  kind  (Ref.  a),  but  is  not  integrable  in  closed 
form. 

The  Q function's  mure  common  application  is  perhaps  in  offset 
bombing  calculations.  For  a circular  normal  distribution  (0,a),  a 
weapon  radius  R (in  units  of  a),  and  a point  target  at  a distance  r 
(in  units  of  a)  from  the  origin  (the  aiming  point),  the  probability 
of  coverage  is  simply  P(R,r)  = 1 - Q(r,R).  If  CEP  in  feet  is  used 
and  r',  R'  are  in  feet. 


^2  B.n  2 


CEP 


R = R 


, /2  Rn  2 


CEP 


TIk'  damano  probablltlv  prom'am  (.Soo.  9)  can  be  iisi'd  In  )’cl  wi-apon 
radius,  and  this  pri>xram  Is  ihon  ciiiplovi’d  ti>  I'ind  cnllati'ral  ilamap.c 
ti>  ether  point  tarp.ets. 

1 8 . biyiAl  IONS 

'rite  bessel  t unction  given  bv 


l„U)  - E t’5)‘"/(n!)‘  . l.>) 

n-0  ‘ 

I’ut  (2)  ill  (I)  and  i ntercliaiig.e  lln'  order  ol  siiiiunal  ion  .iiui  Integration 
to  get 


'IK.r)  - 1 - Q(r,K) 


ir’^/; 


n=*l) 


K (K“/.’) 
n 


l )) 


wile  re 


k Iv)  ■ x''t'  ^ /n! 
n • 


.'ItL 

n 


and  I’u^]  l''n  IncompU'te  gamma  lunctlon. 

The  recursion  relations  t'or  k and  K are 

n n 


kjjly)  - e 


k I V ) ••  ' k , I V ) , n ' 0 


"X 

I' 


K U)  ■ I 


AJ  tov  N itcratU'ns, 


N-l 

l’(.R,r)  - 52  R + « 

n u 

n-  1 


whoro  the  roma liuior 


K(N> 


n-N 


k K 
n n 


CIO 


Rt'l  oronvo  o shows,  ti'r  N '•  rK/»0.  Ih.U 


R(,N  + \)  ^ k^(r-’/2)  K^(RV2)  . {7) 

hot.  ”•  rR/v^2  .iiul  lot.  A ho  tho  doslrod  .u'oiitMoy . Thou  llio 
itotMtlon  oan  ho  totminalovi  at  n » N,  It 

N " N,,  . k„Cr"/21  K^,tR"/2)  S A . (S) 

tl  N N 

It'  r atui  R aro  small,  tho  oi'twoi>\oiu'i'  Is  rapUl.  lUil  as  r a>id  R 
inoroaso,  tho  mimhor  ol  torms  roqiiirod — anil  thus  oi'mput  at  1 ou  t imo — 
hooonu-  oxoosstvo. 

Howovor,  for  i,  R I,  Rot.  li  provUlos  au  oxoollottt  approximation 
in  torms  of  tlto  ovamil.it  ivo  itansstan  fnnot  ton 


rCK,r) 


iln  + K.CX,  O 


t'O 
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The  approximation  is  improved  by  adding  two  more  terms  to  (10), 
taking 

= (R  - r)(l  - 1/Ar^)'-  l/2rf-  l/48r^  - (R  - r)^/H)r^  . (11) 

I 

Then  the  error  E(Z^)  S 0.0005  for  r i 4.25.  It  can  also  be  shown  that 
for 


R a 2.45  + - 
r 


_1 

- 1.7  ’ 


E(Z^)  £ 0.005  . 


(12) 


The  regions  of  interest  for  the  P(R,r)  calculation  are  tlion 
shown  in  Kig.  18.1.  For  region  A,  use  the  Iteration  method,  Eqs.  (3) 
and  (8).  The  maximum  number  of  iterations  required  to  obtain  an 
error  A ^ 0.0005  is  10.  For  region  B,  use  the  approximation  (11)  in 
(9),  with  an  execution  time  of  7 sec.  In  general,  tlie  regions  above 
R - r = -2.8  .and  below  R - r = 2.8  are  without  interest.  However, 
even  in  these  regions,  the  approximation  method  tmiy  be  used. 

18.4.  PROGRAM  NOTES 

1.  For  an  accuracy  A S 0.0005,  enter  R,  r,  and  A and  press  B. 

The  program  will  choose  the  better  method  to  compute  P(R,r).  "Better" 
is  defined  as  entailing  the  least  computer  time  to  obtain  at  least 
accuracy  A. 

2.  For  an  accuracy  A better  tlian  0.005,  again  enter  R, 

A but  press  A. 


r , and 


18.6  fROGRAM  LISTING 

COMMENTS  STEP  KEY  ENTRY  KEY  CODE  COMMEI 


STEP  KEY  ENTRY  KEY  CODE 


1 SET  STATUS  1 

FLAGS 

TWO 

OtSP 

ON  OFF 

ODD 

DEG  D 

FIX  D 

) D D 

GRAD  D 

SCI  D 

1 D D 

RAD  D 

ENG  D 

3 D U 

0 
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19.  LINEAR  PROGRAMMING  AND  3x3  MATRIX  GAMES 


19.1.  REFERENCES 

a.  G.  B.  Dantzlg,  Linear  Pvograjnmitvj  arui  Extensions , The  Rand 
Gorporation,  R-366-PR,  August  1963  (published  by  Prince- 
ton University  Press). 

b.  A.  M.  Glickman,  An  Introduetioyi  to  Linear  Prograrmiruj  ami 
the  Theory  of  Gaines,  John  Wiley  and  Sons,  New  York,  1963. 

c.  R.  W.  Metzger,  Elementary  Mathematiaal  Programming,  John 
Wiley  and  Sons,  New  York,  1958. 

d.  J.  D.  Williams,  The  Compleat  Strategyst,  McGraw-Hill,  New 
York,  revised  edition,  1966. 

e.  M.  Dresher,  Games  of  Strategy:  Theory  and  Applications, 

The  Rand  Corporation,  R-360,  May  1961  (published  by  Prentice- 
Hall). 


19.2.  DISCUSSION 

. The  3-activity  linear  programming  problem  may  be  formulated: 


find  ^0,  X.,  & 0,  S 0 satisfying 


a^l  x^  -b  a^2  >^2  + ‘^3  ^3  ^ ^ 


®21  ^1  ^22  ^2  ‘^23  ^3  ^ ’^2 


^31  ^32  ^2  ^33  ^3  ^ *^3 


to  maximize 


M = Cj^Xj^  -b  <^2^2  *^3^3 


The  recipe  for  solving  this  problem  by  the  pivot  method  is  very 
simple,  although  the  result  of  deep  and  extensive  analysis.  Set  up 
Tableau  1: 


TABLEAU  1 


1.  The  pivot  aolurm  corresponds  to  the  most  negative  of 

i "‘^2’  this  is  -c^. 

2.  The  pivot  row  is  found  by  finding 

' bi/ai2  , ^2/^22  • ^3^®32 

s 

J 

for  only  those  b's  that  are  positive,  and  then  selecting  the  minimum. 
Say  this  is  b2/a22*  Then  the  pivot  is  222' 

3.  Replace  the  pivot  by  its  reciprocal,  and  divide  all  other 
entries  in  the  pivot's  column  by  the  negative  of  the  pivot. 

4.  For  an  entry  other  than  those  in  the  pivot's  row  and  column, 
add  to  that  entry  the  product  of  the  entry  in  the  same  column  to  the 
left  or  right  of  the  pivot  in  its  row  and  the  entry  in  the  (new) 
pivot's  column  to  the  left  or  right  of  the  entry  in  its  row. 

5.  Now  modify  the  pivot's  row  by  dividing  all  entries  other 
than  the  pivot  by  the  pivot. 

, 6.  Interchange  X2  and  u^. 

These  operations  produce  the  new  Tableau  ?: 


X 

1 

*3 

X 

3 

11 

(-*12/*32) 

"*12'^*32 

*13  *33 

(-312/332) 

^ ^3  ^"^12^^32 

21  *31 

(”^22/ ^32^ 

■*22^*32 

*23  *33 

(—322/ ^32^ 

^2  ^3  ^■'^22^'*32 

21^*32 

l/*32 

*33^*32 

'’3'^*32 

■=1  *31 

02/3^2 

-C3  + a^j 

0 + bj 

TABLEAU  2 


In  the  numerical  example  below,  the  initial  pivot  is  {not 


a^2  above) . 


Example  (Ref.  b,  Sec.  5) 


^2  =^3 


1 

1 

1 

100 

3 

2 

© 

210 

3 

2 

0 

150 

-5 

-6 

0 

=^2 

*2 

1/4 

1/2 

-1/4 

47.5 

3/4 

1/2 

1/4 

52.5 

3 

© 

0 

150.0 

-1/2 

-1 

3/2 

315.0 

(Tableau  1) 


(Tableau  2) 


7.  Repeat  the  above  process  until  all  entries  in  the  bottom 
row  are  nonnegative-  In  the  rightmost  column  read  the  optimal  x's 
and  maximum  M.  Any  u left  in  that  column  indicates  that  the  corres- 
ponding X is  0. 

Continuing  the  example. 


.\i 


*1 

“3 

^2 

-1/2 

-1/4 

-1/4 

10 

0 

-1/4 

1/4 

15 

3/2 

1/2 

0 

75 

1 

1/2 

3/2 

390 

= 0,  X2 

= 75,  X3 

= 15. 

M = 390 

(Tableau  3) 

The  previous  example  in  equation  form  is: 

find  Xj^  i 0,  X2  ^ 0,  x^  S 0 subject  to 


Xi  + X2  + x^  S 100 


3x^  + 2X2  + Ax^  s 210 


3x^  + 2x, 


s:  150 


to  ma.-rzrnzc  M = 3x^  + 4X2  + bx^ 


Suppose  the  problem  were: 

find  x^  i 0,  X2  i 0,  x^  ^ 0 subject  to 


Xf  + X2  + Xj  a 100 


3Xj^  + 2x2  + 4x^  i 210 


3Xj^  + 2X2 


S 150 


to  minimize  m » 5x^  + 4X2  + 6x^. 


( 
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We  set  up  tl»e  dual  problem  (tlie  pattern  should  be  clear): 

find  ^ 0,  y,  S 0,  i 0 subject  to 
y^^  + 3y^  + 3y^  S 5 
>’1  + -y-y  + 2y^  S 4 
V,  + Ay,,  S 6 

to  rruiximise  M = lOOy,^  + 210y2  + 150v^,  and  solve  this  problem  by  the 
previous  method.  The  values  for  x variables  are  now  read  off  on  the 
bottom  row  under  the  interchanged  u's  with  the  same  subscript. 

The  relation  of  the  duals  is  clarified  by  using  the  program  to 
solve  3x3  matrix  games.  The  procedure  is  readily  understood  by 
following  through  a specific  example. 

BLUE  is  the  maximizing  player,  RED  his  minimizing  opponent.  It 
is  desirable  to  add  a constant  to  all  entries  in  the  payoff  matrix  to 
make  all  entries  positive  if  nece.ssary. 

This  increases  the  value  of  the  original  game  by  that  constant 
but  does  not  change  the  proportions  in  which  the  strategies  are  played. 
The  new  example  is 

RED 
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Look  at  matters  from  RED's  point  of  view.  Against  each  of 
blue's  f’urt:  strategies,  RED  must  expect  to  pay  BLUE 


2r^  + 6r^,  Sr^^  + 3r2  + 6r^,  5r^  + 4r2  + Sr^ 


Let  p (unknown)  be  the  greatest  of  these  three.  Then  putting 

Yl  = r^/u. 


2yi  + 6y2  S 1,  5yj^  + 3y2  + 6y^  Si,  5y^  + hy  ^ + 3y^  S 1 , 


and  RED  wants  to  minimize  u by  maximiziyuj  M = 1/p, 

M = + y2  + y3  • 

This  is  the  standard  Linear  programming  problem  with  the  right-upper 
border  all  +l's  and  the  left  lower  border  all  -I's.  The  first 
tableau  is: 


^1 

^2 

^3 

2 

© 

0 

1 

5 

3 

6 

1 

^2 

5 

4 

3 

1 

'^3 

-1 

-1 

-1 

0 

M 

Any  column  could  be  the  first  pivot  column.  We  choose  the  circled 
number  as  the  pivot  because  1/6  < 1/5  (a  variation  on  step  (2)  above), 
The  second  tableau  (by  the  program)  is: 


U Is  >4lvoii  ii\  tlK'  hluli-lovol  lanmiago  appropriate  to  ovir  minds,  I'lit 
only  witU  oaro  and  time  do  wo  avoid  errors. 

The  llF-(i/  understands  iinly  a relatively  low-level  lannuage  even 
tliou>;li  it  is  quite  advanced  over  early  coding  in  machine  language. 

So  wh.it  we  see  without  et't'v'rt  (the  Mciit  a!  t ) requires  many  steps  of 
programming. 

1.  I'he  programming  prohU'iii  is  compounded  in  difficulty  because 
the  iMiiniit  work  with  the  diuiliU'  indices  of  matrix  elenu'nts.  Wi> 

must  " 1 i near  I .’e"  t lu'  m.itrix  ,ind  ilevise  an  "adilress  arithmetic"  tor  this 
problem,  .tiul  then  exploit  t lu'  indinst  .uldres-.ing  e.ip.ib  i 1 i t v I’l  the 
1 I eg  I St  el  . 

tonstdei  item  (,1  ot  the  le.  ip,'.  I'veilav  the  lineal  plogiamming 
I I » w I f li  t b,  mat  I ( \ ol  i eg  l ■. I ei  iddi  es',,“.  ot  ent  I i es  ; 


iiji  iTliini 


(SUBROUTINE) 


PIVOT  ADDRESS 
= 0,  I,  2 ? 


RECIPE 
STEPS  3,  5 


PROCESS  ENTRY 
ADDRESSES  IN  ORDER 


SAME  ROW 
AS  PIVOT  ? 


RtCifl  STEP  4 
All  INTRIIS 
I VRMINED  7 


4 

m 

4 

i r: 

s 

SAME  COLUMN 

f 

INCREMENT 

AS  PIVOT  ? 

ADDRESS 
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19.5  USER  INSTRUCTIONS 


19.  LINEAR  PROGRAAAMING 
AND  3X3  MATRIX  GAMES 


mnrr 

OATA/uurrs 


STEP 


INSTRUCTIONS 


I f CL  REG,  f P— 5,  f CL  REG 


2 LOAD  BOTH  SIDES  OF  CARD 


5 03),  032  STO  8,  9 


9 » P — S( IMPORTANT) 


(PRESS  E TO  REVIEW  ENTRIES) 


0 CM  TERMINI  PIVOT 


11  IN  PfV(iT  AOOtfSS 

C it  » S • • ••) 


IS  BEST  QUALITY  FRACTICABB8 
f UKNiijHiiD  rO  DDQ  ^ 


]9.6  LINEAR  PROGRAMMING 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


STEP  KEY  ENTRY  KEY  CODE 


COMMENTS 


901  *LBL9 
092  STOI 

003  STOS 

004  RCLi 

005  ]/X 
B0c'  CHS 

007  STO0 

008  RCLB 


$L8L4  2i  $4 


PIVOT  ADD(RESS) 


PIVOT  ADD  -4 


1 / PIVOT  VALUE 


PIVOT  ADD  +4 


PIVOT  ADD=  0,1,2 


PIVOT  ADD  + 8 


PIVOT  ADD 


PIVOT  ADD -8 


FIRST  ENTRY  ADD 


PIVOT  ADD 


PIVOT  ADD  -4 


ENTRY  ADD 


PIVOT  ADD + 4 


z (PRGM  NOTE  2) 


(RECIPE  STEPS  3,  5) 


COL.  ENTRY 
DIVIDED  BY-PIVOT 


TEST  FOR  ENTRY 

IN  SAME  COL  AS  PIVOT 


PIVOT  ADD  » 4 

(E1^0J_CQU 


TEST  FOR  ENTRY 

IN  SAMI  HOW  AS  PIVOT 


PIVOT  ADO  » 8 
(PIVOT  COL) 


PtVO!  AOO  17 

nvt  'T  It 


INTIRC  HANOI 


|1HS  FAdE  IS  B£ST  QUAMlY  FRACaClCABlJ 
OOPI  f URIUSHED  10  JIDC 


19.6  PROGRAM  LISTING 


COMMENTS 


STEP  KEY  ENTRY  KEY  CODE 


STEP  KEY  ENTRY  Kf  V CODE 


COMMENTS 


rqv/jnir: 


$LBL2  21  12 

RCLB  36  12 


(RECIPE  STEP  4) 


14  24  44 

14  24  44  ROW  ADD  TO  RIGHT 

. -1.3  H 

14  24  44 

23  N LAST  ROW  ENTRY 


ALL  ENTRIES? 


increment 


ALL  ENTRIES? 


RECIP.  OF  PIVOT 


PIVOT  ADD=0,4,8 


PIVOT  ADD  = 1,5,  9 


1 ST  ROW  ENTRY 


2ND  ROW  ENTRY 


ROW  ENTRY  TO  RIGHT 


12S 

RCLC 

34  13 

13? 

1 

«1 

J30 

-55 

13; 

STOC 

35  13 

132 

CTOD 

22  14 

DIG 

GNM' 

at' 

n 

MT  ITATVt 

ElAM  TMO 

oma  1 

/ 

i 
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20.  FOURTH-ORDER  DIFFERENTIAL  EQUATIONS 


20.1.  REFERENCE 

a.  R.  W.  Hanuning,  Numeriaal  Methods  for  Saientists  and  Engineers, 
McGraw-Hill,  New  York,  1962. 


20.2.  DISCUSSION 

Systems  of  four  first-order  (frequently  nonlinear)  differential 
equa^^ns,  sometimes  in  the  guise  of  two  second-order  equations,  occur 
more  commonly  than  one  would  think:  The  basic  beam  deflection  equa- 
tion is  of  the  fourth  order;  chemical  kinetic  systems  with  four  (or 
more)  equations  are  common;  reentry  trajectories  are  specified  by  two 
second-order  equations  (see  Program  4);  Lanchester  equations  with  two 
force  components  *For  ..^J^^jj'side  and  with  variable  coefficients  occur; 
problems  in  cf^Aal  control  theory  and  in  differential  game  theory 
lead  to  such  systems. 

Since  the  HP-67  has  a limited  number  of  storage  registers  (26) 
and  of  program  steps  (224),  and  since  program  space  is  needed  to  de- 
fine the  functions  oti'che  system,  we  seek  an  alternative  to  the  rather 
complex  Runge-Kutta  "standard"  formulas,  an  alternative  that  is  miserly 
of  program  space  and  yet  has  good  accuracy  for  relatively  large  time 
intervals — that  is,  good  relative  stability,  defined  as  the  rate  of 
growth  of  the  error  relative  to  the  growth  of  the  solution. 

Moreover,  programming  must  fully  exploit  the  indirect  control 
afforded  by  the  powerful  I-register.  That  is,  the  number  in  the 
I-register  can  be  the  address  of  a storage  register  or  the  name  of 
a label  (subroutine).  Then  the  instruction  STO  (i)  or  RCL  (i)  or 
GTO  (1)  moves  X-register  data  to  the  right  register  or  recalls  data 
from  the  desired  register  or  sends  the  program  to  the  right  place. 
Hence,  in  conjunction  with  incrementing  and  decrementing  the  I- 
reglster,  serial  treatment  of  all  four  equations  can  he  accomplished 
with  the  same  economical  set  of  processing  instructions. 

The  final  programming  problem  is  to  move  data  around,  like  freight 
cars  in  a marshalling  yard,  so  that  storage  spaces  are  freed  just  in 
time  to  make  space  for  a new  claimant. 
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20.3.  EQUATIONS 

Section  14.3  of  Ref.  a describes  a simple  predictor-corrector 
approach  for  first-order  equations  that  seems  to  have  promise  and  i 
readily  extended  to  a system  of  equations. 

Consider  the  equation 


dx/dt  = X(x,t),  x(0)  = Xq  . 


Let  the  time  interval  be  h.  Suppose  that  at  time  (n  - l)h  w€ 
Then  a good  pvediated  value  for  appears  to  be 


P ..  = X , + 2h  x' 
n+1  n-1  n 


= Vl  . 


since  x^  is  the  slope  at  the  midpoint  of  the  double  interval.  Accord- 
ing to  Hamming,  the  error  term  is  h\*'(0)/3. 

The  value  p^^^  is  now  aorrcoted  by  taking 

"n-Hl  = . 


where  Pn+i^^s  used  to  determine  p^^^  = X(p^^^.t).  The  error  term  for 

^n-l-1  ^ (0)/12.  If  is  approximately  constant  in  the  interval, 

then 


'’..+1  - '„+l  - 5'’ 


■ ‘“n+l  - • 5 


Vi  • '‘'m-i  * 


'‘-r  ^ vi  • 


-188- 


Pn+1  ’ ^^Pu+1* « + 1 • ; 


This  is  done  by 


To  get  stiirted,  we  need  in  addition  to 
expanding  in  a Taylor's  series 


Xj^  = Xq  + hx^  + h^x'/2  + . . . 


3 4 

Hamming  reconunends  carrying  the  series  to  tlie  h or  h terms.  In  our 

■> 

applications  we  have  stopped  frequently  at  h“  because  of  the  labor  in 
computing  Xq  and  x^l  But  note  that  if  h/2  is  used,  the  error  is 
multiplied  by  1/8. 

The  preceding  analysts  is  readily  generalized  to  a system  of  four 
equat ions : 


x'  = X(x,y,u,v,t) 
u'  = U(x,y,u.v,t) 


y'  = Y(x,y,u,v,t) 
v'  = V(x,y,u,v,t) 


We  have 


p = X , + 2hx' 
n+l  n-1  n 


r . , = u , + 2hu' 
n+l  n-l  n 


q , , = V , + 2hv' 
n+l  n-l  n 


s , = V , + 2hv ' 
n+l  n-l  n 


P.;+1  ° ^^Pn+r  Vi*  Vi-  Vi‘  ''  ^ 


Vl  ’ ^ Pn+l  ^ P.;+l^'/^ 


20.4,  PRtX:RAM  NOTES 

1.  For  each  equation  we  need  temporary  storage  for  x^^  x^, 

^n’  Pn+1*  1®  the  register  contents  on  the  first  page  of 


-189- 


the  program  listing  shows  how  the  storage  is  laid  out.  Secondary 
(protected)  storage  must  of  course  be  used.  Fortunately,  in  indirect 
control  the  registers  SO  to  S9  become  10  to  19,  so  that  f P S 
switching  of  primary  and  secondary  registers  need  not  be  programmed. 

2.  Label  numbers  are  assigned  to  the  subroutines  that  compute 
X,  Y,  U,  V to  agree  with  the  register  numbers  where  x^,  y^,  u^, 
are  stored.  That  is,  LBL  3 ~ 3,  LBL  7 ~ 7,  LBL  B ~ 11,  LBL  A ~ 15 

in  the  indirect  control  mode.  Hence  if  we  are  working  with  the  equa- 
tion for  V,  say  with  15  in  the  I-register,  then  f GSB(i)  sends  us  to 
the  correct  routine  for  the  fourth  equation  of  the  set. 

3.  In  evaluating  the  functions  X,  Y,  U,  V,  the  locations  of  x, 

y,  u,  V are  assumed  by  the  programming  to  be  3,  7,  SI,  S5.  Hence  in 

evaluating  p'  ,,  etc.,  we  must  be  careful  to  move  x from  3 to  2,  etc., 
n+1  n 

and  then  from  5 into  3,  etc. 

A.  Extensive  use  of  ISZ  and  DSZ  is  made  to  drive  indirect  con- 
trol to  the  right  addresses  at  the  right  time. 

5.  Time  is  incremented  by  RCL  (}>,  STO  + 1. 

. 6.  To  keep  track  of  the  number  of  equations  in  the  set,  store 
m,  the  number  of  equations,  in  both  D and  E initially.  On  each  cycle, 
as  an  equation  is  processed,  the  value  in  D is  decreased  by  1.  A 
test  tells  if  all  equations  have  been  processed,  that  the  iteration 
is  complete,  that  time  may  be  Incremented,  and  D reset  to  m.  That 
is,  the  program  may  be  used  for  any  number  m of  equations  up  to  A. 

Example.  The  equations  for  the  motion  of  a particular  hydrostatic, 
pendulum  are  d^x/dt^  + 3 dy/dt  + Ax  = 0,  d^y/dt^  - 3dx/dt  + Ay  = 0, 
with  Xq  = 1,  y^  = 0,  Xq  = 0,  y^  = A.  You  can  almost  smell — but  not 
quite — that  the  exact  solution  is 

X = cos  At  y = sin  At  , 

which  is  indeed  the  case. 

The  equivalent  system  of  four  equations  is 

x'  = u y'  = V u'  = -Ax  - 3v  v'  = -Ay  + 3u  . 


-190- 


f 

t 


( 

1 

-t 

1 


The  values  for  t = h are  readily  found  by  expanding,  using 


to  get 


*0 


— 

0 


*0  = 


1 

0 

-16 

0 


^0  = 


y'o  = 


^0  = 


0 

4 

0 

-64 


u*  = 256 


V*  = 


>^0  = 


4 

0 

-64 

0 


= 1 - 16h  /2 


= -16h  + 256h^/6 


y^^  = 4h  - 64h  /6 


= 4 - 64h  /2  , 


which  are  recognized  as  the  leading  terms  In  the  series  for  cos  4t, 
sin  4t,  - 4 sin  4t,  4 cos  4t,  respectively.  For  h = 0.025  we  have, 
to  4 places. 


= 0.9950 


y^^  = 0.0998 


u^  = -0.3993  v^  = 3.9800  , 


but  these  values  are  stored  In  registers  3,  7,  SI,  S5  as  they  are 
computed  to  maintain  full  accuracy.  (Remember  f P +-►  S.)  Also  store 
the  Initial  values  In  2,  6,  SO,  S4,  store  m = 4 In  D and  E,  and  store 
h = 0.025  In  0 and  1.  Now  load  the  program. 

Next  program  the  functions  by  GT0.121,  switch  to  W/PRGM,  and  then 


f 

f 

f 

g 


LBL 

3, 

f P 

^ s. 

RCL  1,  f 

P 

s. 

h RTN; 

LBL 

7, 

f P 

^ s. 

RCL  5,  f 

P 

S, 

h RTN; 

LBL 

B, 

RCL 

3,  4, 

X,  CHS,f 

P 

*<— ► 

S, 

RCL  5, 

3,  X,  -, 

f P 

^ s. 

h RTN; 

LBL 

a: 

RCL 

7r.4, 

X,  CHS,f 

P 

S, 

RCL  1, 

3,  X,  -I-, 

f P S,  h RTN  . 


r 
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Switch  to  RUN  and  press  A.  In  the  following  tabulation  we  will 
record  only  x and  y at  intervals  of  2h  = 0.05.  The  numbers  beneath 
are  the  exact  values,  where  we  first  shift  to  the  radian  mode  by 
h RAD  to  get  cos  4t  and  sin  4t. 

h = 0.025  h = 0.05  h = 0.025  h = 0.05 


_t 

X 

X 

y 

V 

0.05 

.9801 

(.9801) 

.1986 

(.1987) 

0.10 

.9210 

(.9211) 

.9482 

.3894 

(.3894) 

.3158 

0.15 

.8253 

(.8253) 

.8687 

.5646 

(.5646) 

.4938 

0.20 

.6967 

(.6967) 

.7525 

.7173 

(.7174) 

.6573 

0.25 

.5403 

(.5403) 

.6070 

.8414 

(.8415) 

.7935 

0.30 

.3623 

(.3624) 

.4371 

.9319 

(.9320) 

.8981 

0.35 

.1699 

(.1700) 

.2499 

.9853 

(.9854) 

.9670 

0.40 

-.0292 

(-.0292) 

.0528 

.9994 

(.9996) 

.9972 

For  this  example  and  with  a spacing  of  h = 0.025,  accuracy  and 
absolute  stability  are  excellent,  although  the  running  time  is  about 
24  sec  per  iteration  and  21  sec  per  display  of  the  four  variables. 

On  the  other  hand,  a spacing  of  h = 0.05  car icaturizes  the  solution. 
This  points  up  the  wisdom  of  doing  a second  run  with  half  the  original 
interval  to  determine  whether  major  changes  are  occurring. 


6^ 
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21.  _ I'Amvi;  J^KS.  ANP,  mach  numbkrs 


21.1.  AKK1;RKN£KS 

a.  United  States  Air  Force,  Flt-jht  Manual  A-?P  Aivai'aJ’t,  T.U. 
lA-7l)-l.S-32,  29  March  1971. 

h.  United  States  Air  Force,  AFSC,  rlannavt-  uuiilt',  I July 

Idh"!  (For  Official  Use  Only). 

c.  0.  H.  Kaplan,  1..  K.  Donnett,  and  F.  K.  Seldelmann,  Ahnanaa 
i\'n!p}4ti't’a , 1077,  United  States  Naval  Observatory,  tMr- 
cular  No.  155,  I October  197b. 

Military  data  for  analytic  or  operational  use  are  fri'cpientlv 
presented  as  .i  family  c>f  curves  z » f(x,y),  wlieri>  y is  t lie  p.irameter 
namlti);  the  f.imlly's  members.  For  example,  Fij;,  21.1  (taken  from 
Kef.  a)  is  a nomogram  to  determine,  for  the  A-71)  aircraft,  the  Macli 
number  to  maximize  range  for  constant  altitude  cruise,  given  average 
gross  weight  and  drag  Index.  To  use  the  nomogram,  start  witli  the 
average  gross  weight  on  the  upper  left  scale.  Move  lu'r Izontal ly  to 
tlie  pressure  altitude.  Drop  vertically  to  tlie  appropriate  drag  in- 
dex curve.  Move  l»or izontal  ly  to  tl>e  left  to  read  tlie  true  Mach  num- 
ber for  long-range  cruise.  Page  .after  page  of  such  nomograms  appear 
In  the  mission-planning  appendix  of  Ref.  a.  Similarly,  page  after 
page  of  similar  nomogr.ims  app»'ar  in  Ref.  b,  to  be  used  in  pl.inning 
space  missions. 

Neltlier  Ref.  a nor  Ref.  b gives  the  equations  of  the  curve  fam- 
ilies. When  these  equations  c.an  be  found,  and  if  they  are  relatively 
simple,  they  can  be  readily  programmed.  Frequently,  liowever,  as  in 
the  case  of  an  epliemerls  or  almanac  that  t.ibulates  the  coonlinates 
of  celestial  bodies  for  astronomical  and  navigational  use,  the  under- 
lying equations  are  extremely  complex.  To  quote  from  Ref.  c,  these 

A 

The  drag  Index  Is  not  a drag  coefficient.  Its  determination, 
as  explained  In  detail  In  Ref.  a.  Is  a tabulation  of  the  dr.ig  con- 
tributions of  external  stores  by  type  and  station.  The  clean  air- 
craft has  a drag  index  of  0. 
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tabulations  "should  idoallv  be  replaced  by  conoise  nuit  liemat  ica  1 expres- 
sions for  direct  calculations.  Such  expressions  must  take  the  Ivtrm 
of  mathematical  approximations,  however,  since  the  precise  data  . . . 


are  calculated  from  extensive  theories. 


Hence,  the  Observa- 


tory's Circular  No.  135  provides  the  coefficients  for  the  efficient 
expansion  of  tabulations  In  a series  of  Chebyshev  polynomials.  (These 
are  ►k'f  curve  families,  but  functions  of  a slnxle  variable.) 

In  other  cases,  a curve  family  may  arise  from  empirical  observa- 
tions such  as  flight  tests,  or  may  be  the  presentation  of  the  results 
of  extensive  computer  simulations. 

In  all  cases,  the  question  Is:  Can  approx Im.at  Ions  be  found  for 

• * 

the  family  that  arc  of  aooiiracv  to  the  data,  and  that  arc 

well  within  the  capabilities  of  the  Hl’-n7  to  implement? 

Such  approximations  do  not  provide  insight  in  respect  ti’  the 
physical  nature  of  the  underlying  phenomenon  or  model  if  polvnomlal 
.ipproxlmat  Ions  are  used.  I’olvnomlals  are  usually  employed  because 
of  their  simplicity  and  convenience.  If,  hi'wever,  physical  arguments 
or  the  methods  of  dlraenslon.il  .inalvsls  suggest  a functional  rel.it  lon- 
shlp  .imong  the  variables,  this  rel.it  lonship  should  certalnlv  be  used 
In  the  approximation.  More  modestly,  as  we  shall  see  shortly,  it  is 
alw.iys  deslr.ible  to  consider  taking  the  logarithms  of  the  dependent 
V. triable  observables  ami  then  subjecting  these  values  to  a pi'lvnoml.il 


riie  obvious  .ippri'ach  to  finding  a polvnomlal  approx  lm.it  ion  ti'r 

tlie  set  of  observ.ibles  (x  , v , z,,)  is,  in  principle,  bv  the  method 
**  till 

of  least  squ.ires.  For  example,  let 


flx,v'  “ a,,  + a.x  F a ,v  + a .x"  + a,xv  + .i,v 

i*  I •.  » *4  » 


A'./:*  me.ins  .i  c.i  1 ini  l.it  ed  pesult  within  I pein-mit  or  J per- 

cent of  the  result  of  re.idlng  direct  Iv  fri'm  the  curves.  And  one  must 
keep  in  mimi  that  la)  the  width  ot  plotted  curves  and  t h»>  w;iv  in 
which  they  .ire  dr.iwn  c.in  liilri»duce  d Iscrep.inc  ies  of  this  orvler  in 
the  basic  re.idlngs  ot  v.ilues  on  wliich  to  base  an  .ipproximat  Ion,  .iiid 
lb)  the  accuracy  of  the  presented  data  Is  usually  unknown. 

■kik 

In  more  mi>dern  guise,  bv  orthonormal  Iz  lug  innles.  S»‘e  .John 
Tixld  led.),  A'ui!  Uift  Chap.  lO,  Nctlraw-Hi  1 1 , 

New  York,  I'lhJ. 
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be  the  predictor.  Then  choose  the  coefficients  to  minimize 


Take  the  six  partial  derivatives  of  G with  respect  to  tl»e  coefficients 
and  equate  eacl>  to  zero.  Form  ti»e  Indicated  sums  and  Ret  six  linear 

A 

algebraic  equations  for  the  coefficients  a^^,  ....  aj^. 

This  Is  principle.  In  practice  we  have  two  problems: 

1.  Even  this  "best"  approxlmatlon’by  a second-order  polynomial 
may  be  of  unacceptable  accuracy. 

2.  No  method  In  two  dimensions  is  known  that  is  equl.’ilent  to 
the  use  of  Chebyshev  polynomials  in  one-dlmenslonal  fitting 
which  gives  the  "best"  (most  economical)  fit  for  a polynomial 
of  given  order. 

Kr.inkly  experimental  methods  are  used  to  get  an  acceptable  fit 
for  the  curve  f.tmilles  encountered.  As  a working  tool,  the  efficient 
Chebyshev  approx lm.it  Ion  progr.im  available  as  Program  14  in  the  Hl’-b7/ 
lll’-67  Stat  Pac  I is  exploited.  It  Is  also  advisable  to  stare  hard 
and  long  at  the  graphs  ot  the  particular  family  to  get  Ideas  from  the 
geometry  (one  version  of  the  "low  cunning"  approach  to  ad  hoc  computing). 
Generalization  to  vour  problems  of  the  methods  employed  In  the  examples 
below  cannot  be  guaranteed. 

._3_._  TRUE  MACH  ,m'_MBER 

Staring  at  the  upper  set  of  stralght-l Ine  segments  In  Fig.  21.1 
gives  the  Image  that  tliey  might  all  originate  as  a sheaf  from  a cvmi- 
mon  origin.  Use  of  a straightedge  shows  that  this  surmise  is  at  least 
approximately  true.  Hence,  try  as  the  functional  form  for  the  fit, 
the  relation 

* * 

A program  for  the  solution  of  b equations  in  b unknowns  is 
given  In  the  HP-b7/HP-97  Users  l.lbrary  Solutions  book  "High-level 
Math,"  programmed  by  R.  E.  HeBolt. 


- rtizmsr.. 
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X = m(H)  • G - 20  . 

where  m(H)  is  the  slope  for  altitude  H in  kft,  G is  the  averane  gross 
weight  in  klb,  and  x is  the  dummy  vaiiable  for  the  nomogram.  We  now 
build  the  following  table: 


:i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

H 

m(H) 

A 

C?  t»(H) 

A 

exp 

(quad) 

ao  = 1 

0 

1 . 000 

0.198 

0 

0.1807 

1 . 000 

1.000 

5 

1.198 

0.251 

■^.1807 

0.1902 

1.196 

1.196 

10 

1.449 

0.317 

0.3709 

0.1978 

1.442 

1.442 

15 

1.766 

0.348 

0.5687 

0.1799 

1.752 

1.753 

20 

2.114 

0.455 

0.7486 

0.1949 

2.146 

2.147 

25 

2.569 

0.772 

0.9435 

0.2628 

2.649 

2.649 

30 

3.341 

0.892 

1.2063 

0.2366 

3.295 

3.296 

35 

4.233 

1.086 

1.4429 

0.2220 

4.130 

4.131 

40 

5.309 

1.206 

1.6649 

0.2092 

5.219 

5.220 

45 

6.515 

1.8741 

6.645 

6.646 

Column  (2)  is  obtained  by  reading  differences  from  the  figure  and 
dividing.  Considerable  noise  can  be  introduced  by  tliis  process.  (1 
used  an  8X  loupe.  It  is  probably  better  to  plot  the  values  to  a large 
scale  and  fair  a curve  through  the  points,  and  then  read  off  values.) 
The  increasing  first  differences  of  column  (3)  suggest  that  an  ex- 
ponential form  be  used  since  the  exponent  will  be  much  flatter.  This 
is  shown  by  the  differences  of  column  (5). 

A quadratic  fit  to  m(H)  should  be  good.  We  find  from  the  HP~b7 
Stat  Pac  1 program  (six  minutes  to  run)  that 

CTiindi)  - -2.0A  X 10~^  + 0.0351H  + 1.56  x lo"^  . • (3) 


1 


A comparison  of  columns  (1)  and  (b)  which  Is  m(H)  sliows  the  adequacy 
of  the  fit.  Column  (7)  shows  that  the  constant  term  may  safely  be 
put  equal  to  0. 

Tlie  lower  family  of  curves  In  Fig.  21.1  requires  a true  two- 
dimensional  fit.  I’his  family  Is  well  behaved  In  that  a second-order 
pi'lynomlal  fit  looks  promising  (Kq.  (1)).  Moreover,  try  the  fit  with 
the  coefficient  a^  of  the  cross  term  xy  put  equal  to  0,  because  the 
curves  are  so  nearly  parallel.  That  Is,  9f/3x  = a,  + 2a.,x  + a,v  and 
the  dependence  on  y Is  weak. 

P r oc e ed  as  follows. 

1.  Re.td  from  the  curves  the  values  of  M at  33  points,  using 
Increments  of  20  for  x and  30  for  y (the  drag  Index).  Do  not  use 
y • 300,  reserving  It  for  an  extrapolation  check. 

2.  Use  the  HP  Stat  Pac  1 Program  lA  to  get  the  following  direct 
and  cross- fits.  (The  second-order  Chebyshev  fit  Is  used  and  the  time 
per  case  Is  well  under  10  min.) 


f(x.O) 

= 

0.330 

+ 

0.00300X 

- 

1.8133 

10~^ 

f(x.30) 

= 

O.A77 

+ 

0.00327X 

- 

2.034 

\ 

lo”'’ 

f (x, 100) 

0.430 

+ 

0.00539X 

- 

1.987 

\ 

lo"^ 

f (X.130) 

= 

0.398 

+ 

0.0037  lx 

- 

2.232 

\ 

10-'' 

f (x.200) 

= 

0.378 

+ 

0.00343X 

- 

1.87  3 

\ 

u)-^ 

f (O.y) 

= 

0.331 

- 

O.OOl22y 

+ 

2 . 286 

X 

10-^ 

f(20.y) 

S 

0.611 

- 

0. 00091 y 

1.143 

X 

10-^ 

f(AO.y) 

- 

0.690 

- 

0.00083y 

+ 

1.143 

X 

10-" 

f (60, y) 

- 

0.771 

- 

0.00l02y 

+ 

1.837 

X 

10-^ 

f(80.y) 

= 

0.809 

- 

0 . 00090y 

+ 

1.371 

X 

10"*’ 

3.  The  near  constancy  of  the  coefficients  Is  promising.  Using 
simply  their  means. 


f(x,y)  - f(0,y)  + 0.0033bx  - 1.993  x 10~^  x" 
f(x,y)  - f(x,0)  - 0,00098y  + 1.600  x lo"^  y’  . 


From  these 


t(x,y)  - f(0,0)  + 0.00536X  - 1.993  10~^  x^ 

- 0.00098y  + 1.600  x lo“^  . 


where  f(0,0)  « 0.53. 

A.  Programming  the  last  expression  for  f(x,y)  and  checking  its 
output  against  the  35  observed  values  yields  a good  fit.  But  using 
the  HP  improves  one's  "nose  for  numbers."  There  are  some  systematic 
biases  in  the  fit.  The  dependence  on  x is  somewhat  strong,  and  t lie 
dependence  on  y can  be  weakened  sliglitly.  Adjusting  to 

f(x.v)  - 0.53  + 0.0052X  - 2 x lo"^  x‘  - O.OOly  + 1.5  x io“^  y'  , (A) 


the  mean  absolute  error  with  respect  to  the  35  observations  is  0.0076 
and  the  maximum  error  is  0.016. 

5.  It  is  now  trivial  to  program  Kqs.  (3)  and  (A).  Tlie  eight 
constants  can  be  stored  in  primary  registers  0 through  7 and  recorded 
via  f W/DATA  on  side  two  of  the  program  card.  Running  time  is  tliroe 
seconds.  The  mitput  over  tlie  entire  nomogram  agrees  to  ±0.01. 
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22.  TKN-POINT  GAUSSIAN  INTEGRATION 


22.1.  REFKRKNCK 

a.  M.  Abramowltz  ainl  1.  A.  Stegun  (eds.),  !Ljrui}\'ok  of  Mathe- 
Fimo  t i 0*1^ , National  Bureau  of  Standards  Applied 
Mathematics  Series  55,  O.S.  Department  of  Commerce,  3d 
Printing,  March  1965. 

22.2.  DISCUSSION 

This  section  gives  a utility  program  to  evaluate  definite  Inte- 
grals with  high  accuracy.  If  the  Integrand  becomes  Infinite  at  some 
point  within  the  limits  of  Integration,  divide  the  Integral  Into  two 
parts,  using  as  limits  values  slightly  less  and  greater  than  that  point. 
Accuracy  is  checked  by  varying  these  values  and  reevaluating.  In  fact, 
even  If  the  Integrand  does  not  exhibit  this  behavior,  accuracy  can  be 
checked  by  dividing  the  Interval  Into  two  or  more  parts. 

llauss's  formula  for  an  arbitrary  Interval  Is  (.Ref.  a,  P-  S87, 
25.4.30) : 


The  abscissae  ate  the  zeros  of  the  Legendre  orthogonal  polynomial 
P^lx).  The  weights  w^^  are  given  by  a formula  Involving  P^(x). 

The  values  are" 


X 

o 

’‘l 

X, 

X3 

’‘a 


0.148874339 

w 

0 

0.295524225 

(^433395394 

w 

1 

m 

0.269266719 

0.679409568 

W2 

m 

0.219086363 

0^6  5063367 

W3 

s 

0.  149451349 

^73906529 

'"'4 

■1 

0.066671344 
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I'liL'so  are  keyed  into  a data  eard  will)  x.  stored  in  primary  re>;isters 

0 to  4,  and  w.  in  tlie  secondary  registers. 

1 

22.4.  _ 

The  program  is  sLrai^;l)ltorward  and  has  no  features  of  interest. 

I’ROBI.KN 

Write  a program  to  evaluate  elliptic  integrals  of  tl)e  first  kind: 


K(iJ).a) 


2 . 2 -1/2 
sin’’a"sin  0)  “d0. 


(.1)  Load  data  and  program  cards.  The  integral  will  be  evaluated 

in  the  r.adian  fh  IwM))  mode,  but  4)  and  a are  usually  given  in  degrees. 

(2)  STO  0 in  A and  4"  i't  degrees  it)  B. 

•STO  ot  in  degrees  in  K. 

(3)  CTO  A and  switch  to  W/PKGN.  Key  in  the  steps: 

•) 

h R^M).  RCL  B,  g -►  R,  STO  B.  RCL  K.  g R,  f sin,  g x"  STO  K. 

(4)  Switch  to  RUN,  OTO  L,  switch  to  W/l’RCM,  and  key  in  the 

steps : 

f sin,  g x',  RCL  K,  x,  CHS,  1,  +,  f/x,  h 1/x. 

(5)  Switch  to  RUN  and  press  A.  HSP  8. 

4)  = 30°,  a = 40°,  F(30,40)  = 0.533  427  45 
4)  = 85°,  a = 60°,  F(65,60)  = 1.348  H26  43. 

These  agree  exactly  with  the  tabular  entries  (Table  17.5)  of  Ref.  a. 
Running  time  is  about  30  sec. 


t. 
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22.6  TEN  POINT  GAUSSIAN  INTEGRATION 


kEVeNTSV  KtVCOOF 


COMMENTS 


STEP  KEY  ENTOV  KEY  CODE 


INDIRECT  ADDRESSING 


DEFINE  f(y) 


05r  *LBiE 


IS  STORED 


iM-!SS,'433« 
C.  -(33395394 
O.t '9409568 
C. '6506336? 
0.9 '3906529 
0.000000000 
0.000000000 
0. 000000000 
0.000000000 
0. 000000000 
0.000000000 
0.000000000 
0. 000000000 
0.000000000 
0.000000000 
0.000000000 


0.295524225 

i.269266«'19 

0.219086363 

0.149451349 

2.066601344 

0.000800000 


PARTIAL  SUM 


I « 4 ? 
LOOP 


INTEGRAL 
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23.  TRUTH  TABLKS 


23.1,  REKERENCES 

a.  John  E.  Pfeiffer,  "Symbolic  Louie,"  Soientifia  f\mcT*iaan, 

December  19!>0,  pp.  22-24. 

b.  E.  C.  Berkeley,  (/I’.znf  Bminti,  John  Wiley  & Sons,  New  York, 

1949. 

c.  Walter  E.  Cnshen,  "Symbolic  LorIc  In  Operations  Research," 

In  OiH'vatio*uf  h'cdt.\u\'h  fot'  , J.  K.  McCloskey  and 

F.  N.  Trefethen  (eds.).  The  Johns  Hopkins  Press,  Baltimore, 

1954. 

d.  R.  M.  Smldt  and  1.  L.  Reis,  "Symbolic  Louie  and  Plant  Loca- 
tion," of  luLiut'tJ'ial  Bnjtnoori^ui , Vol.  14,  No.  1, 

January-February  1963,  pp.  18-21. 

e.  P.  F.  Strawson,  rntvoduotion  to  U^ijioal  ThooJ'it,  Methuen  & 

Co.  Ltd.,  London,  John  Wiley  & Sons,  New  York,  1952. 

• f.  H.  M.  Semarue,  "Symbolic  Logic  In  Language  Engineering," 
of  tiu'  Wotdct'n  doiut  Corji^ttor  Conforonoe,  May 
3-5,  1960‘,  pp.  61-71. 

23.2.  DISCUSSION 

Tills  section  describes  a calculus  of  propositions  based  on  a binary 
algebra  for  logical  connections  and  operations  particularly  suited 
for  calculator  Implementation  as  well  as  easy  algebraic  manipulation. 

Letters  of  the  alphabet  will  stand  for  propositions.  To  illus- 
trate, choosing  propositions  to  be  used  In  an  example  to  follow,  let 

a stand  for  "A  ituin  Is  a mathematician." 
b stand  for  "A  man  likes  whisky  at  night." 
c stand  for  "A  man  likes  Mozart  In  the  morning." 
d stand  for  "A  man  waits  20  minutes  for  a bus." 

In  the  two-valued  calculus,  a proposition  has  a truth  value  of  0 (false) 
or  I (true),  a - 0 or  a - 1.  For  example,  writing  a - 1 means  "A  man 
Is  a mathematician." 

Propositions  may  be  operated  on  or  connected  by  logical  operations. 

A truth  table  corresponds  to  each  such  operation.  Ordinary  language 

F 
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will  be  used  instead  of  special  symbols.  For  each  basic  logical  oper- 
ation, a truth  table  shows  the  truth  value  of  the  operation  for  all 
combinations  of  truth  values  of  its  component  propositions.  Such 
basic  truth  tables  can  be  summarized  in  binary  algebraic  form.  This 
is  verified  in  the  right-hand  column  of  the  tabulation  below. 


Operation 


Truth  Table 


Binary  Algebraic  Form 


a AND  b 


a AND/OR  b 


a OR  ELSE  b 


IF  a THEN  b 


b 

a AND  b 


a AND/OR  b 


0 10  1 
0 0 11 
0 0 0 1 


0 0 


0 1 
1 1 
1 1 


0 10  1 
0 0 11 


a OR  ELSE  b 0 I 


IF  a THEN  b 


0 10  1 
0 0 11 
10  11 


a + b - ab 


a + b - 2ab 


1 - a + ab 


(Note  that  a false  proposition  can  imply  either  truth  or  falseness. 
Also  observe  that  'NOT  (IF  a THEN  b) ' is  'a  AND  NOT  b'.) 


NOT  BOTH  a AND  b 


1 - ab 


0 0 11 


NOT  BOTH  a AND  b 1 


(This  is  the  negation  of  AND.) 


1 0 


NEITHER  a NOR  b 


b a 0 10  1 

b 0 0 11 

NEITHER  a NOR  b 1 0 0 0 


1 - a - b + ab 


(This  is  the  negation  of  AND/OR.) 


a LIKE  b 


a 0 10  1 

b 0 0 11 

a LIKE  b 1001 


1 - a - b + 2ab 


(This  is  equivalence  and  is  the  negation  of  OR  ELSE.) 
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Binary  alRobra  has  some  InterestliiR  features  not  found  In  ordinary 

alRebra.  Since  the  variables  can  take  on  only  tl»e.  values  0 and  1. 

•) 

a"  " a and  ab(l  - a + ab)  - ab  - ab  + ab  - ab,  which  Is  readily  chciked 
by  ttie  truth  table: 

a 0 1 0 ’ 

b 0 0 11 

ab  0 0 0 1 

1 - a + ab  1011 

ab(l  - a + ab)  0 0 0 1 . 

On  the  other  hand,  a proposition  appearitiR  on  botli  sides  of  an  - 
sign  cannot  bo  cancelled,  because  this  could  be  divisli’n  by  0. 

The  use  and  manipulation  of  these  binary  algebraic  functions  can 
be  Illustrated  by  a problem  that  Waiter  I’ltts  of  tlu*  Massachusetts 
Institute  of  Teclinology  set  long  ago  on  an  exiuiiinat  ion  (Ref.  a). 

Suppose  1 through  4 below  are  known  to  be  true: 

I.  Ij'  ii  luji'c  f.'  :'.0  rv:Kuti't^  /.')• 

>j  then  he  either  likee  Mer,ert  :n  the  mi ev  u'h:ek;i  et 

nujht,  hat  l\>th. 

Ij'  e ffiein  tikee  ;.'h:ek:i  at  >.’/,?/:(,  then  iu'  eithi'r  .'ikea 
Afi’r.iiJ*(  ni  tiu'  mevfityuj  atui  Jeer-  >u't  luu'e  te  '.'0  m fey  a 
hue  oy  Iu;  doee  ni't  like  Me:XiV't  in  tfte  neynhiii  kIiu!  hie  te  u\i:t 
mintitee  fciy  a hue  ey  ehn'  lu'  ie  >u'  rhtthi'nat  feian. 

If  a nan  likt'e  MeHiiyt  in  the  nnmi'n.j  atui  yleee  >u't  lujy,- 
(o  lAiit  "(1  mnintra  foy  a buy,  tlwn  he  likee  idiieku  at  nt\i>it  . 

■!.  If  a mthematieian  likee  Mer.ayt  t>i  the  rknnihuj,  )u’ 
either'  likee  uhieki/  at  ni.jht  ey  fuie  te  u\iit  ninuti-y  fi'y  a 

bun;  eenyeyyelif,  if  he  likey  lyhiykit  iit  ni^jht  a>ki  luiy  fo  uw't  b.O 
mnu«t»’a  fey  a Imy,  >te  iy  a mt  h>'rkitiei<in — if  he  likey  S1e:k;yt  in 
tihi  neynituj. 

Then: 

When  deey  a natiumitieian  u\iit  bO  ninutey  fey  a buy'.' 
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To  solve  this  problem,  first  translate  each  condition  from 

* 

English  to  the  language  of  propositions  and  logical  connectives, 
and  then  express  the  result  in  binary  algebra.  Because  each  of  the 
four  above  conditions  is  true,  each  algebraic  expression  is  put 
equal  to  1 and  then  simplified. 

1.  IF(a  AND  NOT  d)  THEN(b  OR  ELSE  c). 

1 - a(l  - d)  + a(l  - d)(b  + c - 2bc)  - 1 , 

a(l  - d)(b  + c - 2bc  - 1)  = 0 . (1) 

2.  IF  b THEN(((c  AND  NOT  d)  OR  ELSE(NOT  c AND  d)) 

OR  ELSE (NOT  a)). 

To  simplify,  set  A = (c  AND  NOT  d)  OR  ELSE(NOT  c AND  d) . 
Then  A = c(l  - d)  + d(l  - c)  - 2cd(l  - c)(l  - d) 

= c + d - 2cd  . 

That  is,  the  proposition  A is  equivalent  to  the  proposition 
"c  OR  ELSE  d".  Making  this  substitution,  we  obtain 

l-b+b(c+d-  2cd  + 1 - a 

-2(1  - a)(c  + d - 2cd)  = 1, 

b(-a  + (c  + d - 2cd)(2a  - 1))  = 0 . (2) 

3.  IF(c  AND  NOT  d)  THEN  b. 

1 - c(l  - d)  + cb(l  - d)  = 1 , 

c(l  - d)(l  - b)  = 0 . (3) 

4a.  IF(a  AND  c)  THEN(b  OR  ELSE  d). 

1 - ac  + ac(b  + d - 2bd)  » 1 , 

ac(b  + d - 2bd  - 1)  » 0 . (4) 


Ref.  e is  helpful  in  this  respect, 
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4b.  IF  c THEN(IF  (b  AND  d)  THEN  a). 

l-c+c(l-bd+  abd)  - 1 , 

bcd(l  - a)  - 0 . (5) 


To  answer  the  question,  put  a * 1 and  d « 1 in  (1)  through  (5), 
since  these  two  propositions  must  be  true.  The  interpretation  of  the 
question  is  delicate.  The  question  is  rephrased  here  as:  What 
values  of  b and  c are  associated  with  a = 1 and  d = 1 to  make  all 
conditions  of  the  problem  true?  The  set  of  conditions  reduces  to 

be  = 0 . 

This  means  NOT  BOTH  b AND  c (that  is,  1 - be  = 1),  which  may  be  expressed 

A.  When  he  likes  neither  Mozart  in  the  morning  nor  whisky  at  night. 

B.  When  he  likes  whisky  at  night  and  not  Mozart  in  the  morning. 

it 

C.  When  he  likes  Mozart  in  the  morning  and  not  whisky  at  night. 

Of  course,  many  questions  other  than  the  one  above  could  be  asked. 

To  be  exhaustive,  this  means:  Find  all  values  for  the  set  of  proposi- 
tions (a,b,c,d)  that  satisfy  all  given  conditions.  The  program  of  this 
section  is  designed  to  examine  systematically  all  cases,  here  16  in 
number.  Programming  the  conditions  and  running  the  program  yields  only 
eight  satisfactory  combinations  of  (a,b,c,d). 

£ b £ ^ 

0 0 0 0 

0 10  0 
0 0 10 

0 0 0 1 

10  0 1 

110  1 

0 0 11 

* 10  11. 


Something  of  a furor  erupted  in  the  Letters  section  of  the  Sci- 
entific: American  of  February  1951  because  Pfeiffer  had  not  given  an 
answer  in  his  article.  Pfeiffer  and  Morris  gave  the  first  part  of  C 
as  the  answer.  Krause  gave  the  first  part  of  B,  and  Bomgren  gave 
answer  A. 
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This  is  tlie  truth  table  method.  The  lines  indicated  by  arrows 
are  ttie  answers  to  tiie  given  problem. 

The  logical  structure  of  the  problem  may  also  be  captured  by 
using  the  truth  table  to  construct  a Venn  diagram. 

Suppose  a point  placed  on  this  page  represents  a man.  Group 
together  all  men  who  are  mathematicians  so  tliat  this  subset  is  enclosed: 


then  all  inside  men  have  the  value  a = 1,  those  outside  tlie  value 
a = 0. 

From  the  truth  table,  whenever  a = 1,  d = 1.  But  there  are  cases 
where  a = 0 and  d = 1.  Hence  the  set  a is  properly  contained  in  the 
set  d.  Logically  this  means  "if  a then  d",  and  algebraically  "a  = ad". 
Similarly,  c is  contained  in  d.  But  for  a and  c,  by  the  last  four 
lines  of  the  truth  table,  there  are  the  combinations  10,  01,  11.  Hence 
a and  c have  points  in  common — they  intersect. 

Turning  to  the  set  b,  since  be  = 0 in  all  cases,  the  sets  b and  c 
are  disjoint — they  do  not  intersect.  But  b intersects  a and  also 
intersects  d.  However,  because  of  the  line  1101,  and  because  there  is 
no  line  0101,  these  latter  two  intersections  are  the  same  subset. 

Putting  all  of  tliis  together  yields  the  following  diagram: 


i 
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Returnlng  to  the  original  algebraic  formulation,  it  is  readily 
seen  that  all  conditions  are  satisfied  by  appeal  to  this  diagriun.  For 
example,  since  be  = 0 and  abd  = ab,  in  (2)  above  ab  = bd . 

The  original  question  asked  by  this  problem  is  not  clear.  It 
would  be  better  to  ask,  W}iat  can  be  said  about  a vnthematiaian? 

The  answer  to  this  question  is: 

• He  always  waits  20  minutes  for  a bus; 

• He  does  not  like  both  whisky  at  night  and  Mozart  in  the 
morning,  although  he  may  like  one  or  the  other. 

The  calculus  of  propositions,  a branch  of  symbolic  logic,  has 
found  applications  in  optimizing  switching  circuit  design,  in  deter- 
mining insurance  eligibility  (an  ex^unple  from  Ref.  b will  be  given  in 
the  next  section),  in  deciding  on  plant  location  (Ref.  d),  and  in  the 
interpretation  of  contracts  and  law.  Walter  Cushen,  in  a fascinating 
chapter  in  Ref.  c,  discusses  applications  to  production  engineering 
and  to  conflicts  formulated  as  multi-move  games. 

23.3.  EQUATIONS 

None . 

23.4.  PROGRAM  NOTES 

Since  the  program  flow  is  somewhat  complex,  a flowchart  is 
provided . 

Suppose  there  are  n propositions.  Then  there  are  2*^  cases  to 
examine,  which  are  actually  numbered  0,  1,  ...,  2*^-1.  Each  case 
number  in  its  turn  is  reduced  to  its  binary  form,  but  stored  backwards 
in  R,  to  R . Indirect  addressing  is  used. 

Each  condition  is  examined  in  turn,  the  examination  halting  at 
the  first  condition  that  is  not  satisfied.  If  all  conditions  are 
satisfied,  the  binary  case  number  is  displayed.  R/S  continues  to  the 
next  valid  case.  The  end  is  signalled  by  the  decimal  display  of  2*^  - 1. 

As  programmed,  there  is  space  for  7 propositions.  If  a problem 
requires  more  than  7,  make  some  slight  programming  changes.  Primary 
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registers  8 and  9 are  shifted  to  secondary  8 and  9.  Then  up  to  18 
propositions  can  in  principle  be  handled.  For  large  problems,  however, 
there  is  apt  to  be  a large  number  of  conditions  and  there  may  not  be 
adequate  space  available  to  program  them.  Moreover,  execution  time 
will  be  long.  It  is  wise  to  do  as  much  algebraic  manipulation  on 
the  set  of  conditions  as  is  feasible  to  simplify  the  set. 

In  problems  where  some  propositions  are  held  constant,  the  con- 
stant value (s)  (0  or  1)  are  stored  manually  in  some  register (s)  above 
R^,  but  only  if  these  are  needed  in  programming  the  conditions. 


EXAMPLE 

This  is  a group  insurance  problem  taken  from  Ref.  b (pp.  161-165). 
The  rules  (conditions)  applying  to  employees  are: 

1.  Any  employee,  to  be  insured,  must  be  eligible  for  insurance, 
must  make  application  for  insurance,  and  must  have  such 
application  for  Insurance  approved. 

2.  Only  eligible  employees  may  apply  for  insurance. 

3.  The  application  of  any  person  eligible  for  Insurance  without 
medical  examination  is  automatically  approved. 

4.  (Naturally)  an  application  can  be  approved  only  if  the 
application  is  made. 

5.  (Naturally)  a medical  examination  will  not  be  required  from 
any  person  not  eligible  for  insurance. 

The  propositions  are  5 questions  about  an  employee  to  be  answered 
"yes"  (1)  or  "no"  (0).  These  are: 

a:  Is  the  employee  eligible  for  insurance? 

b:  Has  the  employee  applied  for  insurance? 

c:  Has  the  employee's  application  for  insurance  been  approved? 

d:  Does  the  employee  require  a medical  examination  for  insurance? 

e:  Is  the  employee  insured? 
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The  conditions  are  translated: 


1.  IF  e THEN(a  AND  b AND  c) 

e(l  - abc)  “ 0 . 


2.  IF  b THEN  a 


b(l  - a)  - 0 . 


3.  IF  a AND  b AND  NOT  d THEN  c 

ab(I  - d)(l  - c)  - 0 . 


4.  IF  c THEN  b 

cll  - b)  - 0 . 

5.  IF  NOT  a THEN  NOT  d 

d(l  - a)  = 0 . 

The  question  is  "What  are  the  possible  statuses  of  employees  who 
are  not  Insured?"  This  means  that  e must  be  put  equal  to  0.  But  then 
the  first  condition  is  irrelevant  since  a false  proposition  Implies 
any  proposition. 

SOLUTION 

Load  program.  Key  GTO  B.  Switch  to  W/PRGM.  Now  key  in  the  con- 
dition in  the  order  2,  4,  5,  3.  That  is,  the  simplest  conditions  are 
entered  first  to  save  execution  time.  The  program  steps  are: 
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(Nolo  that  If  conditions  J to  5 wore  multii>llod. 


abed (I  - a)(l  - b)(l  - c)(l  - d)  « 0 . 

which  is  true  for  all  cases.  The  mnl t ip  1 Icat ion  lias  destroyed  the 
meaning  of  the  Individual  conditions  by  absorption.) 

Now  switch  to  RbN.  Key  4 (the  niunber  of  propositions)  and  Press 
A.  See  0.  This  is  actually  the  status  0000  for  a,b,c,d.  On  succes- 
sive presses  of  R/S  you  will  see  1000,  1110,  1001,  1101,  1111,  15 
(the  number  of  cases  minus  1). 

Using  the  definitions  of  a,b,c,d,  these  5 statuses  rapidly  trans- 
late into  an  answer  to  the  question. 
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